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The Heat Capacity of Pure Silicon and Germanium 
and Properties of their Vibrational Frequency Spectra} 


By P. Fuusacuer}, A. J. Lrapserrert 
and J. A. Morrison 
Division of Pure Chemistry, National Research Council, Ottawa, Canada 


[Received September 30, 1958] 


ABSTRACT 


Heat capacities of pure silicon and germanium have been measured between 
2-5° and 300°K. The estimated accuracies of the measurements are +0-5% 
for 10°<7<20°K and +0-2% for 7’>20°K. The results for silicon were 
anomalous in the region 7’<7°K; it is suggested that this behaviour resulted 
from adsorption and desorption of exchange gas in the calorimeter vessel. 
No anomalies were observed in the measurements on germanium; their 
accuracy at the lowest temperatures is estimated to be +2%. The results for 
both substances are significantly different from previously published values. 

The temperature variation of @p of both silicon and germanium below 
about Z7’'= @p/3 is consistent with general harmonic theory, but at high 
temperatures ©p decreases with increasing temperature. This appears to be 
a clearly marked anharmonic effect of a type that cannot be explained by 
change of volume alone. The results for 7’ <@p/3 are therefore extrapolated 
to high temperatures to give a heat capacity consistent with harmonic theory 
at all temperatures. From this, properties of an effective harmonic spectrum 
are derived; these include the low frequency expansion and a number of 
positive and negative moments. The apparent anharmonic contribution 
AC,=[C,—C,y (harmonie)] is found for both elements to be approximately a 
linear function of the temperature in the region between ©p/3 and room 
temperature. The possible significance of this result is discussed in the light 
of anharmonic theory. 

The results are compared where possible with vibrational properties 
determined from other experiments and the agreement is very satisfactory. 


§ 1. INTRODUCTION 


Tue diamond structure elements (diamond, silicon, germanium and grey 
tin) present attractive features for an investigation of lattice vibrations. 
Their lattices have a simple crystal structure and their interatomic forces 
should not be very complex. It might be hoped therefore that appro- 
priate models for the crystals could be devised and tested adequately 
with different measured vibrational properties. A model for diamond 
based on short-range forces has been proposed (Smith 1948) but it does 
not seem to give thermodynamic properties in agreement with experiment 
(Desnoyers and Morrison 1958). The same model applied to silicon and 
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germanium (Hsieh 1954) shows an even wider disagreement. Recently 
Lax (1958) has suggested that longer range interactions of a special type 
may need to be considered. 

Although at present a satisfactory model has not been developed there 
is now the possibility of correlating all measurements of vibrational 
properties of silicon and germanium without specific reference to a model. 
Recent work on neutron scattering by germanium (Brockhouse and 
Iyengar 1957, 1958, Pelah et al. 1957) has yielded among other things 
frequency-wave number relations for lattice waves in certain directions in 
the crystal. These have been compared (Brockhouse and lyengar 1957, 
1958) with measured elastic properties (McSkimin 1953) and infra-red 
absorption (Collins and Fan 1954, MacFarlane et al. 1957) and the agree- 
ment is found to be excellent. As a further development Phillips (1959) 
has deduced a lattice frequency spectrum for germanium by a moment- 
singularity method using the dispersion relations found by neutron 
spectrometry. The way is therefore clear to include thermodynamic 
properties in the comparision of vibrational properties since, in the 
harmonic approximation, they depend only upon the frequency spectrum. 

Existing values of the heat capacity of silicon (Pearlman and Keesom 
1952) and of germanium (Esterman and Weertman 1952, Hill and Parkin- 
son 1952, Keesom and Pearlman 1953) cover only limited temperature 
regions and where the measurements have been duplicated the agreement is 
hardly satisfactory. Hence it has seemed valuable to obtain new results 
for as wide a temperature range as possible. In the present paper we report 
measurements for both substances in the temperature region 2-5° to 300°K. 
The specimens used were of such purity that electronic contributions to 
the heat capacity were entirely negligible below 300°K. 

Since it is the heat capacity at constant volume which is directly related 
to the lattice frequency spectrum, this quantity must be computed from 
the experimentally measured heat capacity at constant pressure. For 
both silicon and germanium the difference can be determined accurately 
because both the thermal expansion and the compressibility of the lattices 
are well known for a wide temperature range below room temperature. 
There is, however, an additional factor to be considered which introduces 
unavoidable uncertainty ; this concerns possible effects of anharmonicity 
in the lattices. In the present examples of silicon and germanium we 
shall see that @,, decreases with increasing temperature at higher tempera- 
tures. This behaviour cannot be accounted for by harmonic lattice 
theory. A heat capacity compatible at high temperatures with harmonic 
theory can be obtained by extrapolation of the results from lower tempera - 
tures (7'< @)/3); such a procedure was used previously in the analysis of 
heat capacity results for ionic crystals (Barron et al. 1957, 1959). 

The heat capacity is taken to correspond to the effective harmonic 
spectrum of the lattice and the good agreement of the derived properties 
with those obtained from other work suggests that this is a fruitful concept. 
The detailed interpretation of the anharmonic effect is not simple but a 
general discussion is given in § 4.2, 


i | 
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§ 2. EXPERIMENTAL 
2.1. Calorimetric Specimens 


The silicon (du Pont Hyperpure grade) was obtained in the form of 
eleven single crystal slabs whose average resistivity at room temperature 
was approximately 100ohmem. This corresponds to a carrier concen- 
tration of not more than 10%cem-*. The germanium (supplier: Société 
des Mines et Fonderies de Zinc de la Vieille Montagne, S.A., Belgium) 
was received in the form of rods and was stated to have a minimum 
resistivity of 5|00hmem. The resistivities of the particular rods used for 
the calorimetric work were determined by G. Fischer of the Division of 
Pure Physics and found to be 590hmecm at room temperature. Low 
temperature resistivity meaurements indicated a carrier concentration 
ofnot more than 10'*em~*. For both the silicon and germanium specimens 
the electronic contribution to the heat capacity was negligible down to 
the lowest temperatures. 

The specimens were broken into pieces of about 3mm in size for filling 
the calorimeter vessels. The filled vessels for both calorimeter assemblies 
used were evacuated to a pressure of 10-*mm of mercury and then sealed 
with a small amount of helium exchange gas inside. The amounts of the 
specimens and exchange gas are given in table 1. 


Table 1. Calorimeter Fillings 


Calorimeter I Calorimeter II 
Specimen Helium Specimen Helium 


(g atom) | (10-* g atom) | (g atom) | (10~-® g atom) 


Silicon 6:0318 2° 51658 4-8 
‘7 2-8729 4-7 


ew) 


. 


Germanium 2-8698 


b 


2.2. Calorimetry 

The higher temperature calorimeter (calorimeter II) was used exactly 
as in previous experiments (Berg and Morrison 1957) and no further 
description of it is necessary here. 

The essential features of the lower temperature calorimeter (calorimeter 
I) have been described previously (Morrison et al. 1955). We shall, 
however, discuss briefly the establishment of the temperature scale for 
the measurements with this calorimeter in the temperature range 2° to 
20°k. The working thermometer, a carbon resistor (Allen Bradley 10 ohms 
1 watt) attached to the calorimeter vessel, was calibrated for each set of 
experiments by comparing it against a normal hydrogen vapour pressure 
thermometer at five temperatures and a helium vapour pressure thermo- 
meter at at least one temperature. The same vapour pressure thermo- 
meter bulb, attached to the adiabatic shield, was used for both sets of 
measurements, The heating of the shield was so adjusted that the 


SZ 
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temperature of the calorimeter vessel did not change significantly during 
a comparison. Three of the fixed points (at about 4°, 14° and 20 a8) 
were used to determine constants of an interpolation formula (Clement 
and Quinnell 1952) which gave a working temperature scale; the rest 
of the fixed points served as checks. Liquid hydrogen temperatures 
were calculated from the vapour pressure equation of normal hydrogen 
(Woolley et al. 1948) and liquid helium temperatures from vapour pressure 
data for the 7’, scale (van Dijk and Durieux 1955). 

It remained to relate the working-temperature scale to the thermo- 
dynamic scale; this was accomplished by attaching a gas thermometer 
bulb to the calorimeter vessel and by comparing the gas thermometer with 
the carbon thermometer at about twenty points between 2-5° and 20°K. 
Below 2:5°K the gas thermometer was turned into a vapour pressure 
thermometer. At the same time the two thermometers were compared 
where possible with the vapour pressure thermometer on the shield. 
In this way small temperature differences due to thermal gradients and to 
vibrational energy appearing in the calorimeter vessel were automatically 
taken into account. Although the resistance of the carbon thermometer 
at a given temperature generally changed each time the thermometer 
was cooled from room temperature, the sensitivity did not alter signifi- 
cantly. In principle, therefore, one series of gas thermometer comparisons 
would have sufficed. Actually three sets of comparisons were made 
during the course of the heat capacity measurements. Their consistency 
suggests that the final temperature scale used agreed with the thermo- 
dynamic scale to within +0-02° from 2° to 20°K. The slope of the scale 
is estimated to be correct within +1% at the lowest temperature and 
H 0-0, LOPLI val Ol Ks 

§ 3. RESULTS 


3.1. The Heat Capacity Measurements 


The primary experimental results (the heat capacity at constant 
pressure) are tabulated in the Appendix. The temperature intervals 
used in the measurements were 0-5° (2°<7'<7°K), 1° (7°< 7 < 20°K) and 
3° (7'>20°K). Curvature corrections were applied where necessary. 
The thermochemical calorie (4:1840 abs. joule) is used throughout. 
Results for silicon for 7’<7-7°K are not given in the table because in this 
region anomalous effects appeared which can be attributed to adsorption of 
exchange gas in the calorimeter vessel. The details are discussed in § 3.3.2. 

The internal consistency of the results may be judged from figs. 1 and 2 
where the apparent Debye characteristic temperatures (©) coriesponding 
to the experimental heat capacities of silicon and germanium are plotted 
as a function of temperature. The sensitivity of the graphs are indicated. 

The values of ©, were obtained from C,, which was deduced from Cy 
by means of the usual thermodynamic formula 


aT V 


Cp— Cy= (3.1) 


Heat Capacity of Pure Silicon and Germanium 277 


where «a is the cubical coefficient of thermal expansion, V the molar volume 
and x the isothermal compressibility. Thermal expansion data obtained 
by Gibbons. (1958) for silicon and germanium were used; their estimated 
accuracy is 1%. The compressibilities were calculated from measured 
elastic constants (McSkimin 1953) and are probably reliable to at least 
3%. Hence the values of C,,—C, would seem to be accurate to 5%. 


150 200 250 300 


600 


x 
° 
a 
-) 
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1% t of Cy 
site 50 100 150 
T (°K) 


The apparent Debye characteristic temperature of silicon as a function of the 
absolute temperature. 


e Calorimeter I. 
© Calorimeter IT (12° < 7’<170°K). 
C] Calorimeter IL (170°K < 7’ <300°K). 
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Fig. 2 
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The apparent Debye characteristic temperature of germanium as a function of 
the absolute temperature. 
e@ Calorimeter [. 
© Calorimeter IT (12°< 7 <170°K). 
[] Calorimeter IT (170° < 7 <300°x). 
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Smoothed values of C!, and of C),—C,, for both silicon and germanium are 
given in table 2. Since the actual values of C,,—C', are so small, their 
uncertainty does not affect C, significantly. 


Table 2. C,, and C,,—C,, at Rounded Values of the Temperatures 
| Units — cal/g atom deg | 


see é 
Silicon Germanium 


Cy é C, Cy—-C 


v0 


(0-0000271) 0-0001413 
(0-000219) 0-001199 
0-001845 0-01394. 
0-007300 0-07672 
0-02265 0-2174 
0-05702 0-4140 
0-1149 0-6316 
0-1962 0-8531 
02957 1-067 
0-5272 1-479 
0-7725 1-879 
1-019 2-269 
1-261 2-640 
1-502 2-989 
1-739 3-302 
2-205 3-838 
2-650 4+261 
3-057 4-587 
3-420 4-837 
3-735 5-033 
4-009 5-190 
4-246 00: 5-318 
4-455 00: 5-423 
4-638 . od11 
4-796 5-590 


The accuracy of the measurements is therefore determined primarily 
by uncertainties in the temperature scales. The estimates for these for 
calorimeter I suggest that the heat capacity results are accurate to + 2% 
at the lowest temperatures and +0-5% in the region 10°<7'< 20°K. 
Previous work with calorimeter II (Berg and Morrison 1957, Morrison and 
Patterson 1956) indicates its accuracy to be +0-5% for 12°<T < 20°K 
and +0-2% for 7’>20°x. The results for the two calorimeters agree to 
better than +0:5°% between 12° and 20°K. 


3.2. Comparison with Other Measurements 


There are several previous reports of heat capacity measurements on 
silicon and germanium for more limited temperature ranges than have 
been covered in the present work. Except for low temperature results 
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A comparison of smoothed @p values for silicon. 
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for germanium (Keesom and Pearlman 1953) only smoothed values of the 
heat capacity are given. We have therefore chosen to compare the 
available results through difference graphs using as references smoothed — 
©,-T relations which reproduce the present results to well within their 
estimated accuracies. A comparison with the results of Pearlman and 
Keesom (1952) for silicon is shown in fig. 3, and it is evident that the 
deviations are large. It should be pointed out, however, that the speci- 
mens used by Pearlman and Keesom contained appreciable amounts of 
electronically active impurities. 

The results for germanium are shown in fig. 4. Here, the values obtained 
by Hill and Parkinson (1952) and by Esterman and Weertman (1952) 
deviate for the most part in opposite directions. Over the major portion 
of the temperature region covered the agreement between the present 
results and those of Esterman and Weertman is close to what one might 
expect for totally independent calorimetric measurements. The results of 
Keesom and Pearlman (1953) for the temperature range 1° to 4°K were 
fitted to the 7° law. The present results, however, show that ©, is not 
constant in this region. 

Additional heat capacity results for germanium in the temperature 
range 12 to 270°K have recently been obtained by Clusius and Piesbergen 
(1958). Between 20° and 270°K they agree with the present values 
within about 1%. Below 20°K the difference increases, becoming about 
10% at 12°K. 


3.3. Analysis of the Experimental Results 
3.3.1. General remarks 


The type of analysis to be outlined in the next section has been fully 
described previously (Barron et al. 1957, 1959). Its object is to deduce 
properties of the vibrational frequency spectrum from measured heat 
capacities. An essential factor is the recognition of and correction for 
anharmonic effects. 

In approximate general shape the ©, curves for silicon and germanium 
(figs. 1 and 2) are of the form of ©, curves found for potassium halides 
(Berg and Morrison 1957). Starting from the lowest temperatures ©) 
first decreases to a minimum, then increases to a flat maximum and finally 
decreases slowly. The decrease at high temperatures cannot be accounted 
for by harmonic lattice theory which requires that ©, should increase 
asymptotically to a constant value. We propose to adopt the view taken 
previously that the decrease in @, is to be associated with an anharmonic 
effect of a type which cannot be explained by change of volume alone ; 
we shall further assume that below 7’ = @,/3 the only anharmonic effect due 
to thermal energy which needs to be considered is thermal expansion. 


+ We should like to thank Professor Clusius for sending us these results. 
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3.3.2. Properties of the vibrational spectra 

The extrapolation of the results for intermediate temperatures to 
infinite temperature is illustrated in fig. 5. The solid curves represent 
the fitting of the results to the Thirring expansion which has the form 


ey'=0,!| 1-4 ($2) +B(52) ... |. oe Oe) 


The values of ©, used here have all been reduced to a fixed volume, namely, 
the volume of the crystal at 0°k. The intercepts of the solid curves at 
1/T2=0 give values of ©,, and hence values of the second moments, 1». 
The dashed curves represent estimates of limits on the extrapolations. 
The most probable values of ©, are 674+4°K for silicon and 395 + ak 


for germanium. 
Fig. 5 


4 
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0 5 10 ; (5 20 
lot 2 
The estimation of ©,,. 
©. Silicon. 
e Germanium. 


The coefficients A and B in eqn. (3.2) are related to 4 and ju, (see Barron 
et al. 1957); their values have been estimated graphically, but the 
details are not given here. 

The results obtained for very low temperatures have been used to 
determine the coefficients of the low frequency expansions of the spectra 
of silicon and germanium : 


G(v)=NV (av? + Byttyvet...). 0.0 fn. 0... (3.¢ 
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Fig. 6 
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The corresponding expansion for the heat capacity is 
O;o]* +01 Soll: Pe nb eo) 


the determination of the coefficients a, b and c is illustrated in figs. 6 and 7. 
The fitting of experimental results to eqn. (3.4) has been discussed 
elsewhere (Barron and Morrison 1957). All of the experimental results 
for silicon are included in fig. 6 and those for 7'<7°K show an unexpected 
behaviour. In the experiments with silicon between 2° and 4°K it was 
observed that the equilibrium time after the heating of the calorimeter 
vessel was very long, which indicated that the helium exchange gas was 
all adsorbed because normally the observed equilibrium time was a few 
seconds. If thermal equilibrium was not being attained, the apparent 
heat capactiy would be too small as was found. Desorption of helium 
would be expected to produce too large apparent heat capacities over a 
limited temperature range and it is significant that the observed heat 
capacities lead to a maximum in C/T? at about 6°K (fig. 6). No anomalies 
were observed in the experiments with germanium (fig. 7) where the surface 
area of the specimen was much smaller (less than half) and the total heat 
capacity of the system was some six times larger. It therefore seems 
probable that the anomalous results for silicon are to be accounted for by 
effects of adsorption and desorption of the exchange gas. Hence, it is 
justifiable to represent the true heat capacity of silicon for the region 
T'<7°K by the extrapolated portion of the solid curve shown in fig. 6. 


Table 3. The Coefficients in Eqns. (3.3) and (3.4), ©, and @(elastic) 


Silicon Germanium 

a (cal g atom-! ie 1-73 + 0-04 x 10-6 8-90 + 0-16 x 10-8 

b (cal g atom—! deg) 5+5~x 10-10 2-0 + 0-7 x 10-8 

¢ (cal g atom— deg-*) 6-0+1:4x10- 3-0 + 0-5 x 10-10 

a (sec?) 3°71 + 0-08 x 10-89 | 1-91 +0-03 x 10-38 

B (sec?) 9+9 x 10-65 3-6 + 1-2 x 10-83 

y (sec?) 4-4 + 1-0 x 10-89 2:2 + 0-4 x 10-87 
©, (previous work) 658° + 362 + 6° Bf 
®,) (present work) 645 + 5° 3714+ 2° 

@O(elastic) § 647-8° 374-0° 


if Pearlman and Keesom (1952). 
t Keesom and Pearlman (1953). 
§ de Launay (1956). 


The values of the coefficients occurring in eqns. (3.3) and (3.4) are given 
in table 3. Also shown are values of @, based on the present and on earlier 
results and @(elastic) as calculated by de Launay (1956). For both 
silicon and germanium ©, from the present heat capacity measurement 
and @(elastic) agree within the probable errors assigned to @,. 


Heat Capacity of Pure Silicon and Germanium 285 


Other moments of the frequency spectra, in addition to the even positive 
moments derived from the Thirring expansion, have been calculated 
from expansions involving integrals of the form 


"Co ap 
Jn 


(for details, see Barron et al. 1957). Allofthe computed moments, together 
with their estimated accuracies, are given in table 4. They are also 
illustrated by graphs of the function 


in fig. 8. This function would be a constant for a Debye spectrum; its 
primary usefulness is to serve as a check for internal consistency in the 
calculations. 


Table 4. Derived Moments of the Frequency Spectra 
Silicon Germanium 
pea. (8ee*!?) 1-355 + 0-004 x 10-32 | 5:73 +0-03 x 10-? 
jt» (sec?) 2-450 + 0-006 x 10-26 | 7-82 +0-02 x 10-26 
1.5 (sec?/?) 6-361 0-611 10-7° 1-277 + 0-003 x 10-19 
u_, (sec) 1-304 + 0-002 x 10-13 | 2-314 + 0-004 x 10-13 
y (sec) 8-93 +0-02 x10¥ 5-110 + 0-015 x 10" 
jt, (sec *) 1-003 + 0-004 x 1018 5:82 +0:02 x10” 
ju, (sec-2) 1-182 +0-014x 1028 | 4-06 +0-06 x 1025 
pug (Sec) 2:03 40-13" x 10° asi Danse Ue esoeces RUS 
fig (BOCT*) 30 0-0 10" 17554 0-21 5x 10% 
Table 5. The Entropy and Zero Point Energy 
S° at 298-15°K Zero point energy 
(cal/g atom deg) (cal/g atom) 
Silicon 4-497 + 0-009 1435 +5 


| Germanium 7-432 +0-015 833 +3 


The moment 1, is related to the internal energy of the lattice and vg, the 
geometric mean frequency, to the entropy. 4; gives the zero point energy 
directly (Barron et al. 1957) and numerical values of this quantity for 
silicon and germanium are given in table 5. Also shown are values of the 
entropies at 298-15°K, which have been calculated from Cy. 
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The heights of the vertical bars indicate the estimated uncertainties. Where 
no bar is shown, the uncertainty is less than the diameter of the plotted 
circles. 

$ 4. Discusston 
4.1. Correlation of Vibrational Properties 
It is appropriate now to see if the harmonic vibrational properties 
deduced from the heat capacity results are consistent with the results of 
other measurements. A difficulty here is that neutron scattering and 
some optical absorption measurements have only been made at room 
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temperature where anharmonic effects appear in the heat capacity. Care 
must be taken, therefore, in making direct comparisons. On the other 
hand, measurements of elastic properties (McSkimin 1953) have been 
carried down to low temperatures and a direct comparison of elastic and 
thermal properties is possible through @(elastic) and @,. As we have 
seen (table 3), the agreement between these two quantities is excellent for 
both silicon and germanium. 

Brockhouse and lyengar (1958) have compared their frequency-wave 
number relations for germanium with elastic properties corresponding to 
room temperature and with infra-red absorption results and have shown 
that all of the data are mutually consistent. They have also included 
existing heat capacity results in the comparisons but only to a limited 
extent and without specific consideration of anharmonicity. Concurrently 
with the present work, Phillips (1959) has determined a vibrational 
frequency spectrum for germanium and this allows a better comparison 
to be made with the heat capacity. Phillips used the experimental data 
of Brockhouse and lyengar to determine critical points occurring in the 
spectrum and constructed the remaining parts of the spectrum by his 
interpolation method (Phillips 1956). The resulting spectrum gave 
©, values some 3% larger than those calculated from the experimental 
heat capacities for the temperature range 15° to 125°K (Phillips 1959). 
By reducing all of the neutron frequencies by an amount within their 
probable error, Phillips has constructed a second spectrum for which the 
©,, values are within 1° of those derived from our experiments. 

It is satisfying that Phillips’s spectrum, which has been derived assuming 
harmonic vibrations, fits the heat capacity results rather well in the 
temperature region where we have supposed the behaviour of the lattice 
vibrations to be harmonic. The limiting high temperature behaviour is 
also similar; Phillips’s modified spectrum gives ©, =392°K compared 
with ©,,=395 + 3°K obtained by extrapolation of the heat capacity results 
from the region T < @,/3. 

Marcus and Kennedy (1959)+ have used the neutron results of Brock - 
house and Iyengar (1958) to deduce other features of the vibrational 
frequency spectrum of germanium. In particular, they have been 
concerned with the dispersion of the lattice waves of low frequency and 
have estimated the coefficient B of eqn. (3.3). Their result agrees with 
that given in table 3 within the uncertainty assigned to the latter quantity. 

Rather less can be said about results for silicon because the frequency- 
wave number relations for it have yet to be measured. As is apparent 
from the discussion of the germanium results, these relations contain the 
most important information about vibrational properties. However, 
certain deductions about the frequency spectrum for silicon can be made. 
Brockhouse (1958) has suggested that to a useful degree of approximation 
the spectra of germanium and silicon are related by a constant factor. 
This may be seen from a reduced plot of ©, against 7 (fig. 9), but it is also 


7 We appreciate Dr, Marcus's courtesy in sending us a preprint, 
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evident that there must be small differences between the two spectra. 
Phillips (1959) has shown that apart from the change of scale only small 
changes in the positions of peaks need to be made to obtain a spectrum 
for silicon which is in good agreement with the measured heat capacity. 


Fig. 9 
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We can summarize this section by saying that the present heat capacity 
results for silicon and germanium are consistent with other measured 
vibrational properties. However, full use of the accuracy of the heat 
capacity measurements has not been made in the comparison with other 
results because of uncertainties inherent in dealing with anharmonic 
effects. All of the results are inconsistent with the Smith model (Smith 
1948) as it has been applied to silicon and germanium by Hsieh (1954). 
This inconsistency has been dealt with in detail by Brockhouse and 
lyengar (1958). On the other hand, the behaviour of the properties of 
diamond is quite different, as may be seen from fig. 9. It is not clear at 
present to what extent the Smith model is applicable to diamond because 
Smith’s detailed calculations (1948) were made before accurate values of 
the elastic constants of diamond were available. 


4.2. Anharmonic Effects 


In the foregoing sections we have suggested that the decrease of ©, 
of both silicon and germanium at higher temperatures is an anharmonic 
effect of a type which cannot be explained by change of volume alone, 
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It is unlikely that the effect is due to an under-estimation of C,-C,, 

As we have seen, C,, —C,, for both solids is small and is known sufficiently 
accurately. It seems certain also that the decrease of @, cannot be 
accounted for by the thermal creation of imperfections in the lattices ; 
the energies involved are much too large to give an observable effect on 
C,, in the temperature region below 300°c. Therefore, the only obvious 
remaining possibility is that, in the region 7’> ©,)/3, the thermal energy is 
sufficient to make the vibrations markedly anharmonic. When we 
examine this possibility in more detail, two difficulties become apparent. 

In the first place, one would expect the ratio of the amplitude of the 
vibrations to the interatomic spacing to be the relevant parameter ; 
this ratio is some 30% larger for silicon than it is for germanium, yet the 
departures of the @,,’s from the Thirring expansion occur at about the same 
reduced temperature for both. By contrast it should be mentioned that 
this parameter correlated reasonably well with similar departures found for 
alkali halide crystals (Barron et al. 1957). 

The second difficulty arises from the fact that the crystals have appre- 
ciable zero point energy (table 5). In the region of 7’'~ @,/3 the thermal 
energy is only about a third of the zero point energy. Hence, we must 
suppose that zero point energy also introduces anharmonicity in the 
lattice vibrations even at the lowest temperatures; this point has been 
stressed previously (Barron et al. 1957, 1959). In particular, it has been 
pointed out that while the analysis of experimental heat capacities gives 
properties of an effective harmonic frequency spectrum, this spectrum 
cannot be identified with the spectrum usually derived in theoretical work 
using only second order terms in the potential energy. 

The striking feature of the observed effect at high temperatures is the 
abruptness with which it appears (fig. 5). This is illustrated in another 
way in figs. 10 and 11 where differences between C,, computed from C, 
and C,, corresponding to the supposed harmonic behaviour of the crystals 
are plotted against the absolute temperature. Within the precision of 
the measurements AC, is a linear function of the temperature and dAC,/dT’ 
is nearly the same for both silicon and germanium. 

An extensive theoretical study of the effect of anharmonicity on the 
properties of crystalline lattices has recently been carried out by Ludwig 
(1958). His results predict for high temperatures an additional term in 
the heat capacity which is directly proportional to the absolute tempera- 
ture. A similar result is obtained (Peierls 1955) from a purely classical 
consideration of the problem and Eastabrook (1957), using a simplified 
treatment, has estimated that the magnitude of the anharmonic contribu- 
tion should be less than that of C,—C,. Reference to figs. 10 and 11 
and to table 2 shows that AC, is very much larger than C,—C,. More 
important is the fact that AC’, is not directly proportional to the absolute 
temperature. 

It would indeed be surprising if AC, did agree with the above predictions 
which should hold only in the general temperature region 7'> Op; the 
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highest temperature reached in the experiments was less than ©, even for 
germanium. Quantum effects in the region 7' < @,, will reduce the magni- 
tude of AC, below that given by classical mechanics but the precise tempera - 
ture dependence to be expected is not known. From the experimental 
evidence alone it is difficult to reach a definite conclusion regarding the 
exact temperature dependence. This would require either greater precision 
in the measurements or an extension of the measurements to higher 
temperatures. Ludwig’s work appears to provide the general basis for 
the treatment of anharmonicity but has so far only been applied in detail 
to a model of the face centred cubic lattice. A detailed comparison with 
our results awaits an analysis of Ludwig’s theory for the particular case of 
the diamond structure elements. 
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APPENDIX 


Measured Heat Capacities (units —cal/g atom deg) 


1. Silicon 
1 Ge T CG. T C, 
7-720 0-0008235 48-63 0:4957 183-00 3:470 
8-122 0-0009686 49-80 05241 184-61 3:492 
8-465 0-001090 50-50 05380 186-40 3:526 
9-001 0-001330 51-96 0-5741 188-08 3550 
9-471 0-001541 54:38 0-6315 189-94. 3:582 
9-941 0-001814 55:43 00-6575 191-49 3:607 
10-455 0-002133 58-09 0-:7271 193-25 3°631 
10-894 0-002439 59-01 0-7488 194-86 3°655 
11-501 0-002889 61-66 0-8136 196-69 3-684 
11-908 0-003274 62:71 0:8396 198-18 3:707 
12-435 0-003671 65:10 0-9001 200-08 3:738 
12-592 0-003926 66:27 0-9309 201-45 3°760 
12-649 0-004151 68:67 0:9852 203-42 3°787 
12-927 0-004176 69-74 1-008 204-95 3°804 
12-995 0-004283 72:16 1-069 206-71 3°827 
13-032 0-004467 73°32 1-101 208-67 3:862 
13-618 0-005165 75:57 1-152 210-03 3:877 
14-590 0-006620 76:82 1-180 212-34. 3:921 
14-94] 0:007559 80-40 1-269 213-38 3-924 
15-527 0-008367 83:95 1-359 216-88 3:970 
15-302 0:008035 86:29 1-413 220-56 4-017 
15-420 0-008216 89°73 1-502 224-19 4-063 
16-256 0-01004 93-20 1-579 227-78 4-107 
16-505 0-01050 96-71 1-662 231-26 4-149 
17-221 0-01262 100-07 1-740 231-33 4-147 
17-529 0-01330 101-46 lene, 234-77 4-186 
17-694 0-01380 103-46 1-819 238-24. 4-230 
18-489 0-01641 106-87 1-902 241-68 4-264. 
18-686 0:01737 110-15 1:977 245-09 4-309 
19-543 0-02038 113-53 2-054. 248-64. 4-338 
19-551 0-:02058 117-02 2-137 252-36 4-379 
19-606 0-02094 120-39 2-214. 256-03 4-418 
20-253 0-02416 123-85 2-293 259-68 4-460 
21-967 0-:03348 124-65 2-311 263-29 4-489 
22-101 0:03465 128-05 2-388 266-88 4-521 
22-959 0-04054 131-42 2-463 270-43 4-558 
23-058 0-041 11 134-75 2-537 DYBOV 4-582 
25-17 0-05811 138-00 2-603 277-38 4-621 
26-47 0:07150 141-16 2-677 280-85 4-646 
28-48 0:09537 144-42 2-744. 283-36 4-652 
29-80 0-1121 147-78 2-815 284-30 4-680 
31-62 0-1387 151-05 2-879 287-71 4-696 
32-90 0-1594. 154-26 2-944 287-68 4-712 
34-67 0-1903 157-51 3-010 291-08 4-735 
36-15 Deans 159-72 3°052 292-90 4-733 
BoD 0-2490 160-81 3-076 294-69 4-767 
39-43 0-2841 162-97 3-116 295-52 4-770 
40-88 0-3140 166-35 3:178 296-63 4-770 
42-59 0-3519 169-87 3:243 296-73 4-764 
43-96 0:3827 173-16 3°300 298-41 4-785 
45-75 0:4249 176-39 3°358 300-33 4-790 


46-92 0-4522 Toad 2 3-423 300-45 4-796 


Il. Germanium 


C T 6; - 0, 


p p 


0-0001283 24-548 0-3950 123-84 3°928 
0-0001766 25-86 0-4504 127-36 4-006 
0:0002337 27°22 0-5094 130-81 4-079 
0-0003008 28-46 0-5631 134-18 4-150 
0-0003918 29:87 0-6267 135-93 4-178 
0-0004876 31-16 0-6839 139-22 4-245 
0-0005763 32-68 0-7509 142-42 4-308 
0-0007935 33°91 0-8049 145-75 4-363 
0-0008362 35-72 0-8848 149-23 4-427 
0-001077 36-93 0-9376 152-66 4-475 
0-001170 39-00 1-026 156-05 4-529 
0-001446 40-46 1-087 159-39 4-577 
0-001636 42-43 1-168 162-69 4-630 
0-001915 44-29 1-240 166-06 4-668 
0-002188 45-89 1-312 169-50 4-712 
0-002513 47-87 1-393 172-90 4-753 
0-002802 49-37 1-453 176-26 4-796 
0-003247 51-90 1-557 178-49 4-824 
0-003845 53-59 1-622 183-82 4-880 
0-004668 55-63 1-703 189-31 4-931 
0-005757 57-15 1-765 192-91 4-967 
0-007220 59-02 1-840 196-48 4-999 
0-008950 60-39 1-895 200-02 5-034 
0-01104 62-35 1-974 203-54 5-065 
0-01391 63-76 2-029 207-02 5-091 
0-01679 65-81 Zit 210-49 5-119 
0-02103 67-30 2-169 213-92 5-145 
0-02479 69-22 2-239 217-34 5-173 
0-03039 70-73 2-293 222-82 5-219 
0-03512 72-54 2-363 226-53 5-232 
0-03559 74-07 2-421 230-21 5-256 
0-04187 76-00 2-491 233-87 5-286 
0-04265 7:45 2-546 237-50 5-299 
0-04946 79-38 2-617 241-11 5°325 
0-05829 81-09 2-681 244-70 5:349 
0-06767 82-75 2-743 246-40 5-349 
0-07647 84:36 2-801 249-98 5-374 
0-08112 85-97 2-859 253°51 5°390 
0-08827 87-50 PN! 257-02 5-415 
0-09754 89-06 2-973 262-81 5437 
0-09828 90-70 3-015 268-85 5-461 
0-1127 92-22 3°057 272-84 5-478 
0:1231 93-98 3-120 276-81 5-494 
0-1402 95-46 3165 280-72 5-510 
0-1399 STL 3-216 280-90 o-old 
0-1534 98-76 3-266 284-80 0-025 
0-1675 100-27 3309 286-55 5:537 
0-1701 103-23 3-392 288-66 5-546 
0-1826 106-87 3-502 292-38 D°957 
0:2027 110-39 3-596 292-48 5-560 
0-2256 113-82 3-682 296-23 5-581 
0-2709 i bol br’ 3-768 296-24 5-972 
0-3392 120-44 3-849 300-01 5-609 

300-05 5-600 


Nort: The number of significant figures given reflects the precision of the 
measurements rather than their absolute accuracy. 
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ABSTRACT 


The equilibrium positions of dislocations in a number of types of array 
have been found numerically. Arrays considered are: 

(i) that generated by a dislocation source which is between two locked 
dislocations, 

(ii) a dislocation pile-up against a locked dislocation which is not on the 
slip-plane, 

(iii) a group of positive dislocations attracted to a group of negative disloca- 
tions but prevented from coalescing by locked dislocations, 

(iv) a group of positive dislocations attracted by a group of negative 
dislocations on a parallel slip plane. 

A short table of the zeros of the first derivative of the Laguerre 
polynomials is given. 


§ 1. INTRODUCTION 


A CENTRAL problem in the dislocation theory of solids is the equilibrium 
arrangement of a number of mobile dislocations under their mutual 
interactions and the action of a stress field. As this is a many-body 
problem exact analytical solutions cannot be expected although Eshelby 
et al. (1951) have shown that there are such solutions for a restricted class 
of problem. Their solutions follow from properties of the classical 
orthogonal polynomials which are such that to a polynomial of the nth 
degree there corresponds an arrangement of 7 like dislocations, all on the 
same slip-plane and located at the zeros of the polynomial. These 
dislocations are in equilibrium under an applied stress field related to the 
coefficients of the differential equation satisfied by the polynomial. Some 
of the sets of orthogonal polynomials correspond to dislocation arrange- 
ments of great interest but it does not appear possible to extend this 
method as it is essentially an inverse approach to the problem. It has 
been suggested (Leibfried 1951, 1954, Haasen and Leibfried 1954) that 
approximate solutions can be obtained by replacing the discrete dislocations 
by a continuous dislocation density. In this case the dislocations need not 
be all of the same sign and the equations of equilibrium become an integral 
equation which, if the mobile dislocations are all on the same slip plane, 
can, be solved by routine procedures (Head and Louat 1955). 

There are still many problems for which neither of the above methods 
is convenient and for some of these the equations of equilibrium have 
been solved numerically on the electronic digital computer CSIRAC. 


+ Communicated by the Author. 
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§ 2. GENERAL METHOD 


All dislocations are taken to be of infinite length parallel to the z axis 
with slip plane parallel to the x axis and in most cases the slip plane is 
y=0. For each array a characteristic length is taken as the unit of 
distance and the unit of stress defined as the stress at unit distance along 
the slip plane from an isolated dislocation of the type (screw or edge) 
being considered and with the sign convention that a positive stress moves 
a positive dislocation to the left. In these units the equations of 
equilibrium for positive dislocations at the points x1, %, ...x,, all on 
the same slip plane are 


$1 

j=1 Uji %y 

#j 
where P() is the stress field due to any applied stress, locked dislocations, 
etc. In those cases where there are also mobile negative dislocations the 
positions of these are always related by symmetry to the positions of the 
positive dislocations and their presence can be allowed for by extra terms 
in F(z). 

This set of equations was solved by relaxation. Starting with 

approximate positions for the dislocations, each dislocation in turn was 
moved a distance given by Newton’s rule as 


1 Lo \ Geta 
ree {2 U5 — Uj +Peyh [tS ee) - Pa) 


After all dislocations had been moved, any subsidiary conditions were 
adjusted and the cycle repeated until a sufficiently accurate solution was 
obtained. In the results given here, the least significant figure has some, 
but not necessarily full, significance. Incidentally, this is a possible 
method of computing the zeros of the classical orthogonal polynomials 
but as it gives only first order convergence it is slow compared with any 
of the usual methods for finding the zeros of a polynomial. 


+ P(x,;)=0, fee es oat 


§ 3. DisLocaTion SOURCE BETWEEN LocKED DIsLocATIONS 


The dislocation array of fig. 1 (a) is generated by a dislocation source S 
which, under the influence of the applied stress o, emits pairs of positive 
and negative dislocations. These move off in opposite directions along 
the slip-plane until the leading pair of dislocations become locked. As 
all dislocations are in the same slip-plane the same equations hold for 
either edge or screw. The dislocation source is taken to be midway 
between the locked dislocations and the unit of length such that the 
locked positive dislocation is at «=0, the dislocation source at x=1, and 
the locked negative dislocation at x=2. Since there are equal numbers x 
of positive and negative dislocations then for a positive dislocation at 
«=x, by symmetry there is a negative dislocation at 2=2—a;,. 

The value of o is taken to be equal and opposite to the back stress of 
all dislocations on 8 so that the nett stress at Sis zero. For, if o is greater 
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than this value the source will tend to generate more dislocation pairs and 
if it is less then dislocation pairs will tend to run back into the source and 
be annihilated. The equilibrium positions of the positive dislocations and 
the corresponding value of o are given in table 1 for n=2 to 6. 


Fig. | 


fa\ne Oa ae. enc, 


fom 
(hie c= 


(dd) 2 4 ae 


(e) + +Q@Or 2 
ai ae 
(f) T we ae 
=e i alin 


(€) + Tr ae 


Dislocation arrays for »=3. Locked dislocations are circled. 
(a): Table 1 and § 3. 
(6): Table 2 and § 4. 
(c): Table 3 (a) and § 4. 
(d): Table 3 () and § 4. 
(e): Table 5 and § 5. 
(f). (9): § 6. 
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If this problem is solved by Leibfried’s approximation (which should 
be reasonable for n sufficiently large) it is found that o/n=za. For the 
exact solution this quantity increases from 2-40 for n=2 to 2-79 for n=6. 


Table 1. Positions of Positive Dislocations and Critical Stress for Fig. 1 (a). 
Negative Dislocations at 2—# 


n 2) 3 4 5 6 


o |4-7956 |7-7199 | 10-701 13-714 16-749 


we I) 0 0 0 ) 
0-28459 | 0-10958 | 0-:057832 | 0:035657 | 0-024153 
0:-42426 | 0:20751 0-12448 0:083194 
050672 0:28246 0-18351 
0:56223 0-34076 
060273 


§ 4. Pinn-up Against Lockep DisLocaTION NOT ON SLIP-PLANE 


One of the problems solved by Eshelby et al. was the piling up of n like 
dislocations against a locked dislocation on the slip-plane by a constant 
applied stress. We consider here the related problem where the locked 
dislocation is off the slip-plane. In this case for a pile-up of m mobile 
dislocations, the applied stress o must be less than some critical value. 
otherwise some dislocations will be forced past the locked dislocation and 
be lost from the pile-up. It is necessary to consider screw and edge 
dislocations separately. 

If the slip-plane of the mobile screw dislocations is y=0 and the 
displacement of the locked dislocation from this plane is taken as the unit 
of length, then the position of the locked dislocation can be taken as 
(0,1), (fig. 1(6)). The repulsion from the locked dislocation reaches a 
maximum on the plane y=0 at «=1 (Cottrell 1953) and the critical case 
is when the leading mobile dislocation is at this point. The applied 
stress necessary to produce this critical case and the positions of the 
dislocations are given in table 2 for n= 1 to 6. 

The same problem for edge dislocations has two critical cases due to 
the more complicated variation of the stress field of the locked edge 
dislocation at (0, 1) on the slip-plane y= 0 (Cottrell 1953). As the applied 
stress is increased the n mobile dislocations move towards the locked 
dislocation until the leading dislocation is at ~=1+ y/2, (fig. 1(c)). The 
repulsion from the locked dislocation is a maximum at this point and any 
further increase in stress causes the leading dislocation to move past the 
locked dislocation where it is then held by an attractive force. The stress 
can then be increased until the leading dislocation reaches x=1— 4/2 
(fig. 1(d)) where the attraction is a maximum, any further increase in 
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stress causing the leading dislocation to move off to infinity. The stress 
and the positions of the dislocations for these two critical cases are given 
in tables 3 (a) and 3 (0) for n=1 to 6. 


Table 2. Positions of Dislocations and Critical Stress for Fig. 1 (0) 


0-5 | 0-32340 | 0-24016 | 0-19128] 0-15901 | 0-13610 


1 1 1 1 
6-6625 6-2328 6-0760 5:9989 5:9549 
15-547 14:377 13-895 13-639 
28-045 25-822 24-831 
44-259 40-693 
64-252 


Table 3(a). Positions of Dislocations and Critical Stress for Fig. 1 (c) 


n il 2 3 4 Ip 6 


a | 0-25 0-16103 | 0-11943 | 0-095064 | 0-079005 0-067609 


x | 2-4142| 2-4142 2-4142 2-4142 2-4142 2-4142 


13-654 12-808 12-499 12-348 12-262 
31-509 29-184 28-222 27°713 
56-677 52-220 50-234 
89-324 82-162 
129-58 


Table 3(b). Positions of Dislocations and Critical Stress for Fig. 1 (d) 


n 1 2 3 4 5 6 

a 0-25 0-16702 0-12536 0-10032 0-083617 0-071681 

ES —(-4142 | —0-4142 —(0-4142 —(0-4142 —(-4142 —(0-4142 

+ 11-638 + 10-682 + 10-328 + 10-152 + 10-052 

+ 29-384 + 26-184 + 25-198 + 24-672 

+ 52-251 +47-886 +45-928 

+ 82-950 + 76-050 

+ 120-78 
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When the locked dislocation is on the slip-plane Eshelby e¢ al. have 
shown that for a suitable choice of the unit of length, the mobile 
dislocations are at the zeros of the first derivative of the (n + 1)th Laguerre 
polynomial (table 4). Starting with such a pile-up and keeping the 
applied stress constant suppose the locked dislocation is moved off the 
slip-plane a distance d such that dislocations are about to escape from the 
pile-up. Then dis given by d/D=<c,,/l,, where o,, is the critical stress given 
in table 2 or 3(b), J, is the smallest zero of L,,,,'() (table 4) and D was 
the distance of the leading mobile dislocation from the locked dislocation. 
when the latter was on the slip-plane. As m increases from 1 to 6 this 
ratio increases from 0-5 to 0-516 for screw dislocations and from 0-25 to 


Table 4. Zeros of L,,4,'(%) 


ip tl 2 3) 4 Is 6 


x | 2 | 1-2679 | 0-93582 | 0-74329 |] 0-61703 | 0-52766 
4-7321 | 3-3054 2-5716 2-1129 1-:7963 
77589 5:7311 4-6111 3°8767 

10-954 8:3983 6-9194 

14-261 11-232 

17-648 


0:272 for edge dislocations. ‘These calculations are for ideal elastic 
dislocations and the analysis of Leibfried and Liicke (1949) suggests that 
the finite core of a real dislocation could be allowed for by decreasing the 
above values of d by the core diameter. Stroh (1954) has considered 
dislocation pile-ups in which the leading dislocations are so close together 
that their cores are overlapping. The present work shows that for such 
a pile-up to be possible the locked dislocation must be exactly on the 
same slip-plane as the following dislocations and the pile-up could not be 
formed by a dislocation source which climbs one atomic plane for each 
dislocation loop generated (Thompson 1953). 


§ 5. PosrTivE AnD NEGATIVE DisLocaTIons on Same Siip-PLANE 


We now consider a group of ” positive dislocations attracted to a group 
of n negative dislocations on the same slip-plane but prevented from 
coalescing by the locking of the leading positive and negative dislocations 
(fig. 1(e)). If the origin and the unit of length are chosen so that the 
locked dislocations are at x= +1, then the positions of the positive 
dislocations are given in table 5 for n=2 to 6, and to each positive 
dislocation at «=a, there is a negative dislocation at x= —z,. By 
considering the change in the potential energy of this dislocation array 
for a small movement of a locked dislocation, the argument of Cottrell 
(1949) can be adapted to show that the stress at a locked dislocation, due 
to all other dislocations, is exactly 4n. 
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Figure 1(e) is also the image representation of the case of n positive 
dislocations attracted to a free surface at «=0 but being held up by a 
locked dislocation. For screw dislocations such an image representation 
is exact (Head 1953a) but is only an approximate representation of the 
interaction of edge dislocations with a free surface (Head 1953 b). 


Table 5. Positions of Positive Dislocations for Fig. 1 (e). 
Negative Dislocations at — x 


n 2 3 4 5 6 
a ba? 1 1 ] 1 
4-2361 | 2:6803 |] 2-1427 18676 1-7000 
15-088 7:7678 54101 4.°2593 
44-633 20-300 12-920 


116-88 48-695 
280-44 


§ 6. PostrivE AND NEGATIVE DISLOCATIONS ON PARALLEL SLIP-PLANES 


Louat (1953) and Fujita (1954, 1958) have considered the problem of 1 
mobile positive edge dislocations on one slip-plane attracted by n mobile 
negative edge dislocations on a parallel slip-plane. Fujita argued that 
when the leading dislocations came sufficiently close there would be a 
repulsion between them and so there would be an equilibrium arrangement 
as in fig. 1(f) and that an external applied stress (which he estimated 
approximately) was necessary to cause the dislocations to interlace as in 
fig. 1(9). 

It was intended to calculate this stress more accurately but from the 
results of §5, in particular the magnitude of the stress at the locked 
dislocation, it appeared probable that in Fujita’s problem the attraction 
of the outer dislocations would outweigh the repulsion of the leading 
dislocations and that the dislocations would interlace at least partially in 
the absence of an applied stress. This was confirmed by calculation for 
the case n=2 and so, a fortiori, for n>2. The distance between the slip- 
planes was taken as the unit of length and the origin of coordinates at 
the centre of inversion of the dislocation array so that to a positive 
dislocation at x=2;, y= +4 there was a negative dislocation at v= —2;, 
y= —4. The leading pair of dislocations were locked at various values of x 
and the equilibrium position of the other dislocation pair was found. From 
this the stress field at the locked positive dislocation was calculated. The 
first five entries of table 6 are for the locked dislocations in the neighbour- 
hood of Fujita’s equilibrium position, and it will be seen that although the 
stress on the locked dislocation has a minimum it does not become zero 
and so there is no position of equilibrium. The locked dislocations were 
then moved step by step in the direction of the stress acting on them and 
it was found that after the dislocations had interlaced each positive 
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dislocation paired off with a negative dislocation and that between the 
two dislocation pairs there was a weak repulsion, the last entry of table 6 
being typical. 

Table 6. Fujita’s problem, n= 2 


a ee era 


Locked Positive Mobile Positive Stress at Locked Positive 
Dislocation Dislocation. Dislocation 
c= t= 
0:375 2-9877 0:-47716 
0-3125 28355 0-43506 
0-25 26958 0:41035 
0-1875 2-5695 0-41813 
0-125 92-4572 0-47009 
— 4-0 5:0289 0-00569 


For n=2 then there is no equilibrium arrangement of the dislocations. 
For n>2 it would appear that the dislocations will interlace, at least 
partially, without an applied stress and that if they interlace completely 
they split up into positive-negative pairs which are weakly repellent. It is 
possible that a stable array is formed with the dislocations partially 
interlaced where a number of dislocations pairs would be formed but 
prevented from drifting apart by the attraction between the outer 
unpaired dislocations but the present computer programme is unsuitable 
for examining this question. 
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ABSTRACT 


Stress-strain tests have been made on Li, Na and K at temperatures 
between 293 and 4:2°x. At higher temperatures the behaviour of all three 
metals is very similar. They are soft and after the yield point very little 
hardening is observed. At low temperatures Li and Na transform to f.c.c. 
and h.c.p. structures respectively on deformation, the character of their 
stress-strain curves changes accordingly and they harden continuously 
throughout the test. At 20 and 4:2°x Li deforms with many regions of 
discontinuous slip and at 4-2°x it fractures in a brittle fashion. Whilst K 
remains b.c.c. down to 4:2°xK it also exhibits an unusual type of hardening 
at 20 and 4:2°K which is very temperature dependent. The change in the 
electrical resistivity during low temperature deformation was also measured. 
There was no detectable change when K was deformed at 4:2°xK. This offers 
further evidence that a large proportion of the resistivity change which is 
observed when metals are deformed is due to stacking faults. Both Li and 
Na give considerable resistivity changes when they are deformed at 4:2°K 
and much of this change is ascribed to the scattering of electrons at the 
faulted regions within the transformed phases. 


§ 1. INTRODUCTION 


In the past, experiments on the alkali metals have been very limited, 
mainly because of their high reactivity and the consequent difficulty in 
the handling and preparation of specimens. They do, however, provide 
an extremely interesting field for study. At room temperature they all 
have a body-centred cubic (b.c.c.) structure and are very different from 
the b.c.c. transition metals which have high melting points and which 
are susceptible to low temperature brittleness. The properties of sodium 
and lithium are modified at low temperatures, however, because they 
undergo martensitic transformations and this puts their behaviour in 
marked contrast to that of potassium which remains b.c.c. down to the 
lowest temperatures. 

In this paper we shall describe some experiments on the way in which 
the tensile properties of lithium, sodium and potassium change at 
temperatures down to 4:2°K. In addition to this the change in the electrical 
resistivity during deformation has been measured. 


+ Communicated by the Authors. 


304 D. Hull and H. M. Rosenberg on the Deformation 


§ 2. THE EXPERIMENTS 

The tensile tests were made using a hard beam machine which is shown 
diagrammatically in fig. 1. It is a modification for low temperature use 
of one described by Adams (1958). A constant speed driving motor M 
lowers the cross beam C via a worm drive, W, and the central threaded 
rod, R. This extends the specimen, S, because its lower end is connected 
to the rods P which are in turn fixed to the cross beam. ‘The upper end. 
of the specimen is connected to a load cell, L, by a rod which passes through 
the centre of R. The sensitivity of the load cell was 5-6kg per mv. This 
output potential was fed to an automatic recorder and the specimen was 
extended at a constant rate of 0-325 mm per minute. A dewar of coolant 
could be placed around the whole of the lower section of the apparatus. 
The liquid air container, A, which was attached to the rods P acted as 
a radiation shield when liquid helium was used in the dewar. ‘The moving 
joints were sealed with flexible rubber tubing to prevent loss of helium 
gas. All the low temperature section of the apparatus was made of stainless 
steel. 

The specimens were made by extruding the metal through a circular 
nozzle under liquid paraffin. The metals rapidly recrystallize at room 
temperature so that large grained (>1mm) polycrystalline specimens 
were formed. The specimens were about 5cm long and they had a 
diameter of 0-165cm for the electrical resistivity experiments. The 
stresses required to extend these specimens were too small to measure 
accurately, however, and so for the tensile tests a diameter of 0-47 cm 
was used. Each end of the specimen was supported in a holder made of 
a small block of metal with a conical hole in it (fig. 1, inset), the smallest 
diameter of the hole being slightly larger than the specimen diameter. 
About 2 cm of the specimen was passed through this hole and was packed 
down into the cone, so that when the specimen was pulled the end would 
not slip out. For the electrical resistivity experiments current and 
potential leads of fine copper wire were forced into the specimen, the 
metal being pressed around the wire to ensure good contact. The potential 
leads were inserted into the specimen through small holes drilled in the 
sides of the cones whilst the current leads were pushed into the ends of the 
specimen. Such contacts seemed to be very satisfactory, probably due 
to the relatively large contraction of the specimens on cooling. Currents 
of the order of 1 amp were used and the potential produced across the 
specimen was measured on a galvanometer amplifier. 

The section of the apparatus containing the specimen could be surrounded 
by a dewar vessel and for experiments below room temperature baths of 
dry ice in acetone (195°x), liquid air (~80°x), liquid hydrogen (20-4°xk) 
and liquid helium (4:2°K) were used. 

The specimens were fitted into their holders and into the apparatus 
whilst immersed in a bath of liquid paraffin. For experiments below room 
temperature the bath was removed and the film of paraffin washed off 
with carbon tetrachloride just before the specimen was immersed in the 
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The hard beam tensile apparatus used in the experiments. The inset shows the 
specimen holders used for the alkali metals. 
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coolant. Whilst this technique prevented large scale oxidation, thin 
oxide films could not be avoided. 

The lithium used was 99-7% pure. Both commercial purity sodium 
(B.D.H.) (impurity not greater than 0-1%) and high purity vacuum 
distilled metal (99-9989) were used. All the potassium specimens were 
made from vacuum distilled metal (99-999%). 


§ 3. THe TENSILE TESTS 
The stress/strain curves at the various temperatures are shown in 
figs. 2 to 4 and values of the yield stress are shown in fig. 5. In all cases 


the yield stresses were very small and difficult to measure accurately 
with the present machine. 


Fig. 3 
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The stress-strain curves for sodium taken at 195, 80, 20 and 4-2°K. 
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Fig. 4 
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‘The stress-strain curves for potassium taken at 80, 20 and 4:2°K. Note the 
very sharp increase in the rate of linear hardening between 20 and 4:2°K. 


Fig. 5 
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The variation of the yield stress with temperature for lithium, sodium and 
potassium. 
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The stress/strain curves of the metals at lower temperatures are quite 
different from one another. Both sodium and lithium work harden very 
rapidly to relatively high stresses at 20 and 4:2°K with elongations of 
about 40%. At both these temperatures, and at 4-2°K in particular, 
lithium shows a series of very marked yield drops. At 20° this dis- 
continuous yielding was confined to the second half of the test but at 
4:2°K, yield drops were detected in some cases before the yield point was 
reached. At 4-2°xK the specimens did not neck down on fracture, but they 
snapped in a brittle fashion after an extension of about 30-50% showing 
a bright fracture face. No signs of any discontinuous yielding were noted 
in either potassium or sodium. They were extremely ductile at all testing 
temperatures and the specimens necked down at fracture to a knife edge. 
Nevertheless, whereas in sodium the hardening continued throughout 
most of the test, in potassium the initial work hardening rate was not 
maintained and the maximum stress at 4:2°K was only 20 to 25% that 
of sodium. One result of this is the smaller elongation of potassium. 
There is a marked change in the curves of potassium from 20 to 4:2°x. 
The initial work hardening rate is much higher at 4:2°K than at 20°K 
where there is a long region of almost linear hardening. 


§ 4. THe Evecrricat RESISTIVITY 


The change in the electrical resistivity Ap during deformation was 
measured at 4-2°K for all three metals and in addition at 20°x for lithium 
and sodium. A large part of the resistance change is of course due to the 
increase in the form factor of the specimen as it is being extended. After 
this has been allowed for, however, there was still a considerable increase 
in the resistivity of both sodium and lithium after extension in liquid 
helium. There was, however, practically no change in sodium after 
deformation at 20°K although if the dependence of Ap on strain had been 
of the same magnitude as it was at 4-2°K this could have been observed 
quite easily. In actual fact a detailed examination of four different runs. 
shows that for the first few per cent of strain the resistivity actually decreases. 
by 2 or 3%. It then begins to increase again, reaching its original value 
by about 10% strain and thereafter it continues to increase very slightly. 
No detectable change (less than 1%) in the resistivity of potassium could 
be detected when it was deformed at 4:2°K. Since it was so soft it could 
only be extended by a very small amount (3-4%) before it started to neck 
down and break, and hence very few measurements of the resistivity 
could be taken and these only after very small amounts of damage had 
been introduced into the specimen. Nevertheless, it was quite clear 
that even if there was a slight resistivity change it was very much less 
than for the other two metals. Logarithmic plots of Ap against the 
strain, <, are shown in fig. 6. It will be noted that the relation is approxi- 
mately linear and is of the form Ap=Aec”ohmcm where 

A=2:4x10-7, n=0-92 for lithium at 4:2°K, 
A=1-65x10-?, n=0-92 for lithium at 20°, 
and (Aa 2x10, n=1-1 for sodium at 4:2°K. 
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Whilst there was hardly any change in the resistivity of sodium after 
deformation at 20°K it was noted that if the specimen was warmed up to 
80°K after this extension and then re-cooled to 20°K that its resistance 
had increased by about 10%. This experiment was repeated several times 
with approximately the same result in each case, independent of whether 
the specimen was warmed and re-cooled whilst still under tension or with 
the load removed. On warming to room temperature and re-cooling, 
nearly all this resistivity increase disappeared. The resistivity change 


Fig. 6 
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A logarithmic plot of the increase in resistivity, Ap due to deformation, e¢, 
against the strain for lithium at 20 and 4-2°K, and for sodium at 4:2°K. 
No change in resistivity could be detected when potassium was deformed 


S é : ‘ 
at 4-2°K and only a very slight change (see discussion) when sodium was 
deformed at 20°K. 


in lithium after deformation at 20°K was similar in form to that obtained 
when the metal was deformed at 4-2°K except that the value of Ap for 
a given strain was about 30% less at 20°K. 
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The ratios of the resistivities between 4-2°K and room temperature 
P4-2/Pprp Were 9-3 x 10-3 for lithium, 2-7 x 10-4 for high purity sodium, and 
2-4 x 10-8 for potassium. 

§ 5. Discussion 
5.1. The Tensile Tests 

The general form of the stress/strain curves can be explained in terms 
of the transformations which lithium and sodium undergo at low 
temperatures. Barrett (1948, 1955, 1956) has shown by x-ray diffraction 
experiments that the following transformations occur: 


lithium goes from b.c.c. to h.c.p. below 72°K without cold work, 
from b.c.c. ; ° 
to f.c.c. (faulted) if cold worked below 110°x, 
or h.e.p. 
sodium goes from b.c.c. to h.c.p. (faulted) below 36°K without cold 
work, if cold worked, below 51°K. 


The x-ray data suggest that if the sodium is unworked, about 7% is 
transformed by 5°K. Severe cold work at 5°K transforms about 50%. 
Potassium remains b.c.c. down to 5°K and no transformation is observed 
even if it is severely cold worked at that temperature. 

Thus we can take the results for potassium (fig. 4) as being typical 
of the b.c.c. lattice of the alkali metals. The curve for 80°K shows the 
general behaviour at relatively high temperatures. Beyond the yield 
point there is a short region of slight hardening. After this, little further 
hardening occurs before. fracture. This is the same type of curve as is 
obtained in a b.c.c. transition metal beyond the region of the lower yield 
point. This last effect has not been detected in the alkali metals. Similar 
curves are shown for lithium at 293° and 195°K and for sodium at 80°K 
at which temperatures these metals are still b.c.c. At lower temperatures, 
however, these two metals transform and this is reflected in their 
stress/strain curves by the fact that throughout the test the metal 
continues to harden. The precise nature of this hardening is not yet clear. 
Two main possibilities suggest themselves. Firstly, a hardening due to 
the hold-up of dislocations at the boundaries between the phases and 
secondly, a hardening within the faulted phase itself. 

In sodium there is a considerable difference between the 20 and 4-2°K 
curves. This is because the transformation temperature is relatively 
low and hence much less transformation product will be formed at 20 
than at 4:2°K. 

In lithium the transformation takes place during deformation at 80°K, 
as is shown by the continuous hardening during the test at that temperature. 
By 20°x the transformation to the h.c.p. lattice is practically complete, as 
is shown by the similarity of the work hardening at 20 and 4-2°K. Whilst 
the general form of these curves is similar to those for sodium, there is 
of course the very considerable difference that discontinuous yielding is 
only observed for lithium. In sodium the transformation on cold work is 
from b.c.c. to h.c.p., whereas the transformation for lithium at 20 and 
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4-2°K is from h.c.p. to f.c.c. It is suggested, since no serrations are 
observed for sodium, that the discontinuous yielding in lithium is due to 
a catastrophic change from h.c.p. to f.c.c. 

It should be noted that a strain of a few per cent would be produced if 
the whole of the specimen transformed from h.c.p. to f.c.c., and whilst 
this would account satisfactorily for the integrated strain produced in 
the regions of discontinuous yielding, the large total strain must result 
from other forms of deformation. 

The curves for potassium also show some interesting features. Whereas 
at 80°K very little hardening is observed, by 20°K a very long region of 
practically linear hardening occurs. At 4:2°K this linear hardening has 
become so great that it is practically indistinguishable from the elastic 
part of the deformation. The mechanism for such an extemely temperature 
dependent hardening is not yet clear and it would not seem profitable to 
propose one without more detailed experiments. In the transition group 
b.c.c. metals such hardening (if it exists) is presumably masked by the 
yield point effect. It should be noted that beyond the region of rapid 
hardening at 4:2°K a softening mechanism occurs (presumably due to the 
breakdown or the avoidance of dislocation barriers analogous to stage III 
in f.c.c. metals) so that the stress/strain curve is practically horizontal. 
This softening is probably the same as that which occurs in the horizontal 
region at 80°K and also in the higher temperature curves of sodium and 
lithium. Thus it seems that in a pure b.c.c. metal at low temperatures 
there is a region of intense hardening after the metal has yielded. As 
the temperature is raised this hardening rapidly decreases and at still 
higher temperatures it becomes so small that once the metal has yielded 
the hardness remains approximately constant right up to fracture. It 
should be noted, however, that in spite of the very large increase in 
hardness in potassium between 20 and 4-2°K, a comparison of the ordinate 
scales of figs. 2 to 4 shows that even at 4-2°K the metal is still very soft 
compared with sodium and lithium. 

As the temperature is reduced the yield stress of all three metals increases 
in a similar manner tc that observed in the b.c.c. transition metals. The 
fact that even in the transformation region of lithium and sodium the 
yield stress continues to rise would seem to suggest that we are measuring 
the yield stress of the b.c.c. phase since previous workers have shown that 
the yield stress of h.c.p. and f.c.c. metals is not very temperature dependent. 


5.2. The Electrical Resistivity 

The increase in the electrical resistivity during deformation at 4:2°x 
is usually considered to be due to three main causes: the scattering of 
the conduction electrons by (1) the point defects produced during the 
deformation, (2) the dislocations, and (3) the stacking faults if the 
dislocations can dissociate into partials with a ribbon of fault between 
them. Considerable experimental evidence is now available (e.g. Christian 
and Spreadborough 1956, Seeger 1956, Lomer and Rosenberg 1959) to 
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bear out the suggestion first made by Broom (1952) that in metals in which 
the dislocations are able to separate into partial dislocations nearly all 
the extra resistivity associated with the dislocations is due to the scattering 
of the electrons by the stacking faults, rather than scattering by the 
dislocations themselves. This hypothesis is further supported by the 
present experiments in which no resistivity increase could be detected 
when potassium was deformed at 4-2°K, because stacking faults are not 
thought to occur in the b.c.c. lattice and hence little or no resistivity change 
should be expected. The fact that no resistivity change is observed would 
also suggest that in the b.c.c. lattice very few point defects are produced 
on deformation unless, for example, they are able to diffuse away at 4-2°K 
in potassium. This would require an activation energy of motion of very 
much less than 0-O0lev, which seems rather unlikely, unless a crowdion 
mechanism is invoked. 

One cannot be certain about the cause of the resistivity increases which 
are detected in lithium at 20 and 4-2°K and in sodium at 4:2°K. The x-ray 
experiments of Barrett (1948, 1956) show that in deformed specimens 
the regions of h.c.p. and f.c.c. are heavily faulted and there would seem 
to be no doubt that in addition to point defect resistivity some of the 
increases must be due to the scattering of the electrons by these faults. 
But there is also the possibility that the electrons are scattered at the 
interface boundaries between the phases, or by the strain field associated 
with those boundaries. In actual fact it is not thought that this would 
be a very large effect since the area of the interface boundary will be very 
small compared with the area of stacking fault. 

It should be noted that the resistivity change for sodium is about ten 
times less than that for lithium at any given strain, but that the strain 
dependence is very similar. This might be due to the fact that very 
much less of the material is transformed, particularly if we assume that 
there is no contribution to the deformation resistivity from the b.c.c. 
parts of the specimen. (Barrett only got a 50% transformation and his 
method of cold work was much more drastic than ours.) Another 
possibility is that the stacking fault probability might be much greater 
in f.c.c. lithium than in h.c.p. sodium. 

The resistivity measurements taken when sodium was deformed at 20°K 
are rather interesting. It will be recalled that there seemed to be very 
little change in the resistivity when the specimen was extended but that 
after it had been warmed up to 80°xK and then re-cooled to 20°K the 
resistivity increased. A tentative explanation of this behaviour can be 
given in the light of the recent work by Gugan and Dugdale (1958) on 
the change in the resistivity of sodium when it transforms. The resistivity 
of a metal can be divided into two parts. The first is the residual 
resistivity, pp, produced by the scattering of the electrons by the static 
defects in the specimen and which accounts for practically the whole of 
the resistivity which is measured at 4-2°K. This is the resistivity with 
which this discussion has so far been concerned. There is also the 
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_ resistivity p, caused by the scattering of the electrons by the lattice 
vibrations. This decreases rapidly as the temperature is lowered but at 
20°, in our sodium specimens, it still accounted for about 90% of the 
total resistance of the specimen. Gugan and Dugdale have shown that 
the value of p, for the h.c.p. phase of sodium is considerably less than 
that for the b.c.c. phase at the same temperature. Now as we deform 
sodium at 20°x we transform more of the metal to the h.c.p. phase and 
hence we tend to decrease the value of p;. At the same time, however, 
we introduce dislocations and possibly other types of defect which will 
increase py. By chance these two effects practically cancel one another 
out, hence giving the impression that no resistivity changes are occurring. 
As has already been stated the results actually show that during the first 
few per cent of strain the total resistivity decreases very slightly and 
thereafter it begins to increase again so that an exact cancellation of the 
two effects does not occur. 

When the specimen is now warmed up to 80°K a great deal of this extra 
h.c.p. phase, which was introduced by the deformation, reverts to b.c.c. 
and even when we re-cool to 20°K again, whilst some will transform back 
to h.c.p., there will be very much less of that phase present than before 
we warmed up. Hence the resistivity will increase and this is what is 
observed. The magnitude of this increase, Ap’, will be approximately 
equal to the decrease in the resistivity due to the change to the h.c.p. 
phase during the original deformation. Since as has been explained this 
decrease in resistivity must be approximately equal to the increase in 
resistivity due to dislocations and other defects we can get a measure of 
this increase from the value of Ap’. 

Some support for this argument is obtained when we look at the values 
of Ap’ which were obtained. Three runs in which strains of 0-16, 0-13 
and 0-12 were introduced at 20°K gave values of Ap’ on warming to 80°K 
and then re-cooling to 20°K of 1-1, 2-2 and 3 x 10-® ohm em respectively. 
These figures are of the same order of magnitude as the extra resistivity 
introduced when sodium is strained at 4:2°K by these amounts. At 
4-2°K we assume that the resistivity changes observed are due only to 
the defects produced by the deformation. The fact that one does get 
this rough agreement between these two sets of results suggests that 
a considerable part of the resistivity increase produced by deformation 
at 4-2°k must be due to dislocations rather than to point defects, because 
otherwise Ap’, after warming to 80°K, would be reduced by the annealing 
out of point defects at that temperature. 

Nevertheless, some of the resistivity increase might be ascribed to point 
defects which are produced during the deformation and in this connection 
we should point out that when copper is deformed at low temperatures 
the corresponding values for the parameters A and n are about 5-5 x 10-8 
and 1-3 to 1-5 respectively, i.e. Ap is about twice that observed in our 
sodium specimens. Of this resistivity change about half can be ascribed 
to point defect scattering (see for example van Bueren 1955). This gives 
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the order of magnitude which might be expected but it is obvious that the 
complications of the phase changes in the present case are such that any 
more detailed comparison would be pointless. 

The resistivity changes in lithium at 20°K were about 30% less than 
those at 4:2°K. This reduction might also be due to a decrease in the 
value of p; for the transformed phases of lithium. Since the Debye @ of 
lithium is much higher than that of sodium only a few per cent of the 
resistivity at 20°K will be due to p;, but it is conceivable that the change 
in p; might be large enough to give the difference between the results at 
20 and 4:2°K which is observed. 
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ABSTRACT 


An electron diffraction study has been made of the growth of evaporated 
silver layers on cleaved mica surfaces, with the substrate at various tempera- 
tures during the deposition. It is shown that good orientation (epitaxy) 
occurs over a wide range of temperatures, but that the nature of the orienta- 
tion depends upon both the substrate temperature and the film thickness. 
Under certain conditions of preparation it is found that the surface of the silver 
deposit is very smooth on an atomic scale. These surfaces are of value as sub- 
strates for other experiments, both on epitaxy and on specimen preparation for 
high resolution electron microscopy of metals, and examples of these applica- 
tions are quoted. 


§ 1. INTRODUCTION 


DETAILED studies of the oriented growth of thin surface deposits in single 
crystal substrates (epitaxy) require the use of substrate surfaces which are 
very smooth on an atomic scale. Systematic work (for a general review see 
Pashley (1956)) has been carried out with cleavage surfaces as substrates, 
but when metal single crystal substrates are required, cleavage is rarely 
possible. Electropolished surfaces may be prepared on metallic single 
crystals, but these are usually quite undulating (Kranert et al. 1944), and 
not always suitable as substrates. Several years ago the author developed 
a method of preparing good smooth single crystal surfaces of silver by 
evaporation on to mica, following upon earlier work by Rudiger (1937) 
and Tull (1951). These surfaces have since been used for a number of 
different experiments on epitaxy (Newman and Pashley 1955, Kehoe et al. 
1956a, b, Newman 1957, Mason 1958, Grunbaum 1958, Grunbaum et al. 
1958). Further, the silver substrates have also been used as the first stage 
in the preparation of detached single crystal films of gold and other metals 
for use as specimens in high resolution microscopy (Pashley 1959, Pashley 
et al. 1957, Bassett et al. 1958). : 
Since no detailed account of the technique of preparing the silver surfaces 
has been published to date, the investigations which resulted in the deter- 
mination of the optimum preparation conditions are presented in this 
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paper. The work consists of an investigation by reflection electron diffrac- 
tion of the structure of layers of silver evaporated on to a mica cleavage 
surface under a variety of conditions. Rudiger (1937) was the first to show 
that a good single crystal orientation of silver may be obtained on mica, 
and Tull (1951) prepared such layers which were also quite smooth on an 
atomic scale (see, however, § 4). 


§ 2. THE Evaporation TECHNIQUE 


Preliminary investigations were carried out by heating pieces of freshly 
cleaved mica inside a small porcelain tube surrounded by an electric heater. 
and depositing various thicknesses of silver at different substrate tempera- 
tures, inside a standard type of evaporating plant. It was immediately 
established that good orientation of the silver could be obtained, and that 
the perfection of the orientation and the smoothness of the substrate surface 
depended upon substrate temperature, rate of deposition and deposit thick - 
ness. An arrangement which allowed these parameters to be controlled 
was therefore constructed. 

The mica substrates were placed on a hotplate consisting of a composite 
cylindrical copper block, heated by an internal nichrome winding, and 
supported on three quartz pillars. The temperature of the hotplate was 
measured by a copper-constantan thermocouple embedded in its surface. 
Once thermal equilibrium is established, the temperature of the mica 
surface should be within one degree centigrade of this temperature. The 
silver was deposited from a heated tungsten spiral. With this arrangement 
it was found initially that only very imperfect orientation could be obtained, 
and this became worse as the substrate temperature was raised. It was 
established that the cause of the discrepancy arose from the distillation of 
low melting point impurity metals from the copper block, which was of 
commercial impurity only. Mica specimens were heated on the hotplate 
and then examined by reflection electron diffraction without any silver 
layer being deposited. It was found that small amounts (probably less 
than 104 in average thickness) of two substances had deposited on the 
mica surface. These substances gave bright and clear diffraction patterns 
consistent with the two alloys SnAs and SnTe, which have structures of the 
rocksalt type and lattice constants of 5-68 A and 6-29 A respectively. Both 
alloys were very well oriented as follows : 


(111) alloy parallel to mica (001), 
with [110] and [110] alloy parallel to mica [100]. 


The occurrence of the metals Sn, As and Te as impurities seems quite likely, 
since such metals are sometimes added to improve machining properties. 
Clearly the presence of these thin alloy deposits was affecting the subsequent 
growth of the silver deposits. (To overcome this difficulty, a second hotplate 
was constructed of ‘specpure’ copper, with a nominal purity of better than 
99-999%. This gave results which were always consistent, good orienta- 
tion being obtained under the appropriate conditions. Attempts . to 
construct a furnace of electrolytic copper resulted in the formation of films 
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with poor orientation. Although no impurity deposits could be detected in 
this case, itis presumed that very small amounts did occur, and influence the 
growth of the silver deposit. 

Tt was found also that even with the furnace constructed of specpure 
copper, the surface smoothness could not be reliably controlled unless certain 
further precautions were taken. It was necessary to minimize the effect of 
the radiant heat from the tungsten spiral, which caused rises in the surface 
temperature of the silver film during preparation. ‘This was carried out by 
using suitable shields, the final arrangement being shown in fig. 1. The 
evaporating distance was about 12cm and the vacuum was better than 
10-4mmHg. This gave good reproducible results in all respects. The 
specimens were always prepared in a separate evaporating plant and 
transferred to an electron diffraction camera for examination. 


Fig. 1 
E 


55. te gate SIS oy ame D 


The arrangement for specimen preparation. Ais the furnace; B is the specimen; 
C is a cylindrical aluminium shield coated internally with graphite; Dis 
a flat shield with an aperture through which the silver passes; E is the 
tungsten spiral from which the silver is evaporated. 


§ 3. EXPERIMENTAL OBSERVATIONS 
3.1. The Effect of Substrate Temperature 


Preliminary experiments indicated that thicknesses greater than a 
few hundred angstréms are required for good surface smoothness and 
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coherence of deposit. Thus a series of specimens of about 10004 in 
thickness was prepared, at various substrate temperatures from 20°C to 
500°c. The results are summarized in the table, and illustrated by the 
photos of figs. 2 to 4 (Pl. 48), the interpretation of which is indicated in 
fig. (5). The orientation referred to as single positioning is 


(111) Ag parallel to (001) mica cleavage, | a) 


[110] Ag parallel to [010] mica. 


Furnace | Random | Single Double Degree of 
temp. orienta- position- position- ao eae surface 
(°c) tion ing ing ng smoothness | 
20 S N S SE Na N 
40 M N S aa ee N 
60 Ww N S Sie aly w 
75 VW N S S Small M 
90 N N S Ww Small M 
250-300 N s vw N ie Complete 
500 N = S Saapgie N 


S=strong, M=medium, W=weak, VW=very weak, N=none. 
For double positioning, the following orientation, occurs also. 


(111) Ag parallel to (001) mica cleavage, | (2) 


[110] Ag parallel to [010] mica. 


Twinning was prominent in many cases, and occurred on all of the {111} 
planes. Asthesubstrate temperature is increased the orientation gradually 
improves, the twinning becomes less prominent, the crystal size increases. 
and the surface becomes smoother, until at 250-300°c a good single crystal 
orientation without twinning or double positioning occurs, and the surface 
of the film is very smooth. Above the substrate temperature of 300°C 
the surface structure gradually changes again. In particular the surface 
is no longer smooth, and prominent twinning reappears. 
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3.2. The Effect of Film Thickness 


The influence of the film thickness on the structure of the deposit has 
been investigated with the substrate temperature at its optimum value of 
about 275°c. In the very early stages the deposit consists of small isolated 
nuclei of a few tens of angstréms in linear dimensions. Many are oriented 
in the usual way, but some are randomly oriented, and prominent twinning 
on all {111} planes occurs also. Figure 6 (Pl. 43) shows the diffraction 
pattern obtained from a deposit of 50A. Diffuse rings appear also at 


Fig. 5 
e ap Pitas ® 
seth + + 6 
oan @ t 
7; e tee 
a @ fe) ae 
@ Co) ay a, re) 
@ ® e 
ny Oa 
ur se .@) ovat 
a Mee 
: er z 
ry 


e SPOTS DUE TO SINGLE POSITIONING. [ITO] AZIMUTH 


© EXTRA SPOTS DUE TO DOUBLE POSITIONING 
+SPOTS DUE TO {III} TWINNING 


The interpretation of figs. 2, 3, 4 and 6. 


earlier stages. As the thickness is increased, the orientation becomes more 
perfect and the rings disappear from the diffraction pattern. The spots 
sharpen laterally, and elongate perpendicular to the shadow edge. Ata 
thickness of about 2004 the diffraction pattern is similar to that of fig. 3. 
Further increases in thickness produce a continuation of these trends 
together with a gradual weakening of the diffraction spots due to twinning, 
The twinning at the surface disappears completely at a thickness of about 
500A. Ata thickness of 1000 4 the smooth well-oriented surface described 
above is obtained. 
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A further technique which has been used to study the silver films is that 
of high-angle Kikuchi lines (Alam et al. 1954). The high-angle Kikuchi 
patterns may be obtained with a much larger glancing angle of incidence 
for the electrons, and the depth of penetration is correspondingly larger. 
By means of the technique the single and double positioning are more 
easily distinguished than from a pattern such as that of fig. 4. It was 
confirmed that with a substrate temperature of 250-300°c complete single 
positioning does occur at a thickness of 10004. However, even with 
thinner deposits, which were partly doubly positioned according to the 
conventional reflection diffraction pattern, prominent single positioning 
was indicated by the high-angle Kikuchi patterns. This demonstrates 
that the double positioning occurs mainly in the surface layers, and not 
throughout the depth of the deposit. Figure 7 (PI. 44) shows a high-angle 
Kikuchi pattern from a 10004 layer deposited at 270°c. In addition to 
indicating the complete single positioning of the silver, the sharpness of 
the pattern also indicates the high degree of perfection of the orientation. 


3.3. The Rate of Deposition 


Since this work is concerned with the preparation of films of at least 
several hundred angstréms in thickness, slow rates of deposition lead to a 
long time of evaporation, with the result that large increases in the surface 
temperature of the silver surface occur due to radiant heat from the 
evaporation source. This affects, in particular, the surface smoothness 
of the deposit. A deposition rate of approximately 500 A/min has therefore 
been used in all this work. 


$4, SumMARY AND DISCUSSION 


It has been established that over a wide temperature range silver is 
well oriented on a mica cleavage surface in the orientations given in §3. 
This is consistent with the work of Tull (1951), but does not agree well 
with that of Rudiger (1937) and Thirsk (1950), who found other orientations 
also. In view of the stringent conditions found in the present experiments 
(see §2) concerning the need for cleanliness of the substrate during the 
deposition, it is concluded that the results of Rudiger and Thirsk were 
affected by slight traces of impurity. The consistent results obtained for 
the structure of the deposit grown as described in §2, in at least three 
different pieces of apparatus over a period of several years, supports the 
view that the orientation described above is the one characteristic of a 
clean mica substrate. This orientation has a fairly good fit, as illustrated 
in fig. 8. Along the mica [100] the misfit is — 4%, and is the same along 
the mica [010], if a fit of three silver atom distances with the corresponding 
mica distance is considered. 

The surfaces obtained with a substrate temperature of about 270°C 
reflect as a good optical mirror, in contrast with those obtained by Tull 
(1951), which were somewhat dull in appearance. The exact conditions 
under which Tull’s specimens were prepared were not specified, so that 
no direct correlationis possible. Itis deduced that a substrate temperature 


P.M. x 
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of greater than 300°c was used, Although the surfaces were not optically 
smooth, they appeared smooth according to the electron diffraction 
evidence, as indicated by the extent to which refraction of the electrons 
occurred at the surface. Tull was able to measure the inner potential for 
silver as 23 volts. The main requirement for this type of measurement is 
that the surface is atomically smooth over areas of linear dimension of 
about a few hundred angstroms. Thus irregularities such as steps, pits, 
etc, may be present, but must be greater than about 10004 apart. Such 


Fig. 8 
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The fit of the main silver orientation on mica. 


surface irregularities could, of course, affect the optical reflectivity. It is 
deduced that the surfaces obtained by depositing 10004 of silver at a 
substrate temperature of 250-300°c as described in § 2 are more generally 
smooth than those prepared by Tull. No accurate inner potential measure- 
ments were made on these specimens, but it was confirmed that the 
refractive effect was consistent with a value of about 20 volts. The 
surfaces are therefore considered to be smooth to near atomic dimensions 
over areas of at least about 10004 in linear dimension. Various kinds of 
irregularity might be present within this limitation, but these are not 
sufficient to affect the optical reflectivity appreciably. 
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The value of the (111) silver surfaces as substrates for other experiments 
has been demonstrated by the applications quoted in §1. In addition to 
these applications, the author has studied the growth of silver bromide 
and silver chloride on the silver surfaces, by direct chemical attack with 
halogen vapour. The orientations found are those expected for growth 
on a smooth (111) surface, as previously determined (Pashley 1952), and 
no evidence for any growth on, a facet oblique to the surface was obtained. 
This result was checked with the high-angle Kikuchi patterns, as well as 
by conventional reflection diffraction. Thus there was no possibility of 
any appreciable amount of halide growing in another orientation, and in a 
place on the surface not accessible for diffraction by the conventional 
technique. From this evidence it is deduced that no facets of appreciable 
size and number exist oblique to the surface. 

It is interesting to note that the perfection of orientation improves as. 
the thickness of the film is increased. Two explanations are possible. 
Hither the well-oriented nuclei grow faster, and gradually predominate as 
the film is made thicker, or the well-oriented nuclei grow at the expense of 
the badly oriented nuclei (including randomly oriented and twinned 
nuclei), so that a ‘recrystallization ’ of the film occurs as it is made thicker. 
The possibility of the latter mechanism is supported by the fact that if a 
10004 film is deposited at room temperature, when it is very poorly 
oriented (see table), a subsequent anneal at 300°C produces a well- 
oriented film. 
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ABSTRACT 


An evaporation technique has been developed to allow thin (approx. 200 A 
thick) detached single crystal films of gold to be prepared, ‘These films are 
very coherent and uniform in thickness. Dislocation lines are observed in 
these films with a density of 10!°-10"/cm?. The dislocations may be observed 
to move inside the electron microscope. A prominent feature of the 
movement is the occurrence of wide ribbons of stacking fault due to the 
splitting of a dislocation in to two partial dislocations. Cross-slip is observed 


also. 

Annealing of the films has produced some regular arrays of dislocations, 
including hexagonal networks and arrays at low angle boundaries. 

Small loops are also found in some of the films and are interpreted as dis- 
location loops resulting from the collapse of discs of vacancies. 


§ 1. IyTRODUCTION 


THE application of transmission electron microscopy to the direct study 
of the microstructure of metals is largely limited by the techniques available 
for preparing suitable thin specimens. Various methods have been 
employed. Beaten gold foil (approx. 10004 thick) may be examined 
directly or after further thinning by a chemical etchant (Hirsch et al. 1955). 
Chemical thinning of a thicker beaten foil (0-5 1 aluminium) is also possible 
(Hirsch e¢ al. 1956). Heidenreich (1949) used an electro-polishing tech- 
nique for thinning aluminium, and this has been extended and improved 
by Bollmann (1956) and Tomlinson (1958) to allow a variety of metals to 
be thinned. Nicholson et al. (1958) have used the electro-polishing 
technique to thin aluminium alloys. 

Thinning can also be carried out by ionic bombardment (Castaing 1955) 
and thin sections of metals may be cut by means of an ultra-microtome 
(Haanstra 1955, Fernandez-Moran 1956). Thin films which are suitable 
for transmission electron microscopy can be prepared directly by electro- 
deposition (Weil and Read 1950, Reimer 1956). 

All of these methods provide specimens which are either non-uniform 
in thickness locally, or which are of a uniform thickness over only quite 
small areas. Since the author has been concerned with carrying out 
experiments on controlled deformation of metal films inside the electron 
microscope (this work will be reported in a later paper), it seemed desirable 
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to produce films which are uniform in thickness over reasonably large 
areas. For this reason the preparation of thin metallic films by means of 
evaporation techniques has been investigated. A further advantage of 
the evaporation technique is that it is possible, by using single crystal sub- 
strates, to produce oriented ‘single crystal’ films. 

Often, evaporated films of metals of only a few hundred angstroms or 
less in thickness are quite granular in structure, and not very continuous 
or homogeneous. An improvement in continuity may be obtained by 
the use of thicker films, but this is often undesirable because the electron 
transmission is then poor. Most attempts at making thin detached films 
by evaporation techniques (see for example Raether 1952, Ogawa and 
Watanabe 1954) have involved the use of non-metallic substrates, 
particularly alkali halides. This paper describes a technique involving 
the use of metallic substrates which are themselves produced by evaporation 
(Pashley 1959). Good coherent single crystal films of gold are obtained, 
and dislocations and their movement are observed in these films. 


§ 2. THE PREPARATION AND STRUCTURE OF THE FILMS 


The main difficulty in the use of metallic substrates concerns the surface 
preparation. Electro-polishing is a possible method, but this normally 
produces a wavy surface, and it would be difficult to avoid small amounts. 
of local surface contamination which might cause gross imperfections in 
the films (see, however, Brame and Evans 1958). Itseemed more promising 
to prepare the initial surface by an evaporation technique, since this may 
be carried out conveniently for the case of a (111) silver surface (see previous 
paper). Such a surface is very smooth on an atomic scale, and is an ideal 
substrate for gold since the misfit for parallel orientation is only 0-2%. 
Also, the gold film is readily detached by dissolving the silver in nitric acid. 

It is found that (111) gold films of 1504 and upwards in thickness may 
be prepared in this way. The deposition of the silver and gold is carried 
out in quick succession by means of two evaporation sources built in to a 
standard evaporation unit. The substrate is heated to about 270°C 
during both depositions, by being placed on a small hotplate constructed 
of ‘specpure’ silver with mica as an insulating material and pure molyb- 
denum wire asa heater. The results described in this paper were obtained 
with gold films of approximately 2004 in thickness. The films are good 
single crystals over large areas (many mm?), as indicated by transmission 
electron diffraction patterns. Examination in the electron microscope 
(Siemens Elmiskop I) indicates that the films are coherent and non- 
granular. When the films are stripped and mounted on a supporting 
grid by lifting from the surface of water, local buckling occurs. This leads. 
to prominent extinction contours on the image. 

An extreme case of local bending is shown in fig. 1 (Pl. 45), where complex 
contours are visible. Similar contours have been obtained previously 
with lead iodide (Mollenstedt 1953) and mica (Rang 1953a, b) and have 
been attributed to the occurrence of bubbles in the films, as shown in 
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fig. 2(a). With the gold films the effect may be explained as due to the 
presence of small ‘dimpled regions’ in the film, as illustrated in fig. 2 (0). 
These would arise if the silver and gold is deposited on the surface of mica 
which contains the bubble-like features shown in fig. 2(a). It is assumed 
that the lattice is elastically strained so that the [111] axis is always 
perpendicular to the local plane of the film. The contours obtained 
represent the loci of points in the film which are diffracting in the same way ; 
i.e. regions for which the same reciprocal lattice point is intersected by the 
sphere of reflection. On this basis, the effect of a dimpled region approxi- 
mated to the shape H = Ae~*” (see fig. 3(a)) has been calculated, and the 
resulting contour pattern is shown in fig. 3 (6), for the case of the electron 
beam directed along the normal to the general film plate. The parameters 


Fig. 2 
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Cb 
(a) The type of bubble observed in a thin film of lead (Mollenstedt 1953) and in 
mica (Rang 1953 a, b). 
(6) Dimpled region occurring in a gold film. 


A and a were chosen such that the maximum slope dH/dr is 0-08 radians 
at a value of r=0:16 4 (see fig. 3(a)). The radius of the dimpled region 
is then a few tenths of a micron, and the value of H at r=0 is 2104. The 
resulting contour pattern shown in fig. 3(b) is symmetrical, whereas the 
features shown in fig. 1 are not. This is because the beam is not exactly 
perpendicular to the general film plane. If a tilt of only 0-05 radians is 
allowed for, the contour pattern takes the form of fig. 3(c). This is 
asymmetrical in a manner similar to that of the features shown in fig. 1. 
The calculated example is intended only as an illustration of the main 
features of the contour patterns from the dimpled regions, particularly 
the effect of tilting. No attempt has been made to calculate in any detail 


100 (a) 


; | 230) \co22) a) 


(a) The shape of the dimpled region used for calculating the extinction contour 
P p g g 
patterns of figs. 3 (b) and (c). 
(b) The {202} extinction contours for the beam perpendicular to the general 
film plane. 
(c) The {202} extinction contours for the beam rotated by 0-05 radians about 
[112] with respect to the orientation of (6). 
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the exact shape of the film corresponding to the contours of fig. 1. A 
method for doing this has been given by Rang (1953 a). 

The use of mica as the initial substrate provides single crystal gold 
films with their (111) planes parallel to the film plane. Attempts have been 
made to produce films in (100) orientation, by using rocksalt in place of 
the mica, and retaining the intermediate silver layer of about 20004. 
This technique produces gold films which are uniform, continuous and 
well oriented, but they contain many grown-in twins, as shown both by 
transmission diffraction and by microscopy (see fig. 4 (PI. 45)). Occasionally 
areas of 0-1mm dimension may be prepared free of twins, but no control 
of this behaviour has been obtained. 


§ 3. THe OBSERVATION OF DISLOCATIONS AND THEIR MOVEMENT 


Dislocations in the gold films may be observed by diffraction contrast, 
following the technique first used by Hirsch e¢ al. (1956) and Bollmann 
(1956). The dislocations are formed in the specimens during the growth 
process. All the observed dislocations lie in the (111), (111) and (111) 
slip planes which make angles of 70-5° with the plane of the film, so that 
the projected width of the slip plane projected on to (111) is given by 
tcotan 70-5°, where t is the thickness of the film. With t~ 200A, the 
projected width is only about 704. The (111), (111) and (111) slip planes 
intersect (111) along [110], [101] and [011], which are three directions at 
- 60° to one another. The dislocations pass from top to bottom of the film, 
and tend to assume their shortest possible length within the slip plane. 
Their distribution throughout the area of the film appears to be random. 
Dislocation lines running approximately parallel to the film plane are not 
observed. The projected length of the dislocation lines is very short 
(also about 704), and in most cases the dislocations appear as dots on the 
electron microscope image, usually black (on a positive print), sometimes 
white and sometimes black-white in character. The white dots appear 
only in regions which are diffracting strongly, i.e. in an extinction contour. 
The visibility or contrast of a dislocation varies with the general diffraction 
conditions in the immediate vicinity of the dislocation, and is often quite 
low. ‘The dislocations in the (111) films are more readily detected as they 
move (see below). 

Dislocations are more easy to detect unambiguously in the (100) gold 
films, partly because of the more favourable geometry. The four {111} 
planes make angles of 54-8° with the (100) film plane, and the projected 
width of the slip plane for a film of 2004 in thickness is approximately 
140A. The dislocations appear on the images as short lines, with pro- 
nounced intensity maxima at the two ends (see fig. 4). Because of the 
greater ease of detection, it is possible to make a reliable estimate of the 
dislocation density in the (100) films. A value of 101°_104/em?. is 
obtained, in good agreement with the estimates of 10!/cm? made on the 
(111) gold films by means of moiré patterns (Bassett et al. 1958). 
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In addition to the dislocations described above, a number of the (111) 
specimens has contained a high density of small loops (see fig. 5 (PI. 45)). 
These are very similar to the dislocation loops recently observed by Hirsch 
et al. (1958) in quenched aluminium, and often have the appearance of 
double lines (see A, fig. 5) as frequently observed by Hirsch et al. (see their 
§3.1). The loops are always less than 1004 in diameter, the smallest 
distinguishable loop being about 304 across. No control of the pre- 
paration of specimens containing the loops has been obtained. Because 
the appearance of the loops is very similar to that of the loops observed by 
Hirsch et al., it is assumed that they are dislocation loops of the same kind, 
i.e. they arise from collapsed dises of vacancies in which a Shockley partial 
is nucleated so that a prismatic $[110] dislocation is produced (Read 1953). 
Presumably large numbers of vacancies become built into the lattice during 
the deposition process. In these circumstances the formation of the loops 
might well depend upon the rates of deposition at various stages of the 
growth. This factor is not well controlled in the present experiments, 
and requires further investigation. 

The dislocations observed in the gold films may be caused to move by 
use of an intense and localized illuminating electron beam, as obtained. 
with the double condenser lens of the Siemens Elmiskop I. Dislocation 
movement induced in this way has already been observed by Hirsh e¢ al. 
(1956, 1957) with aluminium films and Whelan eé al. (1957) with films of 
stainless steel. The detailed mechanism which causes the motion is not 
clear, but it is assumed that the main effect arises from the compressive 
stresses which occur in the plane of the film, as a result of the extremely 
localized heating by the incident electrons. There is a remarkable 
influence of the time of electron bombardment on the behaviour of the 
specimen. When an area is first examined at high local illumination, no 
marked dislocation movement is observed, until that particular area has. 
been illuminated for a period of about one minute, after which pronounced 
dislocation movement is usually observed. This induction period does not 
arise because a finite time is required for a certain temperature to be 
reached, since the illumination may be removed from the particular area 
for a period of two or three minutes, and pronounced dislocation move- 
ment occurs immediately the illumination is restored. Apparently the 
dislocations are in some way locked in the film as prepared, and intense 
localized electron bombardment results in the unlocking of the dis- 
locations. Alternatively the formation of a carbon contamination film 
may be in some way responsible, although deposition of a thin carbon 
film on each side of the specimen before observation by the method 
described by Bradley (1954) causes no change in the general behaviour of 
the specimen. In both cases, the detailed mechanism is obscure. 

Various kinds of motion of the dislocation lines are observed : 


(1) A wandering and extremely slow (<10~°cm/sec) motion not. 
confined to any particular slip plane; these are presumably screw dis- 
locations. 
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(2) A faster and sudden motion confined to one well-defined slip plane ; 
this motion results in the formation, of a slip trace, as observed by Hirsch 
et al. (1956), which fades after a period of about 10-20 sec. 

(3) Motion as indicated in (2) except that single and multiple cross-slip 
occurs. 

(4) The sudden splitting of a dislocation into two partial dislocations, 
with the production of a wide ribbon of stacking faults. 


In addition, the dislocation loops such as shown in fig. 5 sometimes 
disappear suddenly during observation, and sometimes appear equally 
suddenly. No details of their motion have been observed. Hirsch et al. 
(1958) have given details in the case of the movement of loops in aluminium. 


A stacking fault AB in a thin film. 


The first three types of motion have already been described by Hirsch 
et al. (1956) for aluminium foils, and no further description will be given 
here. The fourth motion will be considered in more detail. This process 
shows itself by the appearance of fine lines on the micrograph, as shown in 
fig. 6 (Pl. 45) and fig. 8 (Pl. 46). The lines appear very suddenly, some 
disappear equally suddenly, and the total number increases as the electron 
bombardment is continued. The lines do not disappear as the electron 
illumination is reduced. 

This type of feature was not properly understood until the theory for 
the contrast produced by a stacking fault was worked out by Whelan and 
Hirsch (1957a, b). According to their theory, if a stacking fault occurs 
on a crystal plane AB (see fig. 7) inclined at an angle ¢ to the plane of the 
film, and if the shear corresponding to the stacking fault has a non-zero 
component perpendicular to the lattice planes which are giving rise to a 
Bragg reflection, the intensity of the Bragg reflection from a column of 
the film such as PQ will depend upon the depth z of the stacking fault 
below the surface. When zis equal to an integral multiple of the extinction 
distance s for a particular Bragg reflection, the diffracted beam from the 
part of the column PQ above the faulting plane is zero. Similar conditions 
apply to the lower part. It is shown that as a result of this effect the 
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intensity of the diffracted beam, and hence also the intensity of the un- 
deviated transmitted beam, depends periodically on z and well-defined 
fringes are observed on the image with a spacing of scotan¢. This 
analysis applies when the crystal is set nearly, but not exactly, at the 
Bragg angle: at the exact setting the fringe spacing is }scotandé. With 
the (111) gold films used in the work described in this paper, the main 
reflections which occur are {220}. The estimation of extinction distance 
involves a knowledge of the atomic scattering factor f, for electrons. 
Values of s have been deduced from the f, values given recently by Ibers 
(1958), and these are given in the table, the value for {220} being about 
250A. This is approximately equal to the thickness of the film, so that 
only one fringe, or at the most two fringes, of the stacking fault contrast 
can be expected in this case. Films of much greater thickness do not 
give rise to the wide ribbons of stacking fault during observation in the 
Elmiskop I, so no confirmation of the detailed contrast can be made in 
this way. With the films of 200-3004 in thickness, the contrast in the 


Extinction Distance (s) for Gold Based upon Atomic Scattering 
Amplitudes for Electrons (f,) as given by Ibers (1958) 


Reflection f, (A) s (A) for accelerating potential of 
(hkl) : 


80 kv 


(111) 7-6 155 


6-7 


4:8 


stacking faults can be changed by tilting the specimens. This is arranged 
to allow {111} and {200} Bragg reflections to occur, so that extinction 
distances of 1554 and 1754 are brought into operation and the path 
length of the electrons passing through the film is increased by a factor 
sec 6 (9=angle of tilt) so that more than one fringe appears; by a suitable 
choice of the tilting direction the projected width of the stacking fault 
may be increased to allow the individual fringes to be resolved more 
readily. This technique is illustrated in fig. 8. The tilting is such that 
one set of stacking faults is on {111} planes steeply inclined to the electron 
beam, and two fringes are observed, whereas the other set of stacking faults 
is on {111} planes not so steeply inclined, and the two fringes are not 
resolved, so that the normal fine lines appear. The stacking faults appear 
also during the observation of the (100) gold films, but these have not been 
studied in such detail. With these films the {002} reflections are pre- 
dominant, and two fringes may be observed without any tilting. 

A further observation is that the stacking fault fringes sometimes have 
a twisted appearance, as shown in fig. 9 (Pl. 46). This effect is predicted 


332 D. W. Pashley on the Observation of Dislocations 


by the theory of Whelan and Hirsch (1957a, b). A final piece of evidence: 
that the observed fine lines are due to stacking faults is that when they 
disappear suddenly, a single dislocation may be observed to appear 
simultaneously at one end of the original line. This is as expected since 
one partial dislocation moves and combines with the other partial to 
produce an unsplit dislocation as the stacking fault is eliminated. 

The width of the stacking faults (i.e. the separation between the two 
partials) may be as high as one micron. This is much larger than the 
estimated equilibrium separation of about 354 for gold (Seeger 1955). 
The reason why such wide stacking faults can be produced seems to be 
associated with the effect of dislocation pinning which occurs at the 
surfaces of the film. Such pinning has been discussed previously (Hirsch 
et al. 1956, Whelan et al. 1957), and is considered to be due to the presence 
of a surface oxide film in some cases. Although no surface oxide is present 
on the gold films, the specimens rapidly become covered by a film of carbon 
contamination during examination in the electron microscope. This. 
contamination is likely to provide an effective means of pinning dis- 
locations by hindering the formation on the surface of the steps associated. 
with the movement of dislocation. On this mechanism the pinning will 
depend upon the Burgers vector of the dislocation, since the step is given 
by the component of the Burgers vector resolved perpendicular to the 
surface. Thus of the two partial dislocations into which a dislocation 
splits, one may be more effectively pinned than the other. Further, the 
stress acting on a dislocation constraining it to move also depends, for a 
given stress system, on its Burgers vector, so that the stress acting on one 
partial may be much higher than that acting on the other. As a result of 
these two effects, one partial may be caused to move whilst the other 
remains stationary, a certain stress being required to keep them apart. 
Since the partials do not automatically come together again once the- 
stress is removed, it is assumed that the carbon contamination film which 
forms after the one partial has moved and formed a surface step, hinders 
the removal of that step, and so prevents the reversal of the movement. 

Ciné films of the movement of dislocations have been taken by direct 
photography of the fluorescent screen, with an instrumental magnification 
of 40000x. It is found that the maximum permissible frame speed is. 
about 10/sec. Under these conditions the stacking faults appear between. 
the two successive exposures, so that they appear in less than 10-'sec. 


§ 4. Tue Errect oF ANNEALING THE GoLp Finms 


Attempts at annealing the gold films have resulted in quite uncontrolled 
behaviour. At about 500°c the films often develop many small holes. 
which gradually produce a lace-work structure. Sometimes this is 
accompanied by a considerable increase in the thickness of the remaining 
parts of the films. A complete recrystallization may occur in this way, 
leading to thick crystals several microns in lateral extent in random 
orientation. These crystals are often too thick to allow adequate electron. 
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transmission for microscopy, but thin regions of 20004 and less may be 
examined at an accelerating potential of 100kv. Complicated dis- 
location structures are then visible; these include many long dislocation 
lines which run approximately parallel to the plane of the film (cf. §3 for 
the arrangement of dislocations in the film as prepared). 

Regular dislocation arrays are also found in these annealed specimens. 
Figure 10 (Pl. 46) shows an array of dislocations at low angle grain 
boundaries. Figure 11 (Pl. 46) shows an array which appears to have 
arisen from the interaction of two sets of dislocations on near neighbouring 
slip planes. If it is assumed that during the annealing process the dis- 
locations on the near neighbouring slip planes climb towards one another 
so as to intersect with each other as shown in fig. 12 (a), and that the dis- 
locations have identical Burgers vectors, the dislocations will rearrange 
themselves as shown in fig. 12(6). This rearrangement is illustrated in 
fig. 12(c). At each intersection the slip plane is divided into four areas 
A, B, C and D; over areas B and C the crystal is slipped by one Burgers 
vector unit, and over area D the crystal is slipped by two units with respect 
to the slip over A. The areas B and C therefore join, and the interaction 
occurs as indicated in fig. 12(6)._ Hexagonal networks are also produced 
in the annealed films. 

Fig. 12 


(a) 


The mechanism by which the dislocation array of fig. 11 is formed; (6) is the 
Burgers vector of the dislocations. 


§ 5. SUMMARY AND DISCUSSION 


It has been shown that good coherent and uniform single crystal films 
of gold may be produced by an evaporation technique, with a thickness 
(approx. 200 4) suitable for obtaining bright transmission electron, micro- 
graphs at high magnification. Preliminary investigations suggest that it 
will be possible to prepare films of platinum, palladium and rhodium in 
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the same way. The development of similar techniques for preparing thin 
films of other metals depends upon finding suitable smooth substrates on 
which the deposit metal may be grown as a coherent film, and from which 
the deposit film may be conveniently detached. 

Dislocations are observed in the gold films, their density being 
10!°10"/sm2. The movement of the dislocations resulting from the use 
of a highly localized intense illuminating electron beam has been observed. 
A prominent feature is the occurrence of wide ribbons of stacking fault, 
due to dislocations splitting into two partials which separate by as much 
as one micron. 

Whelan et al. (1957) have compared the nature of the dislocation move- 
ments induced in aluminium and stainless steel. Extensive cross-slip 
but no occurrence of stacking faults was observed with aluminium; with 
stainless steel no cross slip was observed, but the stacking faults were 
prominent. This was interpreted in terms of the very high stacking fault 
energy for aluminium (approximately 200ergs/em?) compared with that 
for stainless steel (approximately 15-20 ergs/cm?). The gold appears to 
be an intermediate case, since both cross-slip and stacking faults are 
prominent. This is consistent with the intermediate stacking fault 
energy of approximately 40 ergs/cm? estimated by Seeger (1955). 

Finally, it must be emphasized that both the original dislocation arrange- 
ments and the behaviour of the dislocations under any form of applied 
stress are characteristic of thin films, and not bulk materials. On the 
other hand the techniques of thinning bulk materials to obtain the thin 
film gives rise to dislocation arrangements which are to a large extent 
characteristic of the bulk material, although any movement of the dis- 
locations is still essentially a thin film behaviour. 
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ABSTRACT 


A new approximate formula is derived for the cross sections for excitation 
of molecular vibration and rotation by impact of slow electrons. Account is 
taken of the fact that in such impacts the incident electron is moving much 
faster than the nuclei in the molecule while its wave function suffers 
considerable distortion by the molecular field. 


§ 1. LyrRopuUcTION 


CoMPARATIVELY little attention has been paid to the theory of the excitation 
of nuclear motion in molecules by electron impact. All of the calculations 
which have been carried out so far have employed Born’s approximation. 
This is certainly not valid in many cases when the electron wavelength is 
comparable with molecular dimensions, though there may be exceptions 
to this. Examples in which use of the approximation would appear to 
be justified have been given, as in the excitation of rotation due to inter- 
action with the permanent electric dipole or quadrupole moment of a 
molecule (Massey 1932, Gerjuoy and Stein 1955). In most cases of 
practical interest, as in microwave breakdown in a molecular gas, or in 
the Luxemburg effect, the electron wavelength is comparable with 
molecular dimensions and to obtain reliable results a more accurate 
approximation is required. 

In seeking this, it is natural to turn one’s attention to the distorted wave 
method which has been recently used with some success in dealing with 
excitation of atoms by electrons with energies close to the threshold. 
This method replaces the plane waves which represent the wave functions 
of the incident and scattered electron waves by waves respectively distorted 
by the mean field of the atom in its initial and final states. The difficulty 
in applying this directly to excitation of nuclear motion in a molecule is 
that, although the electrons are slow, as judged from the relation of the 
wavelength to molecular dimensions, they are nevertheless moving rapidly 
compared with the nuclei and it would not seem to be a good approximation 
to regard the distortion as due to the mean molecular field averaged over all 
orientations of the molecule. It is the purpose of this note to develop an 
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alternative method in which the distortion is produced by the molecular 
field with a fixed orientation and nuclear separation, the rate of excitation, 
of the nuclear motion being then determined by the rate of change of the 
distortion with orientation and nuclear separation. If the molecular field 
is small the formula obtained reduces to Born’s approximation. 


§ 2. DERIVATION OF THE FORMULA FOR THE INELASTIC CROSS SECTION 


We confine ourselves in the first instances to diatomic molecules. 
Electron exchange effects will be ignored initially and it will be supposed 
throughout that no change in the electronic state of the molecule occurs. 
Such questions as coupling between the molecular rotation and electron 
spin will not be considered. 

We denote the vector separation of the nuclei by R, the coordinates of 
the nm molecular electrons relative to the centre of mass of the nuclei by 


ron S=1,...,n and those of the colliding electrons by r. The wave 
equation for the system can be written 
ft ee ve ne VR°—-#,-Hyt+ Vir 
| m(¥ sm) — 5 OM = ris fees) R) 


+0 (PF R)+ SE |B =O, oe 1) 


where V(r,,,,R) is the potential energy of interaction of the molecular 
electrons with the nuclei and with each other and U(r,r,,,,R) that of 
the colliding electron with the nuclei and the molecular electrons. Z,e 
and Z,e are the nuclear charges, Hy is the energy of the initial state of the 
molecule and #,, the kinetic energy of the incident electron. 

In the absence of the colliding electron the molecular wave function is 
represented with good accuracy by 


D = f6(Fom> R)xo(R), ee oe ean (2) 
where ¢, is the wavefunction for the electronic motion and y, for the 
nuclear rotation and vibration. ¢) and y, are proper functions which 
satisfy the equations 


A, 
| Paton) + Veron R)-+ EEE — eR) | dalton R=0, (8) 
h 
E; Vat Hy— «(B)+ saz | $V 2% 08g [x0=0. (4) 
If we write the solution of (1) in the form 
verre ERI A(R YT). Gat ce ten me eee (5) 
then on substitution in (1) we find, using (2) that 


aN E + aa 24 H+ Hy— (Rk) — O(a R) |. 


A [AV ety + 2erady A. gradn $o]=0. Bie os 14) 
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Multiplying by 4,* and integrating over dr,,,, we have, since ois normalized 
to unity for all R, 


E Ve+ iY + [o*Vi8#drom} +E,+E,—«(R) 


= 000s Pom N4olFandlP Aran |A=0. fea) 


The term in grad ¢, vanishes if fy is real because 


| $y tadp by dry =tgerad n| Se Unites 


This equation for A is satisfied approximately by taking 
ASA VHF tr Ri Rye ee) 
where y; is the wave function for the jth state of molecular motion which 
satisfies (4) with H, replaced by H#, and F’, is the wave function for motion 
of the free electron of the appropriate energy in the field of the molecule, 
when the nuclear separation is R, averaged over the coordinates of the 
molecular electrons, Thus 


vr +| ke FE [Ot 5 Pesan) (Pollo )i- dram | F=0. (9) 
The terms which do not cancel when (8) is substituted in (7) are 
h2 
aye (Xi Vn P; + 2 grad, yp grade Fe.” iss) Se eee) 


We may now regard these terms as giving rise to transitions between states 

with different 7, i.e. between states with different wave functions y(R) for 

nuclear vibration and rotation. Applying the usual perturbation theory 

we have that the number of transitions per second from an initial state 0 
of nuclear motion to a final state j is — 

2Qar h? ; 2 

; lear Xj PF i(X0 Ve? ot 2 gradp fy. gradp Fy) drpdr, (11) 

By appropriate normalization of ’) and F’, this yields for the differential 

cross section in which the electron is scattered in the direction (Z, ~) (the 


incident direction being taken as axis) after producing the appropriate 
excitation 


Iy,(6) sin CAL du 
k;| ™ : Z 2 
= im con | oF R)(xo Velo + 2 gradp Xo. gradp Fy) drpdr, 
Sin CC Ct ee ee eee (Es) 
where Fy ~exp (tky.r) + fo (0, R)r— exp (tkgr) 


F',~ exp (—tk,.r) +f;(7— ©, R)r exp (ik,r). 
with cos © = cos 6 cos € + sin @ sin ¢ cos (6 — p). 
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§ 3. RELATION TO BoRN’s APPROXIMATION 


[ue Pom R)| $o(t sn)? AT si = Ol, R). * 4 S (13) 


If Uo) can be treated as small 


We write 


exp (tko|r—r'|) 


[rer] exp (tky.r’)dr, (14) 


E m : 
F,=exp (ik, .r)— saga | Vole ,R) 
so that 


e m fexp(ik)|r—r’ : ' ' 


Hence to this approximation 
Xj F (xo. Va? fot 2 gradp Xo. gradp Fy) dr,drp 


=— oa [| [x exp (— 2k; .r)(X9 Vg? U 99 + 2 gradg yo. gradz Ugo) 
x es (iky.r)dz,dt,drp. . . . (15) 

We have now that 
fexp (tk; .r) 

so that (15) reduces to 


m “ee F 
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exp (iko| r—r’|) 


dr, = 
Jr—r’| Ty 


pee ox (tke 1) tL 8) 
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x gradp y,;.gradpz U9)dtdrp. . . . (17) 
Now [xitxoVn2Uwdrn= | VnP04*x0) Uqodtp: 


| xi" gradp Xo. gradp Uggdtp= — | P Va-(x;* grad yo) rR. 
so, since 
Vn?(x;* Xo) = 2Vp . (x; VrXo) = Xo VR°x;" “< x; VR Xo: 
(15) becomes 


m i . , ! Vy 2 
a maps exp [0(Ky—k;) . r”] Ugo(R, t’)(X0 Vin xy — x3" Vin°Xo) 
MT Oty. melece ewe LS) 
Now, because y, and x; satisfy equations of the form (4), 


ee 2M Me 
Xo Vn?x;* — x3" Vr X0= 7D (By — E;)Xox;* = a. (lig? — &,?)Xox;"> 


so that, on substitution in (12) 


y k, le * % , : 
I= Fl aa [ | x*x0exp [ik —Ky) 10 e(Rrdredrn] (19) 


which is Born’s approximation. 
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On the other hand, when the approximation (14) is not valid it is no 
longer possible to reduce the integral in (12) to the form (18) with the plane 
waves replaced by the appropriate distorted waves. 


§ 4. CONCLUDING REMARKS 


The application of the formula (12) to specific cases is difficult because 
it requires that the functions Ff), /; should be known with sufficient 
accuracy for their derivatives with respect to the nuclear separations to be 
obtainable to a reasonable degree of approximation. An attempt is being 


made nevertheless to carry out the calculations for molecular hydrogen. 
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ABSTRACT 


The quasi-harmonic model of a solid (a harmonic model which takes only 
partial account of anharmonicity) is examined and a general notation 
established. The theory involves the quantities [e¢;]/"l=(V"/v;) (d%v,/aV”), 
where »; are the spectral frequencies. By evaluating averages of [e,]l”! 
appropriate at high and low temperatures it is shown that Griineisen-type: 
assumptions are reasonable for the inert gas solids. This confirms work 
by Barron on Griineisen’s y (in our notation y= — <[e;]JUJ> ay). 

Griineisen’s equation of state pV=—VU,(V)+yU7, where Up and U,, 
are the thermal and non-thermal contributions to the internal energy U, is 
extended to yield an equation of the form V/«p=V?U,)"(V)+yU7—-y2TCy 
where y’(=<[e;]2]> ay) is a ‘second’ Grimeisen constant. It may be 
estimated from measurements of the variation of compressibility with 
temperature in much the same way as the first Grtimeisen constant y is 
determined from measurements on expansivity, specific heat and compressi- 
bility. 

Experimental results for argon are compared with the theory, assuming a 
Lennard—Jones 6-12 potential. Good agreement is obtained except above: 
40°K where anharmonicity must be considered explicitly. 

Our estimates of y are consistent with those of previous writers; from 
the theory we conclude that for argon y’ has a value of about 10 and this is 
confirmed by a direct use of experimental results on compressibility. 


§ 1. INTRODUCTION 


THE most useful approximate theory of the equation of state of simple 
solids is that associated with Mie and Griineisen. It has been summarized 
by Griineisen (1926) and its main feature is the introduction of an 
empirical constant y which is a measure of the anharmonicity of the 
interatomic forces. From the point of view of lattice dynamics, the 
Mie-Griineisen theory is best understood by using a harmonic model of 
the solid and assuming that the dimensionless derivatives of the 
frequencies of the normal modes with respect to the volume, dIny,/dIn V, 
are the same for each mode. They are the negative of the Griineisen y. 
We wish to consider a generalization of this theory in which not only the 
first but also higher derivatives of the frequencies are considered and are 
replaced by appropriate constants. To explain exactly what we shall do, 
it is best to start from first principles. 
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The thermodynamic properties of a simple mechanical system under 
isotropic stress can be obtained from the Helmholtz free energy (/’) when 
it is expressed as a function of temperature (7') and volume (V). Indeed, 
working from BF'(B =1/k7') one has 


BU =pa(BF)/@p. 
pVp=—Va(pFy/av, 


CL/k= — B°0°(BF)/0p?, | eee tt) 
(V/k)(a/«x)= Vo/aV (Bo/oB — 1)(BF), 
VB(1/x) = V?0(BP)/0V®, j 
where U is the internal energy, « is the expansivity and « is the isothermal 
compressibility. 
In a number of systems one can express BF in the form 
Br =] (Bey (ea ey iio sca. og ee 


or else as a sum of such functions, >;/,(B<;). In particular the latter is the 
case for the harmonic model of a solid; the right-hand side of (2) then 
consists of a pure potential term B®—the corresponding function f(x) 
being a—where ®(V) is the potential energy of the undisplaced lattice and 
a sum of ‘ thermal’ terms >;f(B«,) where f(x) = 4a + log [1 —exp(—2)], and 
«(V)=hv,(V) (the v; are the spectral frequencies of the lattice). 

If eqn. (2) is substituted into eqn. (1) we obtain 


BU =/®, 
BpV = —f™[e], 
C, k= —f%, 
(V[k)(a/x) =f? Le}, a (Ss) 
VB/e=fU[eJA+ fae ]t2, 
where ml — yranflda” (a=Be), 
a [el=(V"/e)(One[aV"). | 
For manipulation of these quantities it is convenient to note that 
(x0/Ax) fr = flr + nfird , 
VoloV = [e]@x0/dux, ns Ne} 
(Va/0V)[e]}™ = [ef + n[e] — [ee]. | 


In order to apply (3) to the harmonic model of a solid it must, in accord- 
ance with what was said above, be summed. We therefore need: (i) the 
potential energy of the lattice (®) and its first two derivatives ®" and 
1; (ii) the energies of the normal modes and the distributions of [¢]“ 
and [e]"! over the modes. The thermal f-functions are known. and can 
be obtained by successive differentiation of 


f(x) = 32 + log [1 —exp (—2)]. 


Hence the final result comes from carrying out appropriate summations. 


Griineisen’s Theory of Solids 343 


This programme is not easy to carry through and for numerical 
work it is desirable to refer all the lattice calculations to a reference 
volume (V,) and to expand the thermodynamic functions in powers of 
the dilatation 8=(V—V,)/V». To the second power of 8, the thermo- 
dynamic coefficients may be found from A,(8) in the expansion 


UF ks. BF = A,+ A,8+}A,8°+44,5+54,84, . . . . (5) 
Ay=f, 

A,=f™ <P, 

Ag=f™[e]}2)+ felpe]ae, (6) 


Ag =f™ Le} +23 Le}MLe]+/1[ <P 
Ag =f Le} +f2U(4Le]M Le]! + 3Le}2) +f [e}™L J + f[ eM, 


all quantities on the right of (6) being evaluated at the reference volume 
Ve 

In the simplest form of Griineisen’s theory of the solid state one works 
from eqn. (3), treating [¢] and [e]®! as constant in all the thermal terms. 
([e]™ is the negative of Griineisen’s y and [e]!—which appears, for 
example, in Born and Huang (1954), eqns. (4.55) and (17.15)—can be 
regarded as a second Griineisen constant but it is not usually considered 
in detail.) The potential energy terms must be dealt with separately ; 
for them f(x)=2 and fl =0(n> 2). 

We propose to test on solid argon a simple generalization of Griineisen’s 
theory : we work from eqn. (6) and replace the coefficients of the 
jf-functions in all the thermal terms, viz. [«]™, [e]?!, [e]4?, [e]@. [e], ete. 
by appropriately chosen constant average values (‘ generalized Griineisen 
constants’); the potential energy term is dealt with separately and 
without further approximation. 

Granted the (unjustified) use of a harmonic model, the only arbitrary 
element in our procedure is the choices of averages in forming the generalized 
Griineisen constants. The choices we make here are explained in § 3 and 
are suggested by the work of Domb and Salter (1952), Salter (1954), 
Barron (1955), Barron and Domb (1955). In carrying out the calculations 
we have drawn extensively on numerical work by the same authors. We 
shall work with a Lennard—Jones potential (chosen in § 2) with no further 
fitting of constants. (It is of course possible that another choice of inter- 
molecular potential would yield different results.) The summed ‘/f-func- 
tions’ required are obtained with the help of the frequency distribution of 
Leighton (1948) (§ 4). The final results are compared in §5 with the 
experiments of Clusius (1936), Barker and Dobbs (1955), Dobbs eé al. 
(1956), Stewart (1955, 1956) and with previous calculations by Kane 
(1939), Henkel (1955), Domb and Zucker (1956), Zucker (1956, 1958). 
An experimental estimate of the second Griineisen constant ([e]°'= y’) 
is obtained in § 6. The usefulness of the model is assessed in § 7. 
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§ 2. Tue Iyreratomic PorentiaL AnD Lattice ENERGY OF ARGON 
To fit the two constants in the Lennard—Jones potential 


f(r?) = 4e[(a/r)* — (a/r)*] ox Se ae aan 
we have adopted the experimental values cited by Dobbs and Jones (1957), 
dy = 5°31 40-014, (3) 

— U,=L,= 1850 + 12 cal (g mole), 


for the lattice constant and the latent heat at absolute zero and zero 
pressure. The harmonic theory of Salter was recalculated, taking all 
neighbours into account in the zero point energy as well as in the potential 
term, and it was found that the zero point energy and pressure are given by 


Uy _ 6:0659 _ 14-454 3:295A* pip 
Ne y*4 y*2 Vis ; 

9 

24-264 28-908 4-394 (3-626 — V*?) | 9) 


~ 76 ys 0 Ri2 ? 


(R= 2-072 — V*2) 


where V* = 2(a)/2c)3, A* =h/o(me)'?. (The numbers here differ by about 
4°% from those of Salter.) The elastic moduli at absolute zero were also 
recalculated. 

Substituting the experimental values of (8) into (9) we find for « and o 
the values given in the first entry of table 1. These values are used in 
the remainder of this paper. 


Table 1. The Interatomic Potential of Argon ¢(r?) = 4e[(o/r)!? — (a/r)®} 


Harmonic calculation, all neighbours 
(this work) 
1-69 3-402 Anharmonic calculation, all neighbours 
(Zucker 1956) 
1-65 3-405 Virial coefficient (Hirschfelder et al. 1954) 


Table 2. Values of ®M= V"9rM/dV"” in k cal mole! Obtained from 
Eqn. (9) for the Reference Volume at Absolute Zero (V* =0-9514) 


ol?) ou) Ol?) ols) ol4) 


— 2-026 +0-5489 + 12°37 —119-6 +1018 


The corresponding values of V* and A* are 0-9514 and 0-1865 
respectively; they lead to a zero-point thermal energy of 176 cal. mole-, 
the potential energy of the lattice at absolute zero being 2026 cal mole! and 
the volume at absolute zero being 22-57 cm? mole-. 
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The first two terms of eqn. (9) give the lattice energy (®) and from 
them the successive derivatives of table 2 can be found. The last term 
gives the zero point thermal energy (U,) which is needed § 3. 


§ 3. CALCULATION OF GENERALIZED GRUNEISEN CONSTANTS 


Since the harmonic model which we are using is most likely to be accurate 
at low temperatures, we shall evaluate the averages of [e]", [e]!,. . . ete. 
at low temperatures, and assume that these are temperature independent. 

From the form of the thermal function, f(a) =4$a+log [1—exp (—~2)], 
and the definition, f'=2"0"//dx", it is easy to verify that sums of the 
type >, f'(x,)(v@ =Be,;) have two kinds of behaviour for 7'>0 (8—> oo): 
(A) for n= 1 the sum is proportional to 8 and the high frequency part of the 
spectrum is important; (B) At low enough temperatures the frequency 
spectrum in any element of solid angle is proportional to v?. Therefore 
for n>2 the sum is proportional to 8 and only the low frequency part 
of the spectrum is important. 

The ‘type A’ sums—involving f"!/—are more difficult to deal with and 
such progress as has been made appears in the work on zero point energy 
due to Domb and Salter (1952). Luckily we can draw on their results 
for our purpose, as we shall now show. 


3.1. “Type A’ Averages 


For any quantity 7, distributed over the vibrational modes of the 
crystal we define ‘type A’ averages—denoted by a bar—as 


y= lt {2B [Res malice 


Now by egn. (3) the zero point thermal energy is 
eat > po (Bes Ol wee etn se) , = ots (hl) 
Boo 


Differentiate eqn. (11) successively with respect to volume and make use 
of the fact—noted above— that for B > co sums involving f"™! are of order 
8 whereas those involving f'™ (n> 2) are of order B-°. We then find that 


yn oe = Lt (SfMBedlqIB, . - - » (2) 
ov” Bow 
Finally, combine eqns. (10), (11) and (12), and we have 
——  Yrgru, 
DO eee Oe oS es 8 ee te a LO 
tel = 0, ov" ois 


Equation (13) allows us to calculate the ‘type A’ average values which 
we now choose for the first terms on the right side of eqn. (6). To do 
this we differentiate the expression for U, obtainable from eqn. (9). The 
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results are given in table 3. Since the Domb—Salter approximation 
for U, depends on the second moment of the spectrum, ‘type A’ averages 
are likely to be better at high temperatures. 


Table 3. ‘Type A’ Averages of eqn. (10) Calculated from eqns. (13) and (9) 


feJ41= —3-109 feJ@= + 10-02 [e]3]= —45-36 [e]fI= + 234-3 


‘Type B’ Averages 


We turn now to the other terms on the right of eqn. (6). For these 
we adopt a method due to Born (1923) and used by Barron. He has 
observed that the appropriate low-temperature averages involving f'"(n> 2) 
for a quantity 0,4 depending on the acoustic mode and on the direction of 
propagation of the long waves are defined by 


(6) = ff Giae- 3 dQ/ Ife =8)5 dO ee an) 


where v,¢ is the meloclty, of the ith acoustic mode (7 = 1, 2, 3) in the direction 
Q. We call these the ‘type B’ (B for Born) averages. 

Taking into account only nearest neighbours, it may be shown that the 
velocities of long waves are given by 


5% 2=76"(2)b, +6), . . . . . . (15) 
where 6 is an eigenvalue of 
bs = $075 + COS &, COS a,(1 —46,.), ine SS MDS) 
cos «, being the direction cosines of the propagation vector. It follows 
that : 2 2 
fej | ee Oe eee 


6,OV 4; dr 
If we choose as reference volume that corresponding to minimum inter- 
molecular potentialt, then ¢’(r?)=0 and (15) becomes 


M (p2 
[e|u= a CZ ) +1410} har SATS) 
This quantity and its “type B’ average has been found by Barron. By 
continuing the differentiation of (17) one can form expressions for every 
combination of e« appearing on the right side of eqn. (6); the formulae are 
analogous to (18) but more complex. They are given in Appendix I. 
The main point is that they involve derivatives of the potential of the type 


= rn) f(n)( 72) 
Tes fb" (r?) 
evaluated where ¢'(r?)=0. These are listed in table 4 for the potential 
function of § 2 


+ This is not the same choice as was made in evaluating ‘type A’ averages; 
the effect of this inconsistency is being investigated. 
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In carrying out the averages it is necessary to evaluate the integrals 


I,= Pe G,-raie) ao] {| 3 b,- 92 dQ. 
i=1 i=! 


These can be found with the help of values of 6; tabulated by Barron and 
Domb (1955). The results are shown in table 4. 

With these results and the formulae of Appendix I we may compile the 
‘generalized Griineisen constants’ as the ‘type B’ averages of table 5. 

Comparing tables 3 and 5 we note the good agreement between averages 
derived by methods A and B. This means that to a good approximation 
we may take the same averages both at high and low temperatures. Hence 
the approximation of replacing the individual values of [ec], ete. by 
constants is not unreasonable. 


Table 4. Values of F,, and J,, (see text) Evaluated where ¢'(r?) =0 


F 


n 


- 12 


+ 128 
— 1400 
+ 16240 


Table 5. ‘Type B’ Averages Defined by eqn. (14), Calculated from 
Formulae of Appendix I 


<fePI> = —2- <[e]2I> 
{eJ™®>= 7-692 | <{e]Mfe]> 


<[e]®> = —40-20 | <[e]f4J> =202-8 
<fe}e]Fl> = 111-5 


<[e]J2B> = —21- <fe]2[e]2I> 
<{e]M4>= 61- <[e]P2 > 


§ 4. CALCULATION OF THERMAL FUNCTIONS 


The remaining functions required in eqn. (6) involve successive 
derivatives of the function 


f(a) =4etlog [1—exp (—x)], viz. /™a)= arf) 
These may be written in the form 


ax 
ad 
fe) =s2-+ cee? 


n 
got » CP eot 


ee rg imi ec he ye 2) 
where CC ean ao 1) 


N— WwW 


es (-¥ (0) 9 sy) 
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and s,/ are Stirling numbers of the second kind tabulated, for example, by 
Jordan (1950, p. 170). This result is proved in Appendix IT. 
Some values of ¢,,” are given in table 6. 


Table 6. Values of c,,° Defined by 
(—)m(ee—1yntanpiden4= Sc, where f(x) = 4x + log [1 —exp(—2)] 
11 


o= 


1 
57 1 
ELOL 120 1 
15619 | 4293 247 1 
88234 | 156190 | 88234 | 14608 | 502 | 1 


wCmoIoOnrWW 
ee 


From this table it follows, for example, that 
f(a) = —at{e® + 462% + €3*}/(e7 — 1)4. 
Graphs showing (—)"ft4/n! are drawn in fig. 1. 


Table 7. The Summed /f-functions: Values of S, (@/7') (eqn. (20)). 
Values of f!™ from fig. 1; G(A) from Leighton (1948) with 
y/«= — 0-067 (argon); (© =44-51°K for argon) 


8, 


— c 2-954 

— 4-040 2-974. 

— 1-607 3:049 

— 0-990 3-096 

—0-013 3:216 

+ 0-934 3-411 

+ 1-724 3°655 

+ 2-415 3:938 

+ 3:637 4-609 

+4-750 5:383 

+ 5-803 6-228 

+ 6:828 7:120 

+ 8:84] 8:987 
9-00 +17:75 17-765 


In eqn. (6) we need the sums ¥’,f'"(Be,) taken over the energy eigenvalues 
of the lattice. Leighton (1948) has calculated a normalized distribution 
of energy eigenvalues in terms of the dimensionless parameter A, where. 
the energy is «=}¢maxA(0<A<2) and emax is a parameter determined 
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by the force constants of the lattice. Denoting this distribution by 
G(A), the sums in eqn. (6) involve 


© : OA 
S,(—)= im} | — ) G(A) dA pager (20 
(x)= | te(S) ee) (20) 
where instead of emax we now use the parameter © = €max/2k. €max is in 

fact given by hoi 
€max = 2— (=) 
m7 \2m 
where m is the atomic mass; and G(A) has been obtained by Leighton 
Fig. 1] 
4 
$x) 


‘ lO 1°5 20 
O O-5 GE eas 
xf 5 Sa 


Graphs of the thermal functions (eqn. (19)) and frequency distribution (eqn. 
(20)). Plotted are (—1)”—1fl(a~)/(n—1)! (against «/5) and G(A) (against A). 


for two values of the force constant y/a. For argon w=Cq4, 

y = (a/4)(C4 — C42 — C44) Where, by using the work of Salter (1954), one 

may calculate c,,=3-59 x 101° dyne cm-*, ¢,= 2-04 x 101° dyne cm~ and 

C4q = 2°12 x 10!°dynecm~. Hence, with a, from eqn. (8), y/«= — 0-067 and 

@=44-51°K. The corresponding curve for G(A) is shown dotted on fig. 1. 
Values of S,,(©/7') are given in table 7. 
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§ 5. Toe Free ENERGY AND OTHER THERMODYNAMIC PROPERTIES 
or SoLtip ARGON 


With the help of the summed f-functions (table 7), the generalized 
Griineisen constants (tables 3 and 5) and the derivatives of the lattice 
energy ® (table 2) we may now tabulate the values of the coefficients A, 
(eqn. (6)) which determine the free energy. This is done in table 8. 
Finally, by differentiation in accordance with eqn. (1), we may compute 
the thermodynamic coefficients C,, «and x. These are displayed, together 
with experimental results and the results of certain other calculations on 
figs. 2-5. 


Fig. 2 
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JO 


20 


10 


eo + EXPERIMENT (FALL - Segyine- Clots 
On DAVIES & PARKE ) 
+ / -—--— ZUCKER 
ead ete AAD Bomueee coc co, KANe 
i ce oe DEGYE ( cafe. ZucKeER) 


O 20 4o 60 


80 


Ue K. 


The specific heat of solid argon at constant pressure (C,,) as a function of 
temperature (7'). 


+ experimental points obtained from a smooth curve drawn through 
the results of Hill et al. (see Dobbs and Jones 1957). . 
the generalized Griineisen theory of this paper (calculated from 
eqns. (1) and (5) with the help of table 8). 

— — —anharmonic calculations of Zucker (1958). 

—.—.—.—harmonic theory of Kane (1939). 

ashe calculation based on Debye’s model (Zucker 1958). 
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Table 8. The Expansion of the Free Energy: Values of A,./8 (see eqn. (6)) 
found from Tables 2, 3, 5 and 7. (The Potential Energy Contri- 
bution to A, is BD”, see Table 2.) 


et | T(°K) A,/B Aina Aue Ale 
in keal mole! 

oe) 0) — 1-852 + 0-0067 14:12 —127°5 1059 
9-0 4-94 — 1-852 +0-0063 14-12 —127°5 1059 
4:5 9-89 — 1-853 0 14-08 — 127-8 1058 
3:5 12-74 — 1-854 —0-0103 14-03 — 127-9 1059 
3:0 14-84 — 1-856 —0-0218 13-98 — 128-0 1059 
2-5 17-80 — 1.858 —0:0430 13-92 — 128-0 1060 
2-0 22-25 —1:866 ° —0-0846 13-84 —127-9 1061 
1-6 20-67 —1-884 —0-1728 13-74. —127-7 1060 
1-25 35:60 — 1-905 —0:2547 13-69 — 127-5 1058 
1-00 44-51 — 1-944 —0-3888 13-67 — 127-4. 1055 
0:75 59-34 — 2-028 — 0:6296 13:72 — 127-6 1051 
0-55 80-92 —2-186 —0-9984 13-88 — 128-4 1048 


a 
+ DOBBS ET AL .1956 
—— DAVIES & PARKE 
---— ZUCKER 
DEBYE (cale.ZUCKER) 
O 20 40 60 80 


1/6 

The variation with temperature (Z') of the expansivity («) of solid argon. 
+ experimental results of Dobbs eé al. (1956). 

evaluated from the theory (eqns. (1) and (5) and table 8). 


— ——-—anharmonic calculation (Zucker 1958). 
See harmonic (Debye) theory (Zucker 1958). 
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The variation of isothermal compressibility («7) with temperature (7’) at 
zero pressure. 
+ experimental results of Barker and Dobbs (1955). 
M experimental result of Stewart (1956) at 65°K. 
the harmonic theory of this paper. 


— — —anharmonic calculation (Zucker 1958). 
Debye-type calculation (Zucker 1958). 


Fig. 5 
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The isothermal bulk modulus (1/«7) as a function of pressure (/). 
+ experimental results of Stewart (1956) at 65°K. 


———the generalized Griineisen theory of this paper, 
anharmonic calculation (Zucker 1958). 


5000 /0000 /5000 
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§ 6. A Direct EsTimMatTE OF THE SECOND GRUNEISEN CONSTANT 
FROM EXPERIMENTAL Data 


In the work of §§ 3-5 we have obtained generalized Griineisen constants 
by making reasonable theoretical estimates based on a single inter- 
molecular potential function. It is interesting that, just as one can 
estimate the first Griineisen constant, []" = [(V/e)(d«/aV)],,= —y directly 
from measurements of the expansivity, specific heat and compressibility, 
so can one estimate the second Griineisen constant, 


fe) [(PAfe)(oPe(oV*) |, Sy" 


(say), directly from measurements of the variation of compressibility 
with temperature, together with certain further information. In the case 
of argon we have obtained an estimate of y’ by this method which agrees 
satisfactorily with the theoretical estimates given earlier (tables 3 and 5). 


Table 9. Experimental Values of the Quantities in eqn. (24) for Argon, 
using y=2-44 and taking © from the Intermolecular Potential 
of eqn. (7) 


» 2 


7 
y aV2 


U 


(all in keal mole“) 


0 12-4 —1-85 
0-46 12-1 —1-83 
1-44 2027 —1-75 
2-41 8-5 — 1-62 
3°60 6-5 —1-49 


Taking the potential terms and thermal terms into the right side of 
eqn. (3) and making the generalized Griineisen assumption of treating 
fe]™ and [e]® as constants, we may eliminate the thermal functions 
f™ and f®! to obtain three independent relations. 


BpV = —BV(0®/0V) —[e]*B(U — 9), (21) 
(V[k)(a/«) = — [e]2C,/k, “Oy ee a (22) 
B(V /x) =BV7(e?O/0V?) + [e]PIB(U — ®) — [e]*PC,/k. (23) 


Equations (21) and (22) are familiar results. Equation (23) is less 
familiar and is the one that we wish to use. It may be rewritten for 
convenience as 

Vice “(o-Dia V2) °C Tay)... s  ((24) 
Using the sources cited earlier we have obtained V, «, C, and U from 
experimental results, we have assumed the first Griineisen constant to be 
2-44 (which is the mean of the four values of y given by Dobbs and Jones 
(1957), Appendix III, for 7’=20°K, 40°K, 60°K and 80°K) and we have 
calculated ® and V2(d?@/aV2) from the intermolecular potential used 
earlier. The results are given in table 9. 


PM. Z 
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A graph combining these results is shown in fig. 6. From it we estimate 
y’=[e]2=11-2+1-0. This value may be compared with the theoretical 
estimates given earlier: ‘type A’ averages gave [e]"!= 10-0 and ‘type B 
averages gave [e]#=9-1. The value obtained for y’ by the method we 
have just used depends sharply on the value adopted for y. Thus if we 
had used the value y= 2-8 (which is close to Barron’s—and our—‘ type B’ 
average) we would find y’ = 14-5 + 1-0. 


Fig. 6 


(L+vte,-vF8) = y'(u-4) 


OWS O02 


Os O-4 
(U- OY Rect mote” ' 


V aD 


Test of the relation — = V? av? —yTC,+y'(U—©) for y=2-44. @® is the 
KT my 
: ; Pn eee y V de 
lattice energy (§2). y is the‘ first ’ Griineisen constant { = ( —— =) 
e OV Avy 
3 : Val 0 
and y’ is the second Griineisen constant ( =(— a ) . (See §§1 and 
e OV! ay 


6.) The remaining quantities are obtained from experiments (Dobbs 
and Jones 1957). 
The slope of the line yields y’=11-2+1-0. 


§ 7. CONCLUSIONS 


The main results of the calculations, as displayed on figs. 2-5, show that 
our generalized Griineisen theory is in reasonable agreement with experi- 
mental results on argon from absolute zero up to about 40°K. Above this 
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temperature it is definitely inferior to the theory of Zucker. This is 
perhaps not surprising because our theory gives only a partial account of 
anharmonicity ; in effect we have assumed harmonic modes of vibration, 
with frequencies depending parametrically on the volume. Zucker on 
the other hand deals with an anharmonic, one-atom model so that it is not: 
surprising that his theory is less successful at low temperatures. 

In favour of our theory is that, now that the main calculations of ‘type 
B’ averages and thermal f-functions have been made, it can be extended to- 
deal with other rare gases. Certain other extensions (e.g. to anisotropic 
crystals) are also possible. 

The y values obtained by us are: 3-11 (‘type A’) and 2-76 (‘type B’). 
It is not surprising that the latter is close to Barron’s value of 2:82 because 
our method of averaging is the same as his. Zucker obtains y= 2:8, 
while Dobbs and Jones quote experimental values in the range 2-4-2:8. 

Our three estimates of the second Griineisen constant, 


= [(V#/e)(o%«/0V*) 14,5 


are: 10-0 (‘type A’), 9-1 (‘type B’) and 11-2 (from experiment, see § 6). 
It seems reasonable to conclude that y’ is a new approximate ‘constant’ 
(like y) with a value near to 10 for argon. This ‘constant’ has not 
previously been discussed and estimated. It clearly constitutes a second 
measure of anharmonicity (y is the first) and it would be interesting to 
investigate its value in other systems. 


APPENDIX I 


To carry out the calculations of ‘type B’ averages in § 3 and other 
manipulations of the type discussed here, there are some elementary 
formulae which are of great value. 

Let f(x) be any differentiable function of x. Denote the type of 
differentiation by superscripts as follows: 


f= — fils Gear g 7 ss aa 1 oe ACB 


and denote the type of function to be differentiated by brackets: 


1 anf. 


(f= fM= = ; [f]m= Fee {fy = an a 


(26) 


Thus the symbolism admits nine types of nth derivatives. The con- 
version from ( )to[ ]in either of the rows (25) and (26) is immediate : 


Hie een ea (ims 1.2 ee ts 2 (27) 


Z2 
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The conversion from [ ]to{ }is more complicated and the relevant 
formulae are : 


fu sf 1S, 
f2 = fi —f%, | 
ea oS) 
f® = fi cE BY sb afi, 
f™ = fis = 6f 3 at. lif? = 6f 4, 
ee ef 10) 


f%=a=fr+fy, 
f8=fS8+ 3f2+ fa, 
fMafsofF+7fas+sfa, 
The results of (28) and (29) are true for all forms of the function being 
differentiated. 
(H=1FP, 


£3 = (f1— LF™, 
£f} = [f1 — SLfPPLAI + 2 F1%, aes 
| \ 


(29)T 


£f}O = Lf — ALF Cf — BLf]2? + LAL FIM]? — 6 LF] 


ie aalye, 

EB 

ele seat paid east ty a 

[PI =f} + AL PPOLP}® + BGPO2 + GLPIO2LP@ + [Pr4, 


The results of (30) and (31) are true for all types of differentiation. 

With the help of eqns. (26)—(31) we can easily make conversions between 
any of the nine types of derivative. Thus one can form the expressions 
required for [e]®!, [e]4#?, etc. in eqn. (6) in terms of the derivatives of 
¢ in eqn. (15). 


31)I 


vir WeN yt exer 
To prove the formulae of (19) for the nth derivative of 
f (x) = 4a + log [1 —exp(—2)], 
set e2=u, € “=v so that 


d d d 


de dinu do’ 


{ The coefficients in (28) and (29) are Stirling numbers of the first and 
second kind respectively; see Appendix IT. 

{ Equations (30) and (31) are familiar to statisticians as giving, respectively 
the conversion from cumulants to moments and vice versa (see Kendall (1947), 
p. 62). 
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Now differentiating » times the formula 


1 
f%x)=h+ oa tie —$+ (32) 


we have 


_ d 1 d 1 
(n+) = S| Pie (yee ee § 
oad d(ln uw)” (= i) oe d(In v)" (; +) 2 Pak 
Invoking one of the main properties of Stirling numbers of the second kind 
s,) (see Jordan, 1950), viz. 


dq” d \j 
dna =2 s,Jw i ma) ao ot ate (os) 
we therefore have 
» s,Jw(—)j! (w—1)"-4 
(2+-1) Se ee ee ee 
if iej= ce = ‘ oe yD 


S (= )rs,i0ij! (1-0 
fee era ene) 


— 


or, finally, writing v=1/w in (36), 


3S (-Yentf e—1yu 


fe) = a aa . (37) 


Equation (37) is virtually the same as (19) of the text and comparing it 
with (35) we see immediately the symmetry relation c,°=c,”-°*! among 
the coefficients; because of this relation it is necessary to work out the 
c, only for w=1, 2,... [~+1/2]—as is clear from table 6. 
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ABSTRACT 

Measurements in the form of oscillograms have been made of the electric 
field HZ, the electric displacement D, and the strain 2, on both single crystal 
and ceramic barium titanate over a wide temperature range. Below the 
Curie point results depend strongly on the type of sample, and are given an 
interpretation in terms of domain change and piezoelectric mechanisms. 
Above the Curie temperature all specimens show electrostrictive behaviour 
with a coefficient Q3,; which within experimental error is independent of 
specimen type, temperature, and applied field. 


§ 1. INTRODUCTION 

Iv is well known that barium titanate exhibits strong electromechanical 
effects which may be classed with some difficulty as piezoelectric or 
electrostrictive. A general review has been given by Kay (1955), but 
there have been rather few experimental results published in which the 
strain has been measured directly, the work of Caspari and Merz (1950) 
and Schmidt (1956) being notable here. The investigations to be 
described in this paper have some affinity to those of Schmidt; the strain 
gauge used has already been described (Allsopp and Gibbs 1957). 

Briefly, a single crystal or ceramic specimen, with coated electrodes, is 
mounted in a temperature controlled enclosure so that signals proportional 
to the field #, the electrical displacement D, and the strain x can be 
obtained and displayed on a cathode-ray tube, there being facilities for 
the absolute measurement of each variable. In general the calibration 
factors of the apparatus are known to a precision of better than +1%, 
but in many cases the overall accuracy is worse than this owing to greater 
uncertainty in the dimensions of the specimen and its electrodes, or in the 
case of very small strains, owing to a poor signal to noise ratio. 

Any pair of variables can be selected to give X and Y deflections, thus 
giving three types of figure: D/H, 2/E and x/D, and these can be 
photographed. 

For most observations the voltage applied to the specimen had a 
frequency of 5c/s, this low frequency being desirable mainly to avoid 
mechanical resonances in the strain measuring system. 

The electrodes covered the whole of the large faces, and in preliminary 
experiments the crystal was supported under light spring tension between 
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two metal balls of radius 1mm which were mechanically connected to 
the strain gauge, so that a local, and not average, value of the strain 
was displayed. In contrast to this, the electric field was presumably 
uniform (apart from the effects of twinning), and the displacement signal 
proportional to integrated current, thus being a volume average over the 
crystal. The effect of the thrust of the supports on the crystal’s behaviour 
was thought to have been obviated in later experiments, and some 
averaging of the strain obtained, by the use of soft leather pads as 
supports. The evidence for this is that the figures then remained nearly 
the same over a wide range of spring tensions, which was not the case 
with the metal ball supports. 


§ 2. RESULTS WITH SINGLE CRYSTALS 
2.1. Ideal Model 


In order to interpret the results, it is convenient first to recall what 
might be expected from the simplest ideal crystal. Here we consider a 
‘e-domain ’ crystal in the form of a thin plate with the polarization vector 
everywhere normal to the large surfaces, and assume that the polarization 
suddenly reverses throughout the crystal when the applied field passes 
through the coercive value H,. The D/E figure is the familiar idealized 
hysteresis loop of fig. 1(a). Assuming that the crystal behaves in a 
simple piezoelectric manner and that the piezoelectric coefficient ds, 
reverses its apparent sign with domain reversal, it is clear that the v/E 
and a/D figures shown in figs. 1 (6) and 1 (c) would result. 


Fig. 1 


(0) 


Ideal D/E, x/E and x/D figures for c-domain crystal. Paths during switching 
are shown dotted. 


2.2. Results at Room Temperature 

Crystals which, as grown, approached the ideal ‘all c-domain’ 
configuration were selected by examination under polarized light. These 
will be referred to as ‘ group IT’ crystals. A typical set of figures obtained 
with these is shown in fig. 2 (Pl. 47), in which the peak applied field 
increases progressively from (a) to (e). 

The major qualitative deviations from the ideal of these group II 
crystals may be listed as follows: 

(i) At low fields, where the peak value of EH is less than the value of 
HE, obtained from a large field loop, there is often marked asymmetry in 
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the figures, though the D/E figures do not resemble the straight lines 
expected of a single c-domain or a mixture of antiparallel c-domains. 

(ii) Although the values of saturation polarization may be high 
(typically 20-24 ,.c/em*), hysteresis persists out to the tip of the D/H 
figures. 

(iii) The cross over at the base of the 2/D figures is not as clear-cut as 
in the ideal. 

(iv) Even at the highest fields which it is possible to apply without 
breakdown, the tips of the a/# figures are still curved and have not turned 
into the straight lines expected of a simple single c-domain. 

When crystals which optically show a greater degree of twinning 
(group I) are examined, these deviations become on the whole more 
marked, as typified in fig. 3 (Pl. 48). In addition: 

(v) The strain obtained when the peak applied field is less than 2, is 
extremely small. 

(vi) The a/D figures are more rounded at the base and less steep at the 
sides, though the sides remain approximately straight. 

It thus seems that even the best available group II crystals, at least 
without special treatment, show marked qualitative deviation from the 
simple ideal of §2.1. There are also quantitative deviations, and since 
these throw light on the performance of the crystals, we now consider 
three possible simplified mechanisms for which a rough value of the slopes 
of the x/D figures can be estimated. 


2.3. 
2.3.1. Single domain piezoelectric model 
This is what has been assumed in the saturated regions of fig. 1. We 


have Ox a 
3p °38 
by definition, and hence Ox D 
ie Aly. 33/€3> 


where e, is the c-axis dielectric constant. These equations should thus 
give the slopes of the appropriate parts of the x/e and «/D figures if 
complete saturation is reached. A value of 2-5 x 10~%e.s.u. for dj; may 
be deduced from the measurement of d,, by Caspari and Merz (1950) if, 
following these authors, we take ds3/d3;= —2-°3. Substituting «,=200 
(Merz 1953) we then have 


Ox = 1-25 x 10-8 e.s.u. 


0D 


This process will be referred to as the p-mechanism. 


2.3.2. 180° rotation model 

This corresponds to the switching regions of fig. 1. Since d33 reverses 
sign under a field H,, the resultant change in strain is 2d,,H,, and the 
corresponding polarization change is 2P, or twice the saturation 
polarization. Since D=4nxP+£, the slope of the x/D figure which 
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corresponds to this mechanism is given by 


Au _ dg3He 
(55) a ~ 4rPs° 

The same result is obtained if a fraction m only of the crystal is assumed 
to be built from reversible c-domains, the rest being unchanged, provided 
that a ‘natural’ strain x of the (unclamped) c-domain component is 
assumed to result in an average strain ma of the whole crystal. 

If we substitute H,=3e.s.u. (typical of group II), d3; as before and 
P,=7-5 x 10%e.s.u. (Merz 1953), we obtain 


(<3) = 8x 10-" e.5.u. 
AD] 180 


This process will be called the 180-mechanism. 


2.3.3. 90° rotation model 


When a domain initially polarized at right angles to the field direction 
(a-domain) turns into the field direction and thus becomes a c-domain, 
an elongation of approximately 1° results from the tetragonality of the 
unit cell, and an increment of polarization equal to the domain volume 
times Ps occurs. If the crude assumption is made that it is legitimate 
to make a simple volume average of these quantities over the whole 
crystal, it is clear that the slope of the x/D figure is given by 


Az\ _ Strain due to rotation 
(=) ne 4nPs 

The increment of D will in general be greater than that stated owing 
to the increased dielectric constant of the a-domain, so that the slope is 
liable to be rather lower than the above figure, especially if the switching 
occurs at high fields. 

This process will be called the 90-mechanism. 


=1-1x 10-8 e.s.u. 


2.4. Comparison of these Mechanisms with Experiment at Room Temperature 


It will be observed that the slopes of the w/D figures to be expected from 
the p- and 90-mechanisms are nearly the same, and about one thousand 
times greater than that to be expected from the 180-mechanism. Thus 
the p- and 90-mechanisms would be difficult to distinguish by examination 
of the x/D figures alone, but the rate of change of « with Z is likely to be 
very different for the two mechanisms, being a simple straight line of slope 
ds, for the p-mechanism, and probably non-linear and hysteretic for 
the 90-mechanism. When the 90-mechanism is dominant in producing 
strain, we may therefore expect a large amount of hysteresis in the a/H 
figures, but little in the x/D figures. The 180-mechanism is likely 
to give hysteresis in both figures, and the p-mechanism in neither. 

Values of (Az/AD),,5 calculated from the switching part of the x/D 
figures range from about 5 x 10-!! to 5 x 10-e.s.u. for group II crystals, 
i.e. from 5 to 50 times the value calculated in § 2.3.2. It therefore appears 
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that even in the most ideal untreated crystal available, some 90-mechanism 
occurs, though the amount of a-domain necessary may be only about 
0:5%,. 

Values of the apparent d3; calculated from the limiting slopes of the 
x/E figures of group IT crystals also always gave results greater by a factor 
of between 5 and 10 than the value of Caspari and Merz. We conclude 
that the 90-mechanism persists even to the highest fields. 

The average of dx/0D over eight crystals in the high-D region was 
1-0 x 10-8e.s.u., in fair agreement with either p- or 90-mechanisms, and 
of these crystals only one gave a value exceeding the calculated (Ax/AD) 
though this high value was over twice the theoretical. 

Some further evidence on these interpretations is given by the results 
on two special crystals shown in fig. 4 (Pl. 49). Here (a), (b) and (c) refer 
to a group II crystal, (a) shows the effect of initial application of a field, 
(>) shows the effect of maintaining the field for two hours (little further 
change is brought about by field alone) and (c) shows the result after 
having raised the temperature to near the Curie point and again having 
cooled to room temperature, the field being maintained throughout. This 
set of figures seems to indicate a close approximation to the ideal with 
only the p- and 180-mechanisms in operation, nevertheless a quantitative 
examination still indicates discrepancies. The saturation polarization 
measured from fig. 4(c) is 24+2yc/cm?, in satisfactory agreement with 
the accepted value, but the e«, from the saturated part of the D/E figure 
exceeds 2000, compared with the low figure of 200 obtained by Merz. The 
value of d3, calculated from the tip of the v/# figure is 5-2 x 10-*e.s.u., 
still over twice that obtained by Caspari and Merz, while that obtained . 
from the switching region (if assumed to be 180-mechanism) is about 
25x 10-e.s.u. If the Caspari and Merz d,, value is accepted, we must 
conclude that the 90-mechanism is still not completely suppressed, though 
the residue of a-domain is only about one part in 1000. Probably one of 
the factors contributing to its greater suppression in the work of Caspari 
and Merz is the much slower rate of change of field employed by them. 

Results on a crystal showing characteristics of the opposite extreme 
are shown in fig. 4 (d), (e), (f), (g) (Pl. 49). This group IT crystal is probably 
highly atypical but is of interest because of the large degree of 90- 
mechanism displayed. It had been subjected to two heating cycles when 
these results were obtained, after subsequent cycles it reverted to a 
condition very similar to that which produced results in fig. 4(@), (bd), (c). 
The saturation polarization in fig. 4(g) was 23+2c/cm?, yet the 
mechanism of reversal is evidently very different from that in fig. 4 (c). 
The large pendant loops in the «/H figure of fig. 4(g) are believed to 
correspond to the 90-mechanism, and may be compared with fig. 3 (a) of 
Kay (1955). About 7% of the crystal must be converted to a-domain 
to account for the observed strain. 

In the a/D figure, the steeply sloping sides will be seen to be made up of 
two parts, a heavily exposed top part corresponding to near saturation, 


90° 
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and a lightly exposed lower part corresponding to switching. The fact 
that these are nearly collinear provides some confirmation of the near 
equality of dx/AD and (Ax/AD),), though the confirmation is less striking 
than at first appears both on account of the steepness of the slope and 
because the upper region probably involves a mixture of p- and 
90-mechanisms. 

To account for the small amount of hysteresis in the 2/D figure of 
fig. 4(g), it seems necessary to assume that the 90-mechanism and 
180-mechanism always operate at essentially different parts of the cycle. 
An idealized scheme which gives results resembling those obtained is 
shown in fig. 5. Here it is assumed that the crystal is wholly +¢-domain 


Fig. 5 


Ideal figures for crystal of fig. 4 (d), (e), (f), (g) (Pl. 49). Switching paths 
shown dotted. 


from ‘a’ to ‘d’, and converts to 79% a-domain at ‘e’. From ‘e’ to 
‘f’ the slope is greater on the D/# figure but less on the w/D figure owing 
to the higher dielectric constant of the a-domain component. From ‘f’ 
to ‘g’ the c-domain component is assumed to reverse, while from ‘h’ to 
‘i’ the a-domain component also converts to reversed c-domain. In 
practice the various transitions are not so sharply separated. 

As with other group II crystals, the strain is large even when fields less 
in peak value than H, are supplied, as in fig. 4(d). The crystal is here 
highly asymmetrical presumably being polarized mainly in the negative 
direction, the a-domains being formed when the field is in the positive 
direction. 

In contrast, little strain is observed in group I crystals with fields 
lower than #,, but apart from this feature, it appears that group I crystals 
can be understood as having the 90-mechanism dominant but spread over: 
a wider range of field values. 


2.5. Experiments on Single Crystals at other Temperatures 
2.5.1. Hapervments below room temperature 


It has been shown by Merz (1949) that the apparent saturation 
polarization of a crystal suddenly decreases at each of the lower transition. 
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temperatures. The one group II crystal investigated here showed similar 
effects, but in contrast to Merz’s results the saturation polarization also 
diminished continuously with decreasing temperature as shown in fig. 6. 
In the case of a group I crystal the decrease was continuous with no 
evident discontinuities, as also shown in fig. 6. It thus appears that the 
continuous decrease can be attributed at least in part to the gradual 
inhibition of domain switching mechanisms with decreasing temperature, 
and it is perhaps surprising that the continuous decrease was not found by 
Merz, whose crystal showed a room temperature Ps of 16 c/em?, and 
must therefore have been twinned. 


Fig. 6 


4 A group I Xtal 
o A group IT Xtal 


Psat (¢/ em?) 


Temperature(°C) 


50 100 


-150 ~100 


Saturation polarization of crystals as a function of temperature. 


Observations of the x/H and x/D figures have also been made and these 
did not change greatly in shape down to a temperature of — 40°c, though 
the attainable polarization and strain varied widely. In particular, there 
was little change in the slope of the steep part of the x/D loop, which 
seems odd in view of the change through 45° of the polarization vector at 
the 0°c transition, and the consequent necessary alteration of domain 
change mechanisms. 

Below —40°c, x diminishes more rapidly than the apparent Ps, and 
below about —60°c the x/D figures no longer have their characteristic 
U shape. From —60°c to —150°c the total strain at ‘saturation ’ 
polarization is about 10-4, but the figures undergo complex changes which 
will not be pursued here as it may not be wise to generalize from the 
behaviour of one crystal. 


2.5.2. Huperiments above room temperature 


The temperature range of greatest interest is that close to and on both 
sides of the Curie point, but observations in this region are sometimes made 
difficult either by excessive conductivity or breaking of the crystal. 
Results on a crystal of each group are shown in fig. 7 (PI. 50). As would be 
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expected, the v/H and x/D figures both become parabolic in the paraelectric 
region well above the Curie temperature, and a little above the Curie 
temperature ‘double’ D/E figures appear, where the crystal becomes 
ferroelectric only at high fields, as observed by Merz (1953). This field 
induced transition is more sudden in the group II crystals and is accom- 
panied by a switching strain in the a/H figure bridging the parabolic 
(paraelectric) region and the more linear (ferroelectric) region. 

The electrostrictive constant Q,, defined by xv33=Q33P?32 was measured 
from the a/D figures above the Curie temperature and with one exception 
(a very thick crystal) found to be the same within experimental error for 
all crystals, and to be invariant with field and temperature in the range 
130°-150°c. The weighted mean value was calculated as 

Qeq= (6:4 + 0:3) x 10-4 e,8.u. 

If it is assumed that Q;, varies slowly and continuously with temperature 
as the Curie point is traversed, and that it is not necessary to include 
higher order terms in the 2/P relation even when ferroelectric polarization 
sets in, it is possible to relate the value of Q,, to the saturation polarization 
and lattice parameters (Merz 1949) and also to relate it to the piezoelectric 
constant d;, (Caspari and Merz 1950). Taking the lattice measurements 
of Kay and Vousden (1949) and a Ps of 18 wc/em?, we obtain in this way 
a Qs, of 1-1 x 10-“e.s.u. just below the Curie point, in poor agreement 
with our directly determined value. We believe that the underlying 
assumptions rather than any of the measurements are at fault, most 
probably it is necessary to include higher order terms in the x/P relation. 
If this is so, calculations of d3, also become subject to considerable 
revision, quite apart from the higher value of P; which has now become 
accepted. 

In the region just above the Curie temperature where transitions between 
para- and ferroelectric phases can be field-induced, it is possible to make 
fairly good estimates of the switching strain Ax, from the x/H and 2/D 
figures of group II crystals, where the transition is sharp. Values of 
Az, of 2-2, 1-7, 1-3 and 1-:3x 10-3 were obtained in four experiments, 
which values may be compared with 3-3 x 10-* calculated from the x-ray 
data of Kay and Vousden (1949). The apparent saturation polarization 
in this region has also been estimated by linear extrapolation of the 
extremes of the D/H figures to zero field. Values of Ps were also 
somewhat variable and lower than expected, being respectively 15, 11, 10 
and 6 .c/em*. ‘The discrepancies are thought to arise from the switching 
of only part of the crystal, the frequent shattering of crystals under such 
severe stress being understandable. Accurate temperature control is also 
difficult when such relatively large amounts of electrical energy are being 
dissipated in the crystal. 


§ 3. RESULTS ON CERAMICS 


A wide range of ceramic samples has been investigated, mostly in the 
form of discs about 1 cm in diameter and 1mm thick. Some were as pure 
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as possible, and others contained intentional additives. It is difficult to 
obtain accurate and reproducible results in many cases, partly because the 
behaviour of the specimens depends so much on previous history, partly 
because different specimens made as nearly as possible under the same 
conditions may differ considerably, and partly because of electrode 
difficulties. Preliminary experiments showed that many specimens do 
not undergo uniform strain when a field is applied, but may warp or ‘dish’ 
(as has been noted by Schmidt (1956)), hence the ball ended support was 
always used. The strain often varied by as much as 25% at different 
parts of the same specimen. 

Owing to the high dielectric constants involved, minute gaps between. 
the electrodes and the specimens have a very large effect. Electrodes of 
fired-on silver paste are generally satisfactory on dense ceramics, but in 
some cases the silver penetrates the body of the material appreciably, 
making the field non-uniform and uncertain. This difficulty was. 
prominent with the more pure materials. Evaporated films of metal are 
initially good, but easily damaged. Further, their resistance is liable to. 
become very high after heating; presumably they are broken up by 
relative movements of the grains when a phase change occurs. A 
technique of sputtering. platinum, although producing good electrodes, 
was found to overheat and reduce the specimens. 

A silver paste for cold application (‘ Elargol W4’) was found generally 
satisfactory, but tended to separate from the ceramic below —60°c when 
strong fields were applied. A mixture of silicone grease and. precipitated 
silver was found useful in this low temperature range, though inconvenient 
and of higher resistivity than the other electrodes. 


3.1. Results at Room Temperature 


Results on a number of different samples are shown in fig. 8 (PI.51). It 
will be seen that the x/D figures are roughly parabolic in form, but may differ 
from the ideal in that they may possess asymmetry (prepolarization), 
flattening of the vertex, and hysteresis. Values of Qs,’ in the relation 


Xz =Q33 P53" 


were calculated from the parabola of best fit to the x/D figures, and 
ranged from 3-0 to 7-9 x 10-8e.s.u., in rough agreement with the results 
of Schmidt (1956). Unlike the latter author, however, we find considerable 
variation with specimen material and considerable departures from the 
parabolic form, perhaps because larger fields (about 20kv/cm) were often 
used by us. The flattening of the vertices has been illustrated in a 
previous paper (Allsopp 1957) and it is roughly true, as might be expected, 
that Qj3’ is smaller in those specimens which show greater flattening. 

The flattening is thought to be similar in mechanism to that suggested 
for group I crystals in 2-4, and this dominance of 180° changes at the 
vertex of the curve would also account for the hysteresis observed, the 
direction of traverse being invariably in accord with this principle. 
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Although giving a similar hysteresis, conductivity of the specimens does 
not seem to be the cause, both because of the high insulation 
resistance found and because of the insensitivity of the 2/D figures to 
change over at least a 10: 1 range of the frequency of the applied field. 

The coefficient Q,' has also been measured for two samples of pure 
barium titanate ceramic, by suitably modifying the mounting so that 
strain was measured in a direction at right angles to the applied field. The 
values obtained were —0-45x 10-3 and —0-51 x 10-¥e.s.u., while the 
same specimens had Qs; values of 3-6 and 3-8 x 10~1%e.s.u. respectively. 
The ratios Qs,’ to Qs’ were thus — 0-12 and — 0-13, which may be compared 
with — 0-38 and — 0-23 calculated by Mason (1950) and Devonshire (1954). 
The low value of this ratio obtained by us may be the result of the rather 
high porosity of the samples. 


3.2. Prepolarization 


Results on a prepolarized ceramic containing 1-3°% of cobalt oxide are 
shown in fig. 9 (Pl. 52), where the development of the quasi piezoelectric 
behaviour at small fields to the nearly parabolic behaviour at high fields 
may be traced. It must be remembered that the zero of neither D nor x 
can be established on these figures. 

In such materials a large proportion of the domains are thought to be 
locked in the direction nearest to that of the prepolarizing field (we shall 
call these c-domains, and all others a-domains). Kay (1955) has suggested 
that the strain behaviour of a prepolarized ceramic might be accounted for 
by the linear piezoelectric effect on these locked c-domains. It appears 
to us that such a mechanism, although it must contribute, cannot account 
for the magnitude of the strains observed. In fig. 9(b), for example, the 
applied field is 70e.s.u. peak to peak, and the total strain 2-1 x 10-4, 
equivalent to a d;, value in a piezoelectric crystal of 3 x 10-%e.s.u. This 
is a little greater even than the dj, value of a single c-domain crystal, which 
value could not be closely approached even by the most completely 
polarized ceramic. Against this argument it must be admitted that the 
field on the c-domains will be greater than that applied, owing to their 
relatively low dielectric constant, but it is difficult to see how the field and 
strain averaging could be done to produce such a big effect. 

Although the problem is immensely complicated by the large internal 
stresses and fields which must exist when domain changes occur in the 
randomly oriented grains, and although the magnitude of the shear 
piezoelectric constant d,; is unknown, we suggest, in agreement with 
Mason (1950) that much of the strain in these prepolarized ceramics arises 
from 90° domain changes. 

The prepolarized material contains then a preponderance of say 
+ c-domains, with a proportion of a-domains and few —c-domains. Small 
fields in the apparently linear a/H region effect conversion only between 
+c- and a-domains, there being very little conversion to —c-domains. 
Thus under small fields the behaviour depends on the tendency of a 
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reverse field to form or enlarge a-domains, while a forward field tends to 
remove them. The larger slope and greater hysteresis in the reverse field 
region of the D/E figure at moderate applied fields (such as in fig. 9 (b)), 
together with the reduced slope of the x/H# and x/D figures in the same 
region, is thought to indicate the onset of 180° changes. An equivalent 
process would be the conversion of a-domains to —c-domains at a rate 
approaching that of conversion of +¢- to a-domains. 

The low-field asymmetry in the switching process, characteristic of 
prepolarization, is presumably produced by the internal field of the 
excess of c-domains, and it is perhaps surprising that asymmetry in the 
figures persists even when a large degree of switching to —c-domain must 
occur, as in fig. 9(e). Possibly the preference for the +c direction is the 
result of partially compensating volume charges, formed by leakage or 
breakdown, on the polar ends of the original ‘ locked’ +c-domains. This 
may be one of the factors explaining the greater effectiveness of prepolariz- 
ing at a high temperature, where the conductivity is higher. 

If the domain change mechanism is in fact the dominant one in 
producing strain, it would be expected that very strong fields would 
reduce the slope of the x/H figure owing to almost complete elimination of 
a-domains, as in the case of single crystals. There seems a slight tendency 
for this to happen, but it is not nearly so marked as with crystals. This 
is not surprising in view of the small distinction between c- and a-domains 
in many of the randomly oriented grains of the ceramic, and the consequent 
very intense field required to produce alignment. 


3.3. Temperature Effects with Ceramics 


Typical x/D figures for two different ceramics at various temperatures 
from —150°c to +135°c are shown in fig. 10 (Pl. 53). The following 
features are common to all observations made : 

(i) The general form of the x/D figures is unchanged in the temperature 
range between — 150°c and the Curie point. In particular, no changes are 
observed on passing through the two lower transition temperatures. 

(ii) Above the Curie point the x/D figures are parabolic to within the 
limits of experimental error. The Q;,’ coefficient in this range is also 
independent of temperature and specimen type (i.e. to small amounts of 
additives, small variations in BaO: TiO, ratio, and firing temperature and 
time). The mean value of Q3, for five specimens each at several different 
temperatures is (5-8 + 0-3) x 10-e.s.u., and is thus the same (just within 
the estimated error) as Q3, for crystals. 

(iii) Hysteresis in the z/D figures tends to increase below —100°c and 
above +100°c. The latter effect is usually attributable to conductivity, 
and varies widely from specimen to specimen. 

(iv) In high purity material, there is a marked flattening of the vertex, 
especially just below the Curie point, where the decay-effect has a short 
time constant (cf. Allsopp 1957). The dip in Q33' on traversing the Curie 
point, found by Schmidt (1956) may be of the same origin. 


P.M. 2A 
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(v) It appears that the eatent of the flattening increases with decreasing 
temperature, although the situation is confused by time effects. 

Results on the variation of Qs,’ with temperature tended to be erratic 
at very low temperatures, probably because of electrode difficulties, and 
the slowness of time effects. What appeared to be the most reliable 
results indicate a smooth fall of Q33’ with temperature, the coefficient 
changing by a factor of about two over the range + 120°c to — 150°c. 


3.4. Field-induced Transitions in Ceramics near the Curie Temperature 


The figures for a pure ceramic material as the temperature is raised 
through the Curie point are shown in fig. 11 (PI. 54). It will be seen that 
‘double ’ hysteresis loops appear, the figures being comparable with those 
obtained with group I crystals. The paraelectric to ferroelectric transition 
here extends over a larger field range (i.e. a larger fraction of the cycle) than 
with single crystals. Less pure ceramic material failed to show these field- 
induced transitions, which do not appear to have been previously recorded 
in ceramics. 

The strains in a direction at right angles to the applied field for the same 
material are shown in fig. 12 (Pl. 55). As would be expected, the x/# figures 
show transitions, just above the Curie point, between the parabola 
characteristic of the paraelectric phase and the pair of straight lines 
(approximately) of the ferroelectric phase. However, the remarkable 
feature is the sign of the strain, there being an extension in the transverse 
direction, as well as in the in-line direction, on conversion to the 
ferroelectric phase. The result may again be dependent on the comparative 
porosity of the pure material, and it has not been so far possible to check 
this, as the only samples of greater density were also of less purity, and 
did not show the effect. 
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ABSTRACT 


The effect of phonon drag is calculated, taking account of electron—phonon 
Umklapp processes, and anisotropy of the Fermi surface. For Na, the 
observations are perfectly fitted with a spherical surface. For K, Rb and 
Cs, the surfaces must bulge more than half way towards the faces of the 
Brillouin zone, whilst in Li it probably touches the zone boundary. 


ACCORDING to the conventional free-electron theory of a monovalent 
metal, the absolute thermopower, Q, at low temperatures should show the 
following features : 


(a) It should be small in magnitude—perhaps 0-1 wv/degree at 10°K. 


(6) It should be negative, in accord with the sign of the electronic 
charge. 


(c) For all metals it should be approximately proportional to 7’. 


The facts are quite otherwise (fig. 1 from the data of MacDonald e¢ al. 
(1958)): 


(a) Q is often several microvolts/degree. 


(b) It may be positive or negative, or even change sign as the tempera- 
ture varies. 


(c) The temperature variation may show large oscillations, and is quite 
different for each metal. 


That Q may be so large was explained by Hanna and Sondheimer 
(1957) in terms of ‘phonon drag’. This effect, now well understood in the 
case of semiconductors (Herring 1954) enhances the heat current associated 
with a stream of electrons, and hence increases the thermopower. But, 
whilst (a) was explained, the other anomalies remain. Bailyn (private 
communication) has taken a further step, by noting that the sign of the 
phonon drag is not invariable in principle. Consider, for example, an 
electron-phonon ‘Umklapprozesse’. (U-process.) An electron going 


+ Communicated by the Author, 
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from a state of wave-vector k, to state k’, with the emission of a phonon 
of wave-vector q, loses crystal momentum of amount 


K=k—k’=g+q, (nes ae) 


aeols) 


+2.0 


Q pv/°e 


Ge (Jabs) 


-1.0 


Absolute thermoelectric power of the alkali metals, from the data of 
MacDonald e¢ al. (1958), who give details of the specimens. 


where g is a reciprocal-lattice vector. It is quite easy, as we shall see, 
for q to be in the opposite direction to K. As one might say, the phonon 
is emitted “backwards ’ when the electron is scattered, and not necessarily, 
as in the normal ‘N-process’ in which g=0, in the forward direction. 
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A phonon heat-current is thus created in the direction opposite to the 
main stream of electrons (fig. 2). 

We still have to understand the temperature variation of Q, and the 
reason why it should differ from metal to metal In this paper some 
calculations are reported which seem to make this explicable, and which 
show that this phenomenon is a useful and sensitive index to the shape 
of the Fermi surface in the alkali metals. All that we do is to make a 
numerical calculation of the ‘lattice thermopower’, Q,, i.e. the contribution 
of phonon drag, taking account of both N-processes and U-processes 
using a typical form of the electron—phonon matrix element, and allowing 
for distortion of the Fermi surface towards the zone boundary. 


Fig. 2 


<J 
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The derivation of a formula for Q,, will be published elsewhere (Ziman, 
Electrons and Phonons, Oxford University Press, § 9.13). From the 
variational principle (Ziman 1956) 


ba aC re a - 
Sails cae) rae ns 


where 7a is the number of electrons per atom, and C, is the lattice specific 
heat. The other symbols signify scattering integrals, viz. 


Pr=| | fee) Aas dk dq ak aad seat (3) 


P= {{ [rer dk dq dk’+other terms. . . . (4) 


Here ©, is the steady state electron distribution, ®, is the phonon distri- 
bution, and the integrals contain the probability of a process in which 
an electron is scattered from k to k’, with the creation of a phonon. Only 
electron-phonon interactions are here dealt with explicitly, but there 
may be other terms in P,,, arising from phonon-phonon interactions, 
scattering of phonons by impurities, etc. 

To proceed, we must make an assumption about the various 
distribution functions. As usual, we put 


O,=k.u, O,=q.4u, Ge ee ee te ta (0) 
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where u is a unit vector in the current direction. For the phonons, in 
particular, this is not very accurate (cf. the case of a semiconductor, 
where the long wave lattice waves go much further from equilibrium), but 
is not too bad when scattering by electrons is the main source of resistance 
to lattice conduction. With these assumptions, our formula can be 
expressed very simply, at low temperatures : 
4 3 
a=- p= (5) iNaseak Faas git ale foo 
jJe|ma 5 \O (7?) 

The brackets < ) refer to averages over all electron-phonon scattering 
processes, at the rate at which they occur in the steady state. 

If there are only N-processes, g=0 in (1), making (3) and (4) equal. 
This is the case studied by Hanna and Sondheimer, who included in the 
denominator an empirical term as if the phonons were scattered by point 
impurities according to a Rayleigh law. By adjusting the amount of 
this scattering, they were able to fit the curve for Na up to quite high 
temperatures. There are many reasons, however, for supposing that 
U-processes are just as important. For example (Ziman 1954), in order to 
explain the lack of a minimum in the thermal conductivity of the alkali 
metals, we must replace the Bloch constant C? in the electron phonon 
scattering probability by a function @? (K) which is not zero when the 
length of the scattering vector k exceeds the radius, @ of the Debye sphere. 

The problem of averaging for U-processes is geometrically complicated 
but it is not difficult to find a reasonable approximation (loc. cit., § 9.5). 
With the aid of (1) and the geometrical identity 

Ke qs Ke Go) ee 
it can be shown that 
K.g)= [PMCD4 Coola) Me gay (OAC) Srna) 
(e*— 1)(1—e-*) (e?— 1)(1 —e-*) 
(8), (9) 
with z=hv/kT =q0/@T. Here the terms in @?(q) refer to N-processes ; 
if ©7(q) were constant, as in the Bloch theory, each integral would be of the 


familiar type (7/0)? (0/7). The other terms in each integrand derive 
' from U-processes, and are defined by 


a eS 

Cron(Q)= a }(K? +q2—g2)K°@"(K)dK; . . (10) 
299 g—-d 

Cry (Q)= oe Al K°@"K)dK: sn, 4) ba ea 


It is easy to recognize the average of (7) in the first of these. It is taken 
over that part of the Fermi sphere which can be reached by a phonon 
of wave-number g. The lower limit of K shows that it may not be less 
than g—q (when g and q lie in a straight line). Nor can K exceed 2k, the 
diameter of the Fermi sphere. For values of g less than 


Iy=9 — 2k Se ee eee LON 
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U-processes are impossible (see fig. 3), and the functions (10) and (11) 
are zero. The coordination number, z’ of the reciprocal lattice, counts 
the number of places on the Fermi surface where this lower limit on q is 
approached. 


Fig. 3 


Zone Bound ary 


Geometry of the calculation for a non-spherical Fermi surface. 


We see now that @7),,(q) tends to be negative. In (10), the term 
in —g?is much the largest. Because gis larger than 26, most contributions 
to <K .q) come from configurations where K and q are in opposite senses. 
Thus, N-processes and U-processes pull in opposite directions, and the 
sign of Q, will depend on their balance. This, in turn, depends critically 
on the form of the electron-phonon interaction, and on the shape of the 
Fermi surface. 

The best form of @?(K) for describing the transport properties of the 
alkali metals is due to Bardeen (1937). For the sake of simplicity we 
shall use the main variable factor in this formula, i.e. 


€2(K)= {sin (Kr,) —(Kr,) cos (Kr,)}#/(Kr,)§, . . . (13) 


which arises, basically, from the overlap integral of two plane waves 
differing in wave vector by K in Wigner—Seitz cell of radius r,, and which 
is characteristic of all theories which allow U-processes at all (Jones 1956). 
The other factors in the Bardeen formula arise from estimates of the 
deformation potential, and of the way it is screened by the free electrons. 
These factors are much less variable, and are not more significant in the 
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calculation than the corrections for phonon dispersion, and polarization, 
which we have already neglected. 

The formal analysis, up to now, is based on the assumption that the 
Fermi surface is nearly a sphere. Whilst this may be true topologically, 
we know that some distortion is probable in any real metal, and that an 
increase of only 14% in its radius would bring the surface into contact 
with the zone boundary. Even a small deviation in this direction can 
thus have an important effect on q,, the lower limit of U-processes. 

To allow for this possibility, we shall take & in (10 and (11) to be an 
adjustable parameter, capable of assuming larger values than &, say, 
which is what would be for an exact sphere containing one electron per 
atom. That is, instead of deforming the Fermi surface, we simply swell 
it until it is the same minimum distance from the zone boundary as it 
would have been if it had been deformed without change of volume (fig. 3). 
It is just at those crucial points near the zone boundaries that we want to 
get q, correct; precisely how the surface behaves in between does not 
much affect the calculation. We shall call 


e= (R= hs)/k, oo eee ee 


the measure of the anisotropy of the Fermi surface, that is, the amount 
by which it bulges, beyond its average radius, towards the faces of the 
Brillouin zone. 

Equations (6)-(13) now define the computation. The integrations 
were done by simple trapezoidal sums on a desk calculator, and no great 
accuracy is claimed for them. Curves showing Q, as a function of the 
reduced temperature, 7'/@, are plotted in fig. 4, for various values of e. 
For comparison, the experimental data of fig. 1 are plotted on the same 
scale. The reduction of the temperature requires some choice of values 
of ®. Ihave used the following: Li, 430°; Na, 160°; K, 100°; Rb, 68°; 
Cs, 54°, which seem representative ‘specific heat’ values of this fickle 
constant. This at least gets the scale of C, correct, even if it is not the 
best value to put into the scattering integrals P,, and P,,. Where 
there is a very big impurity effect, curves for different specimens are 
shown for some of the metals. 

The first point to be noted is the very precise fit, below @/10, between 
the experimental results for Na, and the theoretical curve for a spherical 
Fermi surface. There are no adjustable parameters in this, which surely 
shows that the theory is on the right lines. It is true that, at higher 
temperatures, the thermopower is not as large, in magnitude, as our 
formula predicts. Here the calculations of Hanna and Sondheimer 
would apply ; that is, the denominator of (2) is increased by the scattering 
of phonons from impurities. The case of specimen IV, a dilute alloy 
of K in Na, where the thermopower is not nearly so large, confirms this. 

The curious curve for K also follows extremely closely the shape 
calculated for e=7%. The difference could well be the standard electronic 
thermopower, 7°k?7'/e¢y, which would just prevent Q from changing 
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Absolute thermoelectric power as a function of reduced temperature. 
lattice thermopower calculated for various degrees of 
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anisotropy of the Fermi surface. 
calculated for N-processes alone. 
calculated for an arbitrary form of electron-phonon 


interaction. 
calculated by Hanna and Sondheimer (1957). 
experimental data from fig. 1. . 
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sign, as observed. The cause of the wiggle is clear. As the temperature 
falls, the term Cou? in (8) becomes dominant over @?, and Q,, becomes 
positive. This is because the typical phonon is only a bit above q, in 
wave number, and K and q are rather rigorously aligned into opposition. 
But then, at lower temperatures, these U-processes are frozen out, and 
the lattice thermopower reverts to its negative value, under the influence 
of N-processes alone. In principle, it should tend to the curve N, which 
is just the limiting form of (6) when both averages are equal. 

The heavy alkali metals, Rb and Cs, are not so well fitted. Qualitatively 
their thermopowers follow the pattern of the curve for «= 10%, Le. rise 
to a high positive peak, and then fall, perhaps to zero and negative values, 
at the lowest temperatures. It is not difficult to fabricate a form of 
electron—phonon interaction which would show such behaviour. For 
example, the curve labelled F (for Fudged) indicates what happens if we 
keep «=10%, but arbitrarily double @?(K) at large values of K. There 
is no mathematical reason why the peak should not be raised as high as 
we want by a suitable choice of @7(K). For Rb we note that the curves 
cut the axis, just as if U-processes were being frozen out at some value 
of q,, corresponding to « lying between 7 and 10%. The same may be 
true for Cs but the extrapolation is not reliable. In these elements it is 
interesting to note the effect of impurities. In Rb, an increase in the 
resistivity by a factor 9 divides the thermopower by about 2-5. Presumably 
this is due to an extra term in the denominator of (6), from the scattering 
of phonons by these impurities. 

Finally consider Li. Here the thermopower is decisively positive, 
and shows no sign of falling below the axis. Is this because, in Li, the 
Fermi surface touches the zone boundary? The calculation is not really 
adapted to this case, and the data are meagre, but it certainly 
seems So. 

To clinch the argument of this paper, we need independent evidence 
on the shape of the Fermi surface in the alkali metals. This has been 
discussed in detail by Cohen and Heine (1958) and Garcia-Moliner (1958). 
In each metal the exact amount of anisotropy is not yet certain, but there 
is a general pattern, into which the present results fit perfectly. Indeed, 
the lattice thermopower is probably one of the properties most sensitive 
to any anisotropy of the Fermi surface, and one of the best indices of its 
shape. 

It also seems as if one could deduce fresh information about the form 
of the electron-phonon interaction. If, as in Rb, one can roughly fix the 
shape of the Fermi.surface, one may then play with the form of @2(K) to 
make theory and experiment agree. It is not difficult to see how it could 
be improved in the right direction. Then, in turn, one could compute 
the other transport properties, and look for complete consistency. Until 
we have a sound derivation of the electron-phonon scattering at 
large angles, this is a promising empirical approach to this difficult 
problem. 
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In two earlier letters to this Journal (MacDonald ef al. 1958 a,b) we have 
summarized some results of note which we have obtained on the Group IB 
and alkali metals in the temperature range below 1°K. We wish here to 
mention some further results of interest on the heavier alkali metals and 
on various transition metals. We deal first with the latter. 

In our first note we reported that the IB metals, and gold in particular, 
showed remarkably large thermoelectric powers—greatly in excess of 
current theoretical predictions—and consequently we wished to examine 
other metals of neighbouring groups. Part of our interest stems from the 
fact that a metal with a sufficiently large thermoelectric power might 
enable thermoelectric refrigeration to be realized at low temperatures 
(cf. MacDonald et al. 1959). We have now made measurements on 
specimens of nickel, palladium, iron, cobalt and platinum; in general 
these metals show much lower absolute thermoelectric powers, S, than 
gold. Platinum of high purity, however, with S ~ — 1pv/°c at about 1°K 
approaches gold most nearly in magnitude of thermoelectric power; it is 
also notable that in platinum S is remarkably constant below say 1:5°K 
down to as low as about 0-3°K. An interesting feature is the dependence 
in sign and magnitude of thermoelectric power on impurity found in some 
of these metals even at these very low temperatures. Figure 1 illustrates 
this behaviour for two samples of nickel and some results for palladium 
and iron are also summarized in table 1. It will be seen that in these 


Table 1 
Metal Approximate residual Absolute 
pave resistance ratio thermoelectric force, E, 
(B4.9°%/ BR W295°K) at about 1-2;°K (volts) 

Nil IG s<¢ WO +6 x10-° 
Ni II 6-4 x 10-3 —2-9x 10-8 
Pd I 43 x 10-3 ~+6 x10-° 
Texel IU 36 x 10-3 —3-5x 10-8 
Fe I 7 se 1k +3 x10° 
Fe II 9-4 x 10-3 —2-2x10-8 
Pt I BSS MO —1:2x10-6 
Pt il 49 x 10-3 —3:3x107 
Co I MO se KF —2-5x 10-8 
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three metals the transition to the purer sample of metal leads to a change 
of sign of thermoelectric force from positive to negative together with a 
quite considerable increase in magnitude. Further experiments are 


planned to examine this process in more detail and to investigate the 
influence of specific impurities. 


+1x 1078 


Els 


ORIOL” 


=3 x 100" 


ABSOLUTE THERMOELECTRIC FORCE ,E (VOLTS) 


Absolute thermoelectric force observed on two specimens of nickel. 
A. Specimen Nil (see table 1); the full line represents the equation 
E=+4-2,~x10-* T? (volts), and thus S= +8-5 x 10-® T (v/°c). 
B. Specimen Nill (see table 1); the full line represents the equation 
E=-—1-8, x 10-8 T? (volts), and thus S= —3-6, x 10=2 7 (vir). 


A suggestion of Dr. C. V. Heer (Ohio State University) that the high 
entropy contribution found in cobalt below 1°K (Heer and Erickson 1957) 
might perhaps lead one to look for a large thermoelectric power, was 
followed up.° A sample of cobalt was measured, but a quite normal 
thermoelectric force was observed as indicated in table 1. 
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Table 2 
Metal Approximate residual Absolute 
: ee resistance ratio thermoelectric force, EL, 
aga es (Rayon I wenn) at about 1-2,°K (volts) 
Rb I 30 x 10-3 ee ea oa! Oe 
Rb IIL 3-5 x 10-3 peo Ome 
Cs I 6-5 x 10-3 Ao Lm 
Cs IT 1-4 x 10-3 alert lew seal (i) 2 
Fig. 2 
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Absolute thermoelectric force observed on two specimens of rubidium. 
A. Specimen RbI (see table 2); the full line represents the equation 
E= +8-4x 10-* T? (volts), and thus S= +16-8 x 10-9 7 (v/°o). 
B. Specimen Rb II (see table 2); the full line represents the equation 
E=—1-6x10-° T® (volts), and thus S= —4-8 x 10-® T? (v/°c). 
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The sensitive dependence on purity of the metals discussed above is algo. 
found in rubidium and the behaviour of two samples is shown in fig. 2 (see 
also table 2). It is interesting to notice again that in the purer sample the 
thermoelectric power, S, becomes negative and also that S then appears 
to become proportional to 7? (rather than to 7’) at the lowest temperatures 
(cf. also MacDonald et al. 1958b). In caesium, however (see table 2), the 
thermoelectric power appears to remain positive even in very pure samples. 
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In a recent paper Frank et al. (1958) present evidence in support of the idea. 
that most ‘trigons’ on natural octahedral faces of diamond are etch-pits. 
They suggest also that the trigons may be centred on dislocations. The 
diamond they examined, a good quality Type I ‘blue-white’ octahedral 
stone, showed clusters of trigons on each face. There was also a centre 
of strain inside the diamond which could be detected by its birefringence. 
The strain nodulus was located approximately at the intersection of 
perpendiculars from the trigon clusters on each face. Hence it appeared 
likely that there were dislocations radiating from the strain nodulus 
outward to the trigons. An X-ray examination of the diamond has been 
undertaken to test the validity of this proposed dislocation distribution. 

It has been shown that individual dislocations in crystals may be 
detected by x-ray diffraction (Lang 1958). A general survey of the 
dislocation content of a nearly perfect crystal can be made using the 
‘projection topograph’ technique (Lang 1959). A pair of projection 
topographs, derived from Bragg reflection in a given plane and its inverse,. 
provide a stereo-pair which enables the three-dimensional distribution of 
dislocations to be seen. Such a stereo-pair is reproduced in fig. 1 (a) 


pee re eee 
+ On leave at H. H. Wills Physics Laboratory, University of Bristol. 
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and (6) (PI. 56). The glancing angle for the 111 reflections used was small 
(eight degrees for Ag Ka radiation) so that the octahedron appears nearly 
as a parallelogram in projection. Both individual dislocations and close 
bundles of dislocations can be seen radiating from the strained inner region 
outwards towards the faces. Elastic strain and the high dislocation 
density prevent resolution of individuals close to this centre. Surface 
damage, especially along the edges of the crystal, accounts for specks of 
enhanced diffracted intensity. About one-half of all dislocations present 
should be visible in this stereo-pair of photographs assuming that the 
various Burgers vectors occur with equal frequency. The hypothesis of 
Frank et al. concerning the general distribution of dislocations is thus 
strikingly confirmed. It is planned to perform a detailed x-ray diffraction 
study of this crystal to establish the direction of Burgers vectors of the 
dislocations and to correlate their outcrop with individual surface features. 
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Probability Density of Electron Separation in a 
Uniform Electron Gas 


By N. H. Marcu and W. H. Youne 
Department of Physics, The University, Sheffield 


[Received November 7, 1958] 


We briefly report here a new method of approach to the correlation 
problem in a uniform electron gas, within the framework of density 
matrix theory. The earlier important paper of Mayer (1955) using this 
method has been shown independently by Tredgold (1957) and Koppe 
(1957) to lead to a variational scheme which is not complete, due to 
difficulties associated with subsidiary conditions. We propose here a 
uniform method of describing the electron correlations throughout the 
entire density range, and show that our method does indeed lead to 
accurate results in the high and low density limits. We hope to report 
detailed calculations for intermediate densities later, but at least we show 
in this note how a reasonably accurate picture of the correlation effects may 
‘be obtained over a wide density range. 
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Following Mayer (1955), we introduce the probability density F',(€,) of 
separation ¢,, between electrons; this being essentially the diagonal 
element of the second order density matrix. We use the variable &, 


rather than r, defined by 
3 \1/3 ; 
rss = (=) Cem aria cs ne) 


where 7s is the radius of a sphere containing one electron. Then, as 
Mayer has shown, the potential energy per particle V may be written 


l 3 1/3 ioe) 
v--3(2) ={ 4n€{1—F(E, rs)}dE 6 Sw. (2) 


s ) 
and to emphasize the importance of /’, we first prove a new result; that 
the total energy per particle # may be written exactly in terms of F,. 
This may be achieved by using the many-body virial theorem for this 
problem in the form given recently by one of us (March 1958). Re-writing 
this theorem as a differential equation for the energy H(rs), given V(rs) 
from (2), it may be shown that 


3 (/97\2/8 1 i) (es 
a ( _— + — V drs, ee Pee OB 
mitea) teal, Le - (3) 


the constant of integration being simply the Fermi energy. 

Thus we focus attention on the calculation of #,, and for simplicity of 
presentation we shall at first neglect spin. With WN particles, let 
W(r,, rg..-."y) be any antisymmetric wave function such that it is 
normalized to unity. We now define corresponding density matrices 
(see for example Léwdin 1955) 


y(ry'|r1) =v | er lola) uC heen lw Oro er a4) 


We il i / \ 
CT Cage A eel Are recital a) erate Oils cane ae 
7 (5) 
Then considering the function ® defined by 


O(r,, rota f (G5) Me re) eee .aeeanr(G) 


it is easily shown that 
[ore Feet Or fo. 2h) Ora. ..Ory—=l. 72. . (7) 


Furthermore ® is antisymmetrical in ry, r3...ry, and can therefore be 
used to construct density matrices of various orders for N —1 particles, 
in coordinates r,...ry. In particular 


—— P(ryre’[ryre) 
y(rylrs) 
is a first order density matrix for V — 1 particles, for any arbitrary fixed rj. 


P.M, 2B 
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We now apply the above general result to the electron gas problem. 
Here y(r,|r,), the mean electron density, is a constant, y say, and the 
‘particle density’ for N —1 particles is 


y 
x P(ryra|ryrs) = — P (rye) 


SS LS 


for fixed r,;. We may write in this case 
2 
oe Mata) =f(tye)s Me=he- PENA E SER) 


and it is sufficient to investigate (2/y)I(0 r.)=f(r.), where we can think of 
axes as fixed at electron 1. Clearly to assume any first order density 
matrix for N —1 particles is not sufficient, for the antisymmetry property 


requires that f(0)=0 (9) 


for particles of like spin. Then f(r) is intimately related to F,(€), the 
latter representing an angular average of f(r). 

Having seen how the second order density matrix may be re-expressed in 
a form related toa first order matrix for N — 1 particles we shall now attempt 
to set up a model in keeping with the requirements put forward here. We 
restrict ourselves to a first order density matrix of the simplest kind, 
constructed from a set of orthonormal single particle functions ¥,(r), 


as N-1 
2 b*(e' )b(r). 


We were essentially helped in our considerations by seeking a unified 
treatment which includes both the Hartree-Fock solution, known to 
be correct in the high density limit, and the low density case (Wigner 
1938). Wigner’s low density calculation assumed that the potential 
energy was minimized when the electrons went on to the sites of a 
body-centred cubic lattice. Now, of course, this can only be interpreted 
as a lattice relative to a given electron we have singled out, as no 
real set of lattice points can occur preferentially in a uniform gas. Thus 
the lattice can be thought of as describing the configuration of all the other 
electrons, once axes are fixed on one particle. Of course, even then, all 
orientations of this lattice are equally probable, and so eventually we must 
average over angles to obtain physically significant results. Wigner 
finds that orbitals centred on the lattice sites can be represented by 
harmonic oscillator functions 


a\ees.) for 1 
= (2) exp ( ; i) aie eae LO) 


for which 7? tends to the delta function in the low density limit. If the 
density is very small, the functions (10) can be treated as orthogonal 
(effectively zero outside their own sphere) and hence are a valid set from 
which to form a first order density matrix of the kind we require. How- 
ever, with NV particles, we have N harmonic oscillator functions, and to 
obtain a suitable first order density for N—1 particles we must remove 
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the density contributed by one of the orbitals ; the one centred on the origin. 
We can then build up an acceptable furineion IF, by averaging over angles. 
The result thus obtained for rg= 100 is eae in curve | of the figure, 
and for such a low density the present method should provide a truly 
accurate representation of the real situation. We see that the probability 
density F’, is violently oscillatory and of a highly complex functional form. 
In order to follow the change in /’, as rs is reduced, we have plotted in 
curve 2 the case rs=4. We should emphasize that extrapolation to such 
a density is not justified fundamentally, as the orbitals are now no longer 


fe) 1-0 2-0 3.0 an v3 4.0 
ce 


Probability density F,(é) as function of electron separation €. Curve 1. Low 
density form for rg=100. Curve 2. Extrapolated low density form 
for rs=4. Curve 3. Fermi hole, correct in the limit rg+0. 


orthogonal, and also do not make f(0) zero for particles of like spin. 
Various estimates of the range of validity of the low density argument of 
Wigner give rg> ~ 10-20. Nevertheless, whichever value we assume for 
the lower limit of approximate validity, it is clear that FP, is still oscillatory 
about the line /’, = 1, to a significant degree, and any acceptable probability 
density must contain such behaviour. 
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We see then that Wigner’s arguments can be fitted into the present 
scheme, and give a good deal of information about the form of 7. How- 
ever, we wish to see how the high density form can be obtained from the 
present method. One thing which suggests itself is to extrapolate the 
low density F,, into rg=0, but doing this we find no Fermi hole, due essen- 
tially to the failure of the Gaussian orbitals to (a) satisfy orthogonality 
requirements and (b) have nodes at the positions of like spin electrons. 
The question then arises as to what the localized functions centred on 
lattice sites must be in this limit. It seemed clear to us that since plane 
waves provide the correct solution in the high density limit, free-electron 
Wannier functions must be intimately related to this case. We conclude 
by showing that this is indeed so, and we demonstrate that for a one- 
dimensional Fermi gas, the exact solution may be obtained in this way. 

If we consider a lattice spacing z/k , then the free electron Wannier 
functions may be written 


1 sink(x—mafky) (—1)"sinkgx 


ee ee ee ee ee eee ih 
PW O= EYE (e—malky) — Grk)™@—mnlk) 
Summing over the sites of an infinite lattice we can write 
2 j 1 (° sink,(x’ —m/ko) sin ko(x—m/ko) 
* gs 0 et Se 
2 bn (Wnt) = i _« (@'—mr|k)  (w—mr]k) 
_ ky sin ko(x' — 2) (12) 
a ky(a' —x) 
Hence, in particular S Ee ee ay 


This latter result may be obtained by direct summation and is exact. 

Now, according to our model, we form the density for (N — 1) particles 
by removing the density contributed by the localized orbital centred on 
the origin, in order to make f(0)=0, and ensure normalization to (N —1). 
Then we find 


% HaMe) —poXa) = “| 1 FEE) Pelee yi 


In this case the exact second order density matrix can be written down 
from a two-by-two determinant (see for example Lowdin 1955) involving 
the first order matrix 

ky sin ko(a’ — 2) 


aw Kko(x’—x) ? 
the result being 
ky sin ko(a’ — a) 


1 1 (a — x) LN fsin k (x' —x)\? 

cS ee il eccrine) 

2| sin ko(a’ — 2) ko le) L | ko(x’ —2) at a 
T 


(a — x) 
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From (8), we must multiply (14) by y/2=ko/2m in order to compare with 
(15). The equivalence then evident proves that we obtain an exact 
description of the Fermi hole in this case, using localized orbitals on a 
lattice; the orbitals being free electron Wannier functions. 

So far, we have not seen how to extend our high density arguments 
exactly to the real three-dimensional gas, with two spin directions. How- 
ever, the following model will lead, at least, to a fairly accurate represen- 
tation of the Fermi hole. In this high density limit we think of Wigner’s 
body-centred-cubie lattice as two interpenetrating simple cubic lattices, 
and we centre on each lattice site a free electron Wannier function for a 
simple cubic lattice, the functions on one lattice being multiplied by a 
spin function a, those on the other lattice by a spin function 8. Then we 
can again form a first order matrix for NV — 1 particles, and the probability 
density we obtain may be written 


. . * 2 
ee! E c€, sinc€, sin cé, ee eee) 
DSS are et 


To obtain F',(€) we must average this over angles, and we have not yet 
achieved any closed result. However, we shall obtain a fair approxi- 
mation to the Fermi hole, as the following argument shows. If we go 
back to the basic form for the simple cubic Wannier functions, 


1 F 
b(r) = 8(aky)®” oe “a exp (vk . r)dk., dk, dk, ° ° (17) 


and replace the cube of volume 8,3 in k space by an equal volume sphere 
before performing the integration, then we obtain the exact Fermi hole. 

Further work is now in progress to determine an acceptable form of 
localized orbital ys(r) in the intermediate density range, which goes into 
the Gaussian function (10) at low densities, and into functions at least 
intimately related to the Wannier functions we have described in the high 
density limit. 
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REVIEWS OF BOOKS 


Solid State Physical Electronics. By ALDERT VAN DER ZIEL. (London : 
Macmillan, 1958.) [Pp. 604.] 37s. 6d. , 


Amona the rapidly growing number of text-books on solid state electronics 
this book occupies a unique place for several reasons. In the first place, the 
value for money is outstanding by present standards, there being just over 
600 pages. The emphasis throughout is on the physical principles involved in 
the working of various solid state devices, rather than on the more academic 
theory of the solid state or on the more engineering details such as circuit 
design. The range of topics is wider than might be expected from the title, 
and includes such things as thermionics, field, secondary and photo-electric 
emission in as much detail as might be expected in a book on vacuum 
electronics. 

With such a wide range, the treatment is necessarily somewhat condensed, 
but most of the important mathematics is worked through in detail and one’s 
feeling after reading a chapter is that the subject has been given a credible 
outline and the minimum necessary has been taken on trust. The number of 
references has been kept down by including among them a large proportion of 
text-books and review articles, but in spite of this the book seems well abreast 
of the times. A number of problems is included, many of a numerical nature 
with values quoted which are typical of current practice. 

Altogether the book is well produced, very readable, and can be strongly 
recommended. 1D; K2G, 


Low Temperature Physics and Chemistry. Editor: J. R. Dmtincer. (The 
University of Wisconsin Press, 1958.) [Pp. 676.] $6.00. 


Dvurine one week in August 1957, 225 papers were presented to 440 delegates 
attending a conference on Low Temperature Physics and Chemistry at Madison, 
Wisconsin. The detailed report now published provides a pretty complete 
picture, albeit somewhat pointillistic, of the situation in almost every branch of 
cryogenic research at the time of the conference. The editor has performed 
his task very thoroughly and produced a volume which will be of value to workers 
in this field for some time to come. It is perhaps a pity we have had to wait 
over a year for it, but that is the price one must pay for elegance of production. 
The monetary price seems very reasonable by modern standards. 


A. B.P. 


Expansion Machines for Low Temperature Processes. By S. C. Couurys and 
R. L. Cannapay. (Oxford: Clarendon Press, 1958.) [Pp.115.] 12s. 6d. 


Tuts short, moderately priced book brings together a large mass of information 
on the design and use of expansion machines. The evolution of the reciprocating 
expansion engine is traced historically from the air expansion machines used for 
refrigerating ships in the latter part of the nineteenth century to the first 
applications to air liquefaction in Claude’s machine. The further development 
of these engines for the liquefaction of hydrogen and helium is then described in 
more detail. This section contains much interesting information about 
Professor Collins’ own important contributions to machine design, and might well 
have been expanded further at the expense of some of the earlier chapters of 
the book. In particular, information about the modifications which have been 
introduced into commercial helium liquefiers and their effect on efficiency and 
maintenance would have been very valuable. ‘ 
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The latter part of the book deals with the design and performance of various 
turbine expansion machines used in large scale air liquefaction, and contains a 
useful comparison of the advantages of different types. If the book is intended 
primarily for low temperature engineers, this section might profitably have 
included a discussion of the changes in the rest of the plant which have to be 
made when a turbine is used instead of a high pressure reciprocating machine. 

The book is well illustrated and has a useful bibliography. 


Principles of Modern Physics. By A. P. Frencu. (New York and London : 
John Wiley & Sons.) [Pp. 355.] 54s. 


MopERN physics spreads its net so wide that it can mean whatever an author 
wants it to mean, and especially the field in which he is himself active. By 
the same token it involves so many trees that it is very easy to get lost among 
them, wandering on with a catalogue of detail. Dr. French knows the subject, 
and the dangers, very well and sets out to carve a clear path from one landmark 
to another. For him modern physics is the discovery of the individual atom 
and the exploration of its structure, by the classical, quantum theoretical and 
finally wave mechanical methods. He is concerned to trace the development 
of the essential physical concepts, rather than to impart a lot of information, 
with the advanced undergraduate or early postgraduate student in mind. 
The jacket (not the Preface) speaks of the ‘thorough examination of a 
limited number of topics’. On the contrary, the text is much more a concise 
introduction to the whole range of subject matter of modern atomic and 
nuclear physics. As the Preface says, ‘ collateral reading ’ is expected, and a 
considerable amount of it (supplemented by seminar discussions) will certainly 
be needed from the average student. The book has the clarity of Oldenberg’s 
Introduction to Atomic Physics at the level of treatment of Born’s Atomic 
Physics, to which it is a worthy successor. More one cannot expect at the 
length and the price. Vien. 


Elektronenbeugung. By Ernst Bauer. (Munich: Verlag Moderne Industrie.) 
[Pp. 240, with 120 figures.] DM 32. 
ELEcTRON diffraction has long since reached the stage of being accepted as a 
research tool with a field of application all its own. Techniques are well 
established and a great number of applications are to be found in the 
literature, but there has been a conspicuous lack of a practical handbook 
explaining the experimental procedures in a way that a non-specialist can 
understand. It is this need which the author sets out to meet, with a 
considerable degree of success. The book is divided into three parts : 
theoretical (56 pp.), description of apparatus and methods (95 pp.) and an. 
outline of applications (47 pp.). He deliberately excludes proofs, quoting 
only the essential relations, and hence is able to cover the geometric theory 
(Bragg) in 15 pp., the kinematic theory (Laue) in 9 pp. and the dynamic theory 
in 16 pp. The treatment of practical topics is similarly compressed, and it is 
doubtful whether a beginner will gain sufficient grasp of the subject unless he 
follows up the suggestions given for further reading. These are, however, very 
full, and 550 references are quoted in all. The survey of applications ranges 
over the chemical, electrical and optical industries, as well as metallurgy. 
The book is likely to find a welcome where it is particularly addressed : the 
industrial scientist who wants to know how’ to use a new technique and is 
more concerned about the how than the why of it. He will especially 
appreciate the list of suppliers of various types of equipment on the last page. 
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Handbook of Physics. Edited by E. U. Conpon and H. OpisHaw. (New 
York: McGraw-Hill Book Company, Inc., 1958.) [Pp. 1473.] 194s. 


Tuts massive volume with short articles by experts on every subject from 
arithmetic to nuclear fission could be very useful to anyone planning an honours 
course in physics, a series of text books or wanting to know the salient factor 
on some field of research outside his speciality. Some of the articles are very 
up to date with references to papers as recent as 1957. N.F.M. 


Progress in Orystal Physics, Vol. 1. By R. S. Krisanan. (Madras : 

S. Viswanathan, 1958.) [Pp. 198.] Rs 20 or 30s. 

Tuts book is both a review of experimental and theoretical work in certain 
selected fields of Crystal Physics and also a compilation of the available data 
concerning these properties. After a short introduction which refers to 
geometrical crystallography and to crystal growth, there are seven chapters 
dealing respectively with Thermal Expansion, Thermal Conductivity, Elastic 
Constants, Photoelastic Effect, Thermo-optic Behaviour, Faraday Effect in 
diamagnetic crystals and dielectric properties of ionic crystals. Each chapter 
has been written by two or three members of the Department of Physics of the 
Indian Institute of Science, and altogether there are seven contributors. 
Professor R. 8S. Krishnan has contributed to four chapters himself and is to be 
congratulated as general editor on the coherence and readability of the work. 
The summaries of the experimental methods used in the measurements of the 
physical properties and the theories by which they can be interpreted are 
clearly described. In many cases the reader is referred to a full review or 
book and little is given here. A valuable feature of the articles is that each 
author has done original work in the subject on which he has written. The 
reviews cover the literature and give a fair account of conflicting theories 
and experimental results. 

The difficulty of finding adequate tables of data on the physical properties 
of crystals is well known. The present book is therefore sure to be welcomed 
for the full, clear, and detailed tables provided at the end of each chapter. 
The collections of references to original papers and books will also greatly 
facilitate the study of the subject by research workers. 

The choice of symbols for representing the many physical quantities discussed 
in such a book presents many problems to the author as well as to the reader. 
Only in one chapter is there a glossary of symbols. In any future edition it 
would be helpful if a general index of symbols were provided. The early 
appearance of succeeding volumes will be much appreciated by all those who 
study the physical properties of crystals. W. A. W. 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents. | 
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Fig. 2 


A 1000 A layer of silver deposited at A 1000 4 layer of silver deposited at 
a substrate temperature of 40°c. a substrate temperature of 75°c. 
Mica [010] azimuth. Mica [010] azimuth. 

Fig. 4 Fig. 6 


A 1000 & layer of silver deposited at a A 504 layer of silver deposited at a 
substrate temperature of 275°c. substrate temperature of 275°c. 
Mica [010] azimuth. The elongation Mica [010] azimuth. 


of the diffraction spots arises from 
lack of electron penetration into 
the smooth (111) surface. 
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Fig. 7 


< Angle of scattering 
approx. 125° 


< Shadow edge 


< Central spot 


A high-angle Kikuchi pattern from a 1000 & layer of silver deposited at a sub- 
strate temperature of 270°c. Glancing angle of incidence approximately 
ten degrees. Mica [100] azimuth. 
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Fig. 1 Fig. 4 


Prominent extinction contours attri- A (100) film of gold showing stacking 
buted to the occurrence of ‘dimpled faults and individual dislocations, 
regions’ in the film (see fig. 2 (0)). such as at A. (x140 000.) 


(x 57 000.) 


Fig. 5 Fig. 6 


Dislocation loops in a (111) gold film. The appearance of stacking faults during 


( x 340 000.) the electron bombardment of the 
specimen. (x 120 000.) 
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D. W. PASHLEY Phil. Mag. Ser. 8, Vol. 4, Pl. 46. 
Fig. 8 


Fig. 9 


Stacking faults in a (111) film as seen with A stacking fault in a (111) film 
the electron beam inclined to the showing fringes with a twisted 
plane of the film. (x 180 000.) appearance. (x 180 000.) 

Fig. 10 


A low angle boundary in an annealed gold A regular dislocation array in an annealed 
film. (x71 000.) gold film. (x71 000.) 
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